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Zusammenfassung 



Das Problem wechselwirkender Elektronen in zwei Dimensionen ist spatestens seit der Entde"ckung 
dcr Hoch-Tcmpcrat\ir-Supralcitung von grosser Aktualitat und Gegcnstand intcnsivcr Forschnngs- 
tatigkeit. In der vorliegenden Arbeit wird ein erweitertes Hubbard-Modell auf einem Quadrat- 
gitter untersucht. Obwohl man weiss, dass die Leitungselektronen in den Kupraten relativ stark 
mitcinandcr wcchsclwirkcn, bcschrankc ich mich hier bewusst auf den Grenzfall einer schwachen 
Wechselwirkung zwischen den Elektronen. 

Selbst eine beliebig schwache Wechselwirkung fiihrt bei geniigend tiefer Temperatur zur Bildung 
von starken Korrclationcn und scliliesslich zum Phascniibcrgang in cincn gcordneten Zustand. Die 
Natur dieser Instabilitat hangt stark von den Eigenschaften der Fermi-Flache ab. Wenn keine beson- 
dere Situation vorliegt, ist es ein Ubergang in einen supraleitenden Zustand. Dies obwohl die Wech- 
selwirkung der Elektronen rein repulsiv ist und keine phononischen Freiheitsgrade berii" cksichtigt 
werden. Eine mogliche Interpretation dieses Resultats ist, dass der Austausch antiferromagnetischer 
Spinfluktuationen auf eine effektive Anziehung zwischen den Quasiteilchen fiihrt. 

Im Spezialfall halbcr Bandfiillung crfiillt die Fermi-Flache die "perfect nesting" -Eigcnschaft. 
Gleichzeitig hat die Dispersionsrelation der Elektronen einen Sattelpunkt am Fermi-Niveau, was 
zu einer Van-Hove-Singularitat in der Zustandsdichte fiihrt. Diese beiden Eigenschaften bewirken, 
dass es zu einem Wcttbewcrb zwischen sechs verschiedenen Instabilitatcn kommt. Neben s- und 
d-Wellen Supraleitung kann eine Spindichtewelle oder Ladungsdichtewelle auftreten. Die zwei weit- 
eren Instabilitaten sind sogenannte Flussphasen, die durch spontan zirkulierende Ladungsstrome, 
respektive Spinstrome charakterisiert sind. 

Um die konkurrierenden Instabilitaten vorurteilslos gegeneinander abzuwagen, benotigt man 
die Technik der Renormierungsgruppe. Der notige Formalismus wird hier auf einem moglichst 
elementaren Niveau prasentiert. Die Idee ist dabci, ausgehcnd von der naiven Storungsrechnung 
die dominanten Terme (d.h. die Klasse der Parkett-Diagramme) konsistent aufzusummieren. Das 
modernere Konzept der Wilson'schen effektiven Wirkung wird ebenfalls kurz diskutiert und zu dem 
hier verwendeten Ansatz in Bezug gebracht. 

Als Resultat erhalte ich ein Phasendiagramm als Funktion der Modell- Wechselwirkung, die 
neben dem iiblichen Hubbard- Term auch eine Dichte-Dichte- und Spin-Spin- Wechselwirkung zwis- 
chen benachbarten Gitterplatzen enthalt. Fiir das repulsive Hubbard-Modell erhalt man wie er- 
wartet eine Spindichtewelle. Dieses Phasendiagram ist im Grenzfall unendlich schwacher Wechsel- 
wirkung vermutlich exakt, denn auf der Grenzlinie zwischen zwei benachbarten Phasen weist das 
System dann jeweils eine besondere Symmetric auf, die die Entartung der beiden benachbarten 
Phasen garantiert. Das physikalische Bild dieser Instabilitaten wird erganzt durch Grossen wie die 
uniforme Spin-Suszeptibilitat und die Ladungs-Kompressibilitat, deren Verhalten bei tiefen Tem- 
peraturen qualitativ vorausgesagt wird. Dabei zeigt sich eine allgemeine Tcndenz zur spontanen 
Verformung der Fermi-Fache (Pomeranchuk-Instabilitat), die jedoch nicht notwendigerweise mit 
einer Zerstorung der "perfect nesting" -Eigenschaft einhergeht. 

Obwohl sich die konkurrierenden Supraleitungs-, Spin- und Ladungskorrelationen gegenseitig 
beeinflussen, kann sich die fiihrende Instabilitat letztlich als einzige ungehindert entwi" ckeln. Dieses 
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ZUSAMMENFASSUNG 



Resultat ist im Widcrspruch zu altcrcn Losungsvorschlagcn, wo die Spindichtewelle jeweils zusam- 
men mit der d-Wellen Supraleitung und der Ladungs-Flussphase auftrat. 

Abschliessend stellt man fest, dass sich zwei wesentliche Eigenschaften der supraleitenden Kuprate 
im schwach wcchsclwirkcndcn zwcidimcnsionalcn Hubbard- Modcll wicdcrfindcn. Namlich cine anti- 
ferromagnetisch geordnete Phase bei halber Bandfiillung und das Auftreten von d-Wellen Supraleitung 
im dotierten Material. Leider ergibt sich aus der vorliegenden Theorie, die von schwacher Wechsel- 
wirkung ausgcht, kcinc Erklaning der ausscrgcwohnlichcn Eigensehaften in der nictallischcn Phase 
des optimal dotierten und unterdotierten Bereichs. Starke Wechselwirkung scheint daher ein unab- 
dingbarer Bestandteil fiir eine adaquate Theorie dieser Materialien zu sein. 

Trotzdem konnte sich cine konscquentc Wcitcrfiihrung dcs hicr vcrwcndetcn storungsthcorctischen 
Ansatzes als durchaus fruchtbar erweisen, wobei unter anderem die Selbstenergie der Elektronen 
wichtig wird. Die dazu notwendige Technik ist jedoch noch nicht ausreichend entwickelt. 



Resume 



Les systemes electroniques bidimensionnels sont d'une grande actualite tout particulierement depuis 

la dccouvcrtc dc la supraconductivitc a haute temperature. Ici, on se restreint a I'ctude d'un 
modele de Hubbard etendu, a la limits d'un couplage faible. En general, le gaz electronique subit 
una instabilite supraconductrice meme sans phonons. Cependant, dans le cas special d'une bande 
demi-rcmplie, la surface de Fermi est emboitcc ct se troiive a une singularitc dc Van Hove. Cctte 
situation conduit a une competition entre six instabilites differentes. Outre la supraconductivite 
en onde s et d, on trouve des ondes de densites de spin et de charge ainsi que deux phases qui 
sont caracterisces par des courants circulaires dc charge et de spin respectivemcnt. Le formalisme 
du groupe de renormalisation est presente en reliant I'idee de la "sommation parquet" au concept 
plus moderne de Taction effective de Wilson. Comme resultat on obtient un diagramme de phases 
richc en fonction de I'intcraction du modele. Cc diagramme dc phase est exact dans la limite 
d'une interaction infiniment faible, puisque dans ce cas les lignes de transitions sont fixees par des 
symetries du modele. Les comportements a basse temperature de la susceptibilite de spin ainsi 
que de la comprcssibilitc dc charge complctent I'image physique dc ces instabilites. II s'avcre que 
la surface de Fermi a une tendence generale de se deformer spontanement, mais I'emboitement 
n'est pas detruit. En resume, le modele de Hubbard a couplage faible reproduit deux proprietes 
essentielles des cuprates: une phase antiferromagnetique a demi remplissage et la supraconductivite 
en onde d dans le cas dope. Mais elle n'explique pas les proprietes inhabituelles de I'etat metallique 
dans le regime sous-dope. Une extension systematique de I'approche perturbative pourrait aider a 
mieux comprendre ces proprietes, mais reste difficile puisque les techniques necessaires ne sont pas 
encore completement developpees. 
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Abstract 



Interacting electrons in two dimensions are of particular interest in relation to high-temperature 

superconductivity. In this thesis, I study a two-dimensional extended Hubbard model in the weak 
coupling limit. Quite generally, the electron gas is unstable towards a superconducting state even in 
the absence of phonons. However in the special case of a half-filled band, the Fermi surface is nested 
and the system is at a Van Hove singularity. In this situation, there arc six competing instabilities: 
s- and d-wave superconductivity, spin-and charge-density waves and two phases with circulating 
charge and spin currents, respectively. The required renormalization group formalism is presented 
on a most elementary level, connecting the idea of the "parquet summation" to the more modern 
concept of Wilson's effective action. As a result, a rich phase diagram is obtained as a function of 
the model interaction. This phase diagram is exact in the weak coupling limit, since the transition 
line between two neighboring phases is then fixed by symmetries. The physical picture of each 
instability is completed by studying the low temperature behavior of the spin susceptibility and the 
charge compressibility. We also observe a general trend towards a Fermi surface distortion, but the 
nesting is not destroyed. In summary, the weak-coupling theory of the Hubbard model reproduces 
two essential features of the cuprates, namely an antiferromagnetic phase at half-filling and d-wave 
superconductivity in the doped material. But it does not explain the unusual properties of the 
metallic state in the underdoped regime. A consequent extension of the perturbative approach to 
sub-leading orders which would imply self-energy corrections could reveal further insight, but the 
required techniques are not yet fully developed. 
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Chapter 1 

Introduction 



The interest in the problem of interacting electrons in two space dimensions has been increased 
tremendously by the discovery of quasi two-dimensional (2D) materials with novel electronic prop- 
erties. The high- Tc superconducting cuprates are certainly the most famous examples^]. Other un- 
conventional superconductors with a quasi 2D electronic structure are strontium ruthenatc ||MRS01 |, 
some graphite-based compounds ||dSTKOl[| , |N+Ol[ , and doped Ceo crystals jSKBOOj PKBOlj . 

Yet another class of quasi-2D materials is given by some transition metal dichalcogenides with inter- 
esting charge-density wave formations |APBL01|. 2D electronic systems have also been experimen- 
tally realized in Si- or GaAs- based semiconductor devices, where the 2D electron gas undergoes a 
metal-insulator transition |K+95, H+98|. 

The metallic phase of the cuprates deviates considerably from the predictions of Landau's Fermi 
liquid theory, i.e. the theory of usual metals based on the concept of quasi-particleg|^. Impressing 
evidence for the breakdown of the quasi-particle concept in a cuprate material has been provided 



very recently by the violation of the Wiedemann- Franz law [H+01 



Most of the unusual properties of the cuprates are likely to be linked to the quasi- two-dimensional 
nature of their electronic structure close to the Fermi energy. Therefore certain single-band 2D mod- 
els of interacting electrons may be able, in principle, to account for at least part of the anomalies 
observed in these compounds [ And88| , |ZR88 |. Angle-resolved photo-emission spectroscopy mea- 
surements support this idea, since they show a single well-defined Fermi surface at least in the 
over-doped regime [T+88, 0+8£, G~'~90| p|. Furthermore the undoped insulating parent compounds 
are well described by the 2D spin 1/2- Heisenberg model [ Man91 |, which is also a single-band model. 

Unfortunately, even very simple models, such as the 2D Hubbard or the 2D t — J model, have 
so far resisted a rigorous analysis. Moreover, the available numerical studies are not yet conclusive 
enough for making definite predictions for, e.g., the zero-temperature phase diagram of these many- 
electron systems. Another difficulty is that fiuctuations (both thermal and quantum) are strong in 
two dimensions so that mean-field approximations cannot be trusted. 

The coupling between electrons in the cuprates, for instance the parameter U of the Hubbard 
model, is large, i.e. of the order of the bandwidth. Therefore it is not clear whether a ground state 
consisting of occupied Bloch orbitals with energies below e^? is a good starting point or whether 
one has rather to think in terms of configurations of singly occupied and empty sites (doped Mott 
insulator). Actually, the successful analysis of the insulating phase in terms of the Heisenberg model 



suggests that the Mott insulator is the appropriate reference state [Man91 



I deliberately choose the limit of weak bare couplings, keeping in mind that this parameter range 

bee | |Wal96t for a. ri iritrn riiint.inn to the cuprates and superconductivity in general and [pMOCf for a recent review 
^See tor example [FNGG 

o — ' 



''For recent reviews on photocmission spectroscopy of the cuprates, see [JFVOl, DLSOl, 3+01 
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CHAPTER 1. INTRODUCTION 



may miss completely some important characteristic aspects of the region of strong bare interactions. 
For example, the spin-density wave (SDW) instability of the half-filled Hubbard model at weak 
coupling, which will be studied in Chapter ^, fulfills the physical picture of a Slater instability 
Sla51| rather than a Mott transition. The Slater instability arises through the formation of a 
SDW, whereas the Mott transition is characterized by the formation of local magnetic moments, 
which are not necessarily ordered^. Nevertheless, it is to my opinion important to understand the 
weak coupling limit, which is already a non-trivial problem. It cannot be excluded that certain 
properties are qualitatively the same over the whole range of couplings, as is the case for the ID 
Hubbard mo del, a Luttinger liquid for all positive values of U and all densities except n ~ 1 
| S092| , |Voi95tl . 

Although it is not excluded that lattice vibrations (phonons) are important to understand high- 
temperature superconductivity, I will completely neglect them in this thesis. 

The most clear picture of two-dimensional interacting electrons has been obtained for the 2D 
jellium model with its circular Fermi surface. A series of rigorous studies has shown that the Landau 
Fermi liquid theory is stable at not too low temperatures such that |J71ogT| < const, where U is 
the strength of the local interaction |Sal98 , DROO , but it breaks down at a critical temperature 
Tc ^ exp — '^'j'^^^ . For lower temperatures, the properties of the interacting system are no longer 
analytically connected to those of the non-interacting system. The usual interpretation is a phase 
transition into a superconducting state, although there is no rigorous statement about T < Tc 
available up to now. 

Strictly speaking, there can be no long-range order in two dimensions due to thermal fluctu- 
ations |MW66|. Tc is thus interpreted as a mean-field critical temperature, which is close to the 
temperature of the Kosterlitz-Thouless transition |KT73|. In reality, there is always some coupling 
of the 2D planes by an electron hopping tj_. At small temperature, the coupling between the planes 
and the strong correlations inside the plane will result in three-dimensional long range order^. The 
same remark applies to the SDW- and all other instabilities which are discussed in this thesis. 

While the interpretation of Tc in terms of superconductivity is very natural for attractive interac- 
tions, it is conceptually more problematic in the repulsive case. The possibility of a superconducting 
state in a purel y repu lsive system without any phonons, which was pointed out already by Kohn 
and Luttinger [ KL65| , is the main subject of Chapter || of this thesis. It has to be mentioned 
however that the interpretation of Tc in terms of Kohn-Luttinger superconductivity is not accepted 
by everybody |And|. 

The discussion of the superconducting instability in Chapter ^ proceeds with a minimum of 
formalism. It is based on the fact that in general there is only one kind of dominant fluctuations, 
the superconducting ones. The problem is more difficult in situations where density fluctuations 
are equally strong (actually diverging) as the superconducting ones. This can arise in 2D, due to 
the effects of the lattice. 

In fact, electrons hopping between the sites of a square lattice yield a spectrum that differs 
in two respects from the parabolic spectrum of the jellium model. First, the spectrum exhibits 
extrema and saddle points in the Brillouin zone. General considerations imply that there are at 
least two saddle points and two extrema (one maximum and one minimum). Obvious points are 
Po — (0,0) and Q = (tt, tt) for the extrema and Pi = (tt, 0) and P2 = (0,7r) for the two saddle 
points, but more complicated patterns are also possible. The density of states has a logarithmic 
van Hove singularity at the saddle points, in strong contrast to the constant density of states of the 
parabolic spectrum of the jellium model. The second difference is the shape of the lines of constant 
energy. These are circles in the case of the jellium model, whereas in the case of the square lattice 
one can easily find portions with almost vanishing curvature. In fact, for the tight-binding model 



*For a review on the Mott transition, includin g the d i stinction between Mott and Slater transitions, see |Gek]. 
^This was shown explicitly for the XY model IJM90l |BBCZ92|. 
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(with hopping restricted to nearest neighbor sites) the Fermi surface for the half-filled band is a 
perfect square. 

The van Hove singularities and flat Fermi surfaces lead to strong fluctuations both of the super- 
conducting and the density type. For the one-dimensional electron gas the same complication has 
been solved successfully by the renormalization group (RG) method. 

RG concepts have been used in very different fields of modern physics and RG can have quite 
different meanings for different people. The notion was first introduced by Stiickelberg and Peter- 



mann |SP53| and independently by Gell-Mann and Low [GML54| in the context of quantum field 
theories (like QED) in order to cope with infinities that appear in naive perturbation theory. 

In the early 1970s, Kadanoff and Wilson have associated the RG to the procedure of mode 
elimination in classical statistical mechanics (or systems of bosons) |Wil71, WK74|. In Wilson's 



formulation, the RG transformation consists of integrating out some degrees of freedom of the system 
and including them in the renormalization of some parameters (for example coupling constants). 
This alternative formulation of the RG idea proved tremendously successful in analyzing the critical 
behavior in the vicinity of second order phase transitioufj^ 

In general, a RG transformation is some change of the length or energy scale and the RG equation 
describes the response of the system as the length or energy scale is changed. In the interacting 
electron problem considered here, the energy scale is given by the temperature or alternatively by 
some cutoff A in the band energy of the electrons. The single-particle states with energies far from 
the Fermi level, i.e. with lattice momenta far from the Fermi surface, are subsequently integrated 
out. The result is an effective theory for the degrees of freedom at the Fermi surface. It turns 
out that the effective coupling constants of the low-energy effective theory have an alternative 
interpretation in terms of two-particle Green's functions with an infrared cutoff A. In this sense, 
the RG equation describes the change of Green's functions as a function of the energy scale A. In 
many cases, the RG flow to low energies produces a singularity at a finite energy scale Ac. This 
is interpreted as a transition into a strongly correlated state. The latter can establish long-range 
order, if stabilized by a hypothetical coupling of the 2D planes in three dimensions. The energy 
scale Ac is then comparable to the transition temperature towards the ordered state. 

It turns out that the one-loop (i.e. the leading order) approximation of the RG is equivalent to 
the so-called parquet approximation. This method was developed by the soviet school | PSTM56 | 
in order to treat different diverging fluctuations on an equal footing. It was successfully applied 



to one-dimensional conductors |BGD66, DK72, GD74], to the Kondo problem |Abr65| and to the 
X-ray absorption edge singularity problem in metals [RGN6E, >JGR69|. A detailed description of 
the method can be found in [RGN6£|. In my derivation of the one-loop RG equations in Chapter 
^, I will follow mainly the parquet philosophy, but the relation to the Wilsonian interpretation in 
terms of an effective theory depending on the energy scale will also be explained. 

The application of renormalization group (RG) ideas to 2D fermionic problems started in the 
mathematical physicists community [FT90, BG91] and has led to a considerable progress during 
the last decade]]. In fact, the rigorous works [ DROC| , ^al9S | mentioned above are based on similar 
ideas. 

A numerical scheme for calculating the complete flow from the bare action of an arbitrary mi- 
croscopic model to the low-energy effective action as a function of a continuously decreasing energy 
cutoff A has been presented by Zanchi and Schulz | Zan96| , ZS97, ZS98, ZSO0[ . The application of 
this method to the Hubbard model near half filling does provide an appealing picture, namely a 
transition from an antiferromagnetically ordered ground state at half filling to a d-wave supercon- 



ductor upon doping |ZSOO|. This result has been confirmed by Halboth and Metzner using a similar 



approach |fIMOO|. Recently, numerical RG calculations have brought up two additional phases, one 



^see |Wil75| for an early review including als o a, rlp t. ailerl H i sniissinp of the Kondo problem 

^Pedagogical introductions can be found in [Pol92, Sha94, |CFS9!:1. A precise r e lations hip between Fermi liquid 

nCSQq . |Dup94 |DupOC| ]. 



theory and the renormalization group has been established in |pS95 



4 



CHAPTER 1. INTRODUCTION 



with a deformed Fermi surface (Pomeranchuk instability) |HM01| and one with suppressed uniform 
spin and charge susceptibiHties ("insulating spin liquid") |HSFR01]. 

The analytical approach presented here is complementary to these numerical RG calculations. I 
start from the same RG equations and analyze the flow in the limit of small A. The RG equation is 
reduced to its leading order terms in the low-energy regime A — > whereas in the regime of higher 
energies, the only control of the flow is by naive perturbation theory. In this sense the analytical 
approach is limited as compared to the numerical studies, which are formally not restricted to 
small values of A. On the other hand, the numerical methods suffer from the need to replace the 
continuous Fermi surface by a discrete set of points. 

It is argued in Chapter |^ of this thesis that the reduction of the one-loop RG equation to the 
leading terms in A ^ is imperative for a consistent treatment. In fact, it turns out that sub- 
leading contributions to the one-loop equation are comparable to higher order terms, which have 
been neglected in the one-loop approximation. 

In our approach - as well as in previous RG calculations carried out to one loop order - self-energy 
effects are neglected. While this can be easily justified for the jellium model, the argument is more 
subtle in the case of lattice fermions. In fact, the second-order contribution to the self-energy is 
infrared divergent in the case of the half-filled nearest-neighbor tight-binding band. Nevertheless, I 
find that also in this case self-energy effects are of subleading order in A, provided that an instability 
(superconducting or density wave) occurs. 

Four distinct cases are considered in this thesis. First the generic case with a finite density of 
states and no Fermi surface nesting. Second, an anisotropic model with different hopping parameters 
tx and ty for electron hopping between nearest neighbors in the x- and y directions. The Fermi 
surface of this model is perfectly nested at half filling. Third, a situation where the Fermi surface 
passes through the saddle points of the dispersion at Pi = (7r,0) and P2 = (0,7r), without being 
nested. Finally the main emphasis is on the half-filled nearest-neighbor tight-binding model, where 
the Fermi surface is completely flat and contains two saddle points. 

RG calculations for a model where the Fermi surface contains flat portions have been performed 
by various authors |VdAD97| , |VdADOl| , |DVdAD02| , |Kwo97| , |ZYD97| , |GGV97a| . They agree in that a 
d-wave superconducting instability occurs for repulsive interactions, due to the interplay of particle- 
particle and particle-hole correlations. 

Early scaling approaches to the problem of van Hove sing ularities |Sch87| , [DzySTj |LMP87| fo- 
cussed on the interactions between electrons at the saddle poi nts, by treating these points in analogy 
to the two Fermi points of the one-dimensional electron gas |^ol7E |. In this thesis I show that, in- 
deed, the logarithmically dominant RG flow at low energies is controlled by the neighborhood of the 
van Hove points. However, in contrast to the one-dimensional case where the scattering processes 
can be characterized in terms of a few coupling constants connecting the two Fermi points, in two 
dimensions the effective couplings are functions of incoming and outgoing momenta, even if these 
are restricted to the Fermi surface. It turns out that this functional dependence plays a crucial role 
in the asymptotic decoupling of competing instabilities. A step in this direction has already been 
made in the parquet approach of Ref. | DY88a , DY88b| . 

When the Fermi level is at a van Hove singularity the system is not renormalizable in the 
traditional sense of field theory. Nevertheless, electrons near a van Hove singularity have been 
treated b y applying the field theory formalism to the particle-hole sector [ GGV96 , GGV97b , |GGV99 , 
GGVOOj. No mixture with particle-particle diagrams can be treated within this formalism. The 
Wilsonian RG used here does not assume renormalizability and may be applied without constraints. 

The case where the Fermi surface contains the van Hove singularity without being nested has 
been addressed in many recent investigations [[KKOl, IKOl, HSOIa, HSOlt]. It turns out that 
(in the case of a repulsive interaction between the electrons) the leading order terms in the small 
parameter A is not sufficient to obtain the full phase diagram. The abovementioned works are hence 
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mainly based on subleading contributions to the RG flow|^. It is not entirely clear to me whether 
this can be done consistently within a one-loop approach. 

In contrast, the analysis of the dominant parts of the RG equations is sufficient for establishing 
a rich phase diagram for the nearest-neighbor tight-binding band with a nested Fermi surface. This 
will be done in Chapter ^. 



e. simply logarithmically diverging terms of the perturbation theory. In contrast, the leading terms diverge as 
log^ A. 



CHAPTER 1. INTRODUCTION 



Chapter 2 



Superconducting instabilities of a 
general 2D Fermi surface 



This chapter is organized as follows. After a brief presentation of some basic notions and defini- 
tions in Section I will illustrate the breakdown of naive perturbation theory due to infrared 
divergences in Section 2.2. Section 2.3 presents the so-called ladder approximation, which consists 
of a consistent summation of the leading diverging terms. It is shown in Section 2.4, how this 
approximation indicates superconductivity in the extended Hubbard model. Finally the Section 
2.5, where it is shown that the ladder summation is not sufficient if the Fermi surface is nested, 
serves as motivation for the more formal investigations of Chapter ^. 



2.1 Interacting electrons on a lattice 
2.1.1 The model 

We consider a general single-band model of interacting electrons on a two-dimensional square lattice. 
Formally, the size of the lattice is L x L with periodic boundary conditions, but all the calculations 
are done in the thermodynamic limit L — > oo. The lattice spacing is put equal to unity. 

The single-particle states consist of one localized Wannier state for each lattice site r and they 
are labeled by a spin index a =t, I- cj^ and c^.^ are the usual creation and annihilation operators. 
The Hamiltonian is of the form H = Hq + Hj, where 

Ho= J2t{r,r')clc,,^ (2.1) 

r,r',cr 

is the non-interacting part and 

"i = l C/(ri,...,r4)^4,.4..'Cr3.'Cr,. (2.2) 

ri,...,r4 a,(T' 

is the most general two-body interaction respecting the global spin rotation symmetry and particle 
number conservation. 

Translation-invariance is also assumed, i.e. t{r,r') = t{r — r',0) and U{ri, . . . ,r4) = U{ri — 
r4, r2 — r4, — r4, 0). After Fourier transformation Ck^o- = 1/^ X^r ^~^^'^'^r.a, Hq is diagonal 

^o=5]ekcLck, (2.3) 

k,cr 
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with the single-particle dispersion ek = J^^^^^^^i'^^^)- "^^^ properties of ek are crucial for the 
behavior of the system at weak coupling and different specific cases will be discussed in this thesis. 
The interaction after Fourier transformation reads 

J2 3(1^1' 1^2, kg) ^4^^4^^,Ck3^,Ck^+k2-k3^' (2-4) 

ki,k2,k3 a, a' 

where 

5(ki,k2,k3)= J2 e^^'^i''i+''=''=-''^''^)t/(xi,X2,X3,0). (2.5) 

Xi,X2,X3 

The momenta k are only defined modulo a reciprocal lattice vector (i.e. mathematically speaking 
k G ^Z^/27rZ^) and throughout this thesis it is implicitly assumed that sums, equations and 
Kronecker symbols involving momenta are taken modulo 27rZ^. For example the equation 2p = 
has four solutions, p = (0, 0), (tt, 0), (0, tt) and (tt, tt). The summations are by convention over 
the first Brillouin zone ] — tt, tt] x ] — tt, tt] . The hopping matrix as well as the interaction is assumed 
to be short-ranged in real space, which leads to smooth functions ek and (?(ki, k2, ka) in momentum 
space. The Hamiltonian was written in a very general form for the following reason. Even if one 
starts with a simple model such as the Hubbard model, the effective models which are generated 
by the renormalization group procedure are of the general form (^^) and ( ^.4| ). 

Nevertheless the effective model will always satisfy the symmetries of the original model, at 
least if these symmetries are not spontaneously broken. It is therefore worthwhile to list some 
basic symmetry relations in addition to particle number conservation, spin rotation invariance and 
translation invariance already mentioned above. 

Time reversal invarianceQ requires g{ki , k2 , ka) = $((— ki , — k2 , — ka) , where 5 means the complex 
conjugate of g. Together with the spatial reflection invariance this means, that g is real. In addition, 
we require 

(7(ki,k2,ka) = g(k2, ki, ki + k2 — ka) permutation symmetry (2-6) 
= g(ki -f k2 - ka, ka, k2) Hcrmiticity (2.7) 
= 3(i?ki,i?k2,i?ka), (2.8) 

where i? is a point symmetry operation of the lattice. 

There is no symmetry with respect to the operation Xg(ki, k2, ka) = 5(k2, ki, ka) = g(k-i, k2, ki-f- 
k2 — ka). The symmetric part g^ = ^{1 + X)g and the antisymmetric part g^ — X)g describe 

scattering of singlet and triplet pairs, as becomes clear if we write the interaction as 

^ kjl |/(ki,k2,ka)4(k2,ki)0s(ka,ki + k2-ka) 

kl,k2,k3 

+g^(ki,k2,ka) V </)l(k2,ki)^Jka,ki+k2-k3)l, (2.9) 



a=0,±l 



with 



0s(k,k') = ^(ckTCk'i -CkiCk't), (2.10) 
0o(k,k') = ^ (ckTCk'i +CkiCk'T) , <?!'i(k, k') = CkjCk't , 0-i(k, k') = CkiCk'x . 
Sometimes one separates the term with parallel spins {a ~ a') from the term with anti-parallel 



spins in Eq. (2.4). It is clear that the coupling function for electrons with anti-parallel spins is 5, 
whereas for parallel spins it is g"^. 



'^See appendix 
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Alternatively, via the combinations g'^ — (2 ~ X)g and g'^ = —Xg one can write Hj in terms of 
normal ordered charge-charge and spin-spin interactions 

Hi = 2^ ^ ^g"" 0^1,^2,^3) ■ C'kakaC'ki ki+k2-k3 : k2, ka) : Skakg'Ski ki+ks-k., :| : 

ki,k2,k3 

(2.11) 

where normal ordering is given by : c\c^c\c^ : — c\c\c^c^. C and S are given by 

Ckw = ^^cl^c^,^, Skk' ^ ^^cl^ff„„,cj^,^,, (2.12) 



where a are the Pauli matrices. 

For example a general interaction restricted to on-site and nearest-neighbor terms consists of 5 
independent terms. The usual Hubbard term U, a nearest-neighbor density-density interaction V, 
a Heisenberg term J, a pair-hopping term W and the Coulomb- assisted hopping K 

Hi = U^ni^rnir + V ^ HrUr' + J ^ SrSr' +W ^ (cl^cl^c^,^c^,.^ +h.c^ 

r (r,r') (r,r') (r,r') 

+K J2 (cLK-o-+nr'-a)Cr'„ + h.C.), (2.13) 

where = Y,a^^<y = J^a'^ia^ra and Sr = ^J2aa' '^tcr^cc'^^ra' ^"^^ usual charge- and spin 
operators on the lattice site r, the sums J2{r r') '^^^^ nearest- neighbor bonds and U, V, J, W 
and K are parameters. After Fourier transformation, this leads to Eq. (2^) with the coupling 
function 

5(ki,k2,k3) - ;7+(l^-J/4)/k3-k.-J/2/k3-ki+Wki+k.+i^(/ki+/k.+/k3+/ki+k2-k3), (2.14) 
where /k = 2(cosfc2: -I- cosfcj^). 

2.1.2 Green's functions 

Since we are concerned with weak coupling instabilities, the first idea is perturbation theory in Hi. 
The diagrammatic technique using the imaginary time (or Matsubara-) formalism is well established 
and explained in various textbooks (for example |N088i |AGD61| , |FW7l| , |Mah8l[ ). 
The partition function in the grand canonical ensemble is given by 

Z = TTe-^^"~^^\ (2.15) 

where /3 is the inverse temperature, /i the chemical potential and N — X^ro- 'ia^ra ^^e particle 
number operator. The trace is to be taken over the entire Fock space. 
Thermal averages of physical observables are calculated by the formula 

(O) = ^Tre-'3(^-''^)0. (2.16) 

The aim is to calculate the imaginary time Green's functions, which are not directly observable but 
they are analytically related to various physical response functions. 
They are defined by 

G"(xi, . . . ,x„|x2„, . . . ,a;„+i) = (-l)"(T(c,^ • • • c^J^^^, ■ ■ ■ 4,J), (2.17) 
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where X = (r, r, a) is a multi-index containing the imaginary time t G [0, /?] in addition to the labels 
of the single particle states r and a, the imaginary time Hcisenberg representation of operators is 

0{t) = e^iH-f.N)Q^-r{H-f.N)^ (2.18) 

and T is the usual time ordering operator with respect to r. Note that the Hermitian conjugate of 
0(r) is Ot(-T). 

Because of the translation invariance in space and imaginary time, it is convenient to work with 

Fourier transformed quantities. This means passing from (r, r) to the 2 + 1-dimcnsional frequency- 
momentum vector k = (fco,k). It contains the Matsubara frequency Uq which runs over the odd 
multiples of 7r//3. 

The annihilation operators in frequency-momentum presentation are defined by 



Jo 



(2.19) 



The Fourier-transformed one- and two-particle Green's functions are then 

G{k) = / dr Ve^('=°^-'''')G(rr(7|0Oa) 
Jo 

= -{Tc^kclk) (2.20) 

(independent of a) and 

r0 ^ 

G"" {ki,k2,k3) = / dndradra ^ e'('="^^+'=''=^=-'=''=*^^-''i'"i-'*=''=+''=*''^)G(Tiria,r2r2CT'|00CT,T3r3£7') 

^ ri,r2,r3 

= /3i'(rc,,^c,-,,4„34^+,^_,3) (2.21) 
The single-particle Green's function of the non-interacting system is given by 

^^^'^ = -ir-^^ (2.22) 

where = fk — A*- It is also called free electron propagator or covariance of the fermion field. 
The full Green's functions are readily extracted from their one-particle irreducible (IPI) parts. For 
example the self-energy S(fc), defined as the IPI part of G{k) is related to the latter by Dyson's 
equation 

The IPI vertex function F"'"' (fci, fc2, k^), defined as the IPI part of the two-particle Green's function 
G""' (fci, fe, fcs), fulfills the relation 

4 

G'^'^'(/ei,fc2,A;3) =/3i'(^fe,fc3 -'^a,<.'<5fc.,fc3)G(/ei)G(fc2)+r'^'^'(fci,fc2,fc3)n^('^^)' (2-24) 

where k^ = ki + k2 ~ k^. 

At this point it is worthwhile to investigate the symmetries of the vertex function. A simple 
SU{2) transformation shows that for a spin rotation invariant system, 

rTT = (l-X)FT^ (2.25) 
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where XT{ki, k2, fcs) = r(/c2, ki, k^). Since the whole information is contained in the function F^^, 
we will only consider the vertex of anti-parallel spins from now on and omit the spin indices (i.e. 

r = rTi). 

The permutation symmetry gives 

T{ki,k2,k3) = r(fc2, ki,ki+ k2 ~ ks). (2.26) 
Time reversal invarianceQ, spin rotation invariance and parity imply 

T{ki,k2,k3) = r{ki + k2-k3,k3,k2). (2.27) 

From the behavior under complex conjugation one obtains 



I](fco,k) = S(-fco,k) (2.28) 

and 



T{ki,k2,k3) ^T{ki + k2 ~ k3,k3,k2), (2.29) 
where k ~ {—ko, k). It follows that both E and F are real if the frequencies are put to zero. 

2.2 Naive perturbation theory 

A calculation of the vertex function to second order in perturbation theory gives 

F(/ci, fe, fca) = -g(ki, k2, kg) + PP + PHI + PH2 (2.30) 

= E , k2, k - p)C{p)C{k - p)g{p, k - p, kg), 
p 

PHI = ^^5(ki,P + qi,k3)C(p)C(p + gi)5(p,k2,p + qi), 



p 



™2 = ^^C(p)C(p + g2)[-2.g(ki,p,p + q2).g(p + q2,k2,k3) 

p 

+g{p, ki, p + q2)5(p + q2, k2, kg) + ^(ki, p, p + q2)5(p + q2, k2, p)] , 
where k = ki + k2, qi = k^ — ki, q2 = ks — k2. The representation of the second order contributions 



in terms of Feynman diagrams is shown in Fig. 2.1. The internal electron lines stand for free 
electron propagators C and the wavy interaction lines stand for coupling functions g, that depend 
on the incoming and outgoing momenta. The convention is that the spin index is conserved along 
the fermion lines. The minus sign in the first of the three PH2 diagrams comes from the fermion 
loop and the factor 2 in the same diagram from the sum over the spin index in the fermion loop. 
PP is referred to as the particle-particle (p-p) diagram and PHI and PH2 as particle-hole (p-h) 
diagrams. 
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PH2: 





Figure 2.1: Second order diagrams for the IPI vertex T{ki, k2, k^). k^ is determined by energy- 
momentum conservation. 



2.2.1 Divergences in the particle-particle channel 

Suppose for a moment that g is constant (= U in the Hubbard model). The p-p diagram is then 
proportional to the p-p bubble 



B''P{k) = ^Y.^{p)C{k-p), 



(2.31) 



where fc = fci + ^2- The p-h diagrams are proportional to the p-h bubble 



(2.32) 



where g = fcs — fci for PHI and q = k^ — k2 for PH2. This section is concerned with the p-p bubble. 
Summation over the Matsubara frequencies leads to 



(2.33) 



where n(^) = (1 -I- exp /3^) ^ is the Fermi distribution function. This quantity is diverging at fc = 
as the temperature goes to zero. At a finite temperature T = j3^^ one finds 



1 r 



W 
Y 



tanh^ 
'finite", 



(2.34) 



where W is the bandwidth and iy{e) — l/L"^ Si_^i^k — e) is the density of states, and "finite" is any 
contribution with a finite limit T — > 0. Eq. ( |2.34 ) can be shown by taking the zero temperature 



limit of T^^-^y^^. The only two assumptions are reflection symmetry — ^-k and that the density 
of states v{^) is an analytic function of ^ at the Fermi level. 
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Under the same assumption it can be shown that the biggest term of order n in the couphng 
diverges as \og"--\W/T). It is clear that naive perturbation theory breaks down at large enough 
values of the couphng or low enough temperature, such that |5log{VF/r)| ~ 1. Alternatively, one 
can organize the perturbation series in the form 

r = ai{g\og{W/T)) g + a2{g\og{W/T)) g^ + ■ ■ ■ . (2.35) 

We only calculate the first term in this expansion, which means that diagrams of every order in g 
have to be calculated to the leading logarithmic precision and all the subleading contributions such 
as the term "finite" in ( 2.34 ) are neglected. In other words, g is considered small, but g\og{W/T) 
is not. 

At any finite frequency, the zero temperature limit of the p-p bubble can be taken and the result 

is 

pW ■ p 



lY 2^ - iko 

W 



1/(0) loi 



iko 



+ "finite". (2.36) 



The finite frequency just replaces the temperature in the logarithmic expression. 

For a more general coupling function g{\ii, k2, ks) the p-p diagram at fci + fc2 = is of the more 
complicated form 

1 ^ 



PP = ^ E C{p)C{-p)F{p), (2.37) 



p 

where the function F(jp) comes from the coupling functions and depends on the external momenta. 

Any higher order diagram which is obtained by replacing the two coupling functions in PP by 
a more complicated sub-diagram is of the same form, with a frequency dependent function F{p). 
The frequency summation leads to 

pp^ 1 J2 ^i^p)F{p)Up„=^,-{l-n{^p))F{p)Up,,=_^^ ^^^^^ 

We can now perform a change of variables from p to and an additional variable ^p, which 
parameterizes the curves of constant band energy in the Brillouin zone. Assuming that F is an 
analytic function of po and ^, one can Taylor expand 

F{co, 9) = FooiO) + Fw{e)iu; + Foi{9)C + ■■■. (2.39) 

It is clear, that only the first term contributes to the logarithmic divergence and the others can be 
neglected to leading logarithmic order. 



2.2.2 Cutting off the infrared divergence 

We have seen that naive perturbation theory breaks down at zero temperature. One can regularize 
the zero temperature perturbation theory by introducing artificially an infrared cutoff A. There are 
several ways to do this. One method to avoid the singularity is to replace the bare propagator by 
C^ik) = 6(|fco| — A)C(fc), where 6 is the Heavyside step function. This construction leads to 

J-W ^ 

w 

= log-r- "finite". (2.40) 
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We see that Bpp{T, k — Q) — B^{T — 0,k — 0)|a=t to leading logarithmic precision. Alternatively 
the cutoff can be introduced in the band energy via CA(fc) = 9(|'?k| — A)C(/c), leading to 

BHO) = i (/> + />) f (2.41, 

w 

= iy{0) log^ + "finite". (2.42) 

The logarithmic term is always the same. One could also use the quantity \/fc§ + to introduce the 
cutoff or replace by a smooth function without changing the result to leading logarithmic order. 
We are thus free to calculate at zero temperature, introduce the cutoff in the most convenient way for 
calculations and in the end of the calculation replace A by T. Useless to say that the correspondence 
A ^ T is only correct to the leading logarithmic order. In a more elaborate calculation including 
subleading contributions, the result depends on how the cutoff is introduced. 



2.3 The case of a general Fermi surface: Ladder diagrams 
and BCS theory 



We now calculate the IPI vertex to leading order in the expansion (2.35). In a general situation 
without any fine tuned parameter, the only diverging terms occur in the p-p diagram. The other 
second order terms have a finite limit T, A — > exce pt for special situations, which will be discussed 
later. The leading term in the development ( 2.35| ) is then given by the sum of all the p-p ladder 
diagrams, shown in Fig. 2.2. 



ki+k2- k3 



Figure 2.2: The class of p-p ladder diagrams. They give the leading logarithmic contribution to the 
vertex in general situations. 

The analytical expression of this series is 

rA(fcl,/c2,/C3) = -g(ki,k2,k3)-H^^5(kl,k2,k-p)Z?PP,(p)(?(p,k-p,k3) (2. 

p 

^ ^ ^ ' P)^'''^ ^^P' ^ - P' ^ - p')DTAp') 9{P', k - p', ka) + . 

P:P' 

where k = ki+k2 is the total frequency-momentum and D^i^{p) — CA{p)C\{k—p) is the propagator 
of a pair of particles. 



In principle, the internal electron propagators in Fig 2.2 could be dressed by self-energy cor- 



rections. Such corrections have three main effects. First they change the shape and location of 
the Fermi surface, second they modify the properties of the single particle dispersion (the Fermi 
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velocity) and third they lead to a reduction of the quasi particle weight. The deformation of the 
Fermi surface requires in general the introduction of counter-term^ I will ignore this effect for 
the moment and assume, that in the weak coupling limit, the interacting Fermi surface is not so 
different from the non interacting one such that the difference is not crucial for the study of the 
weak coupling instabilities. In fact, perturbative corrections to the Fermi velocity VpE(p) and the 
quasi particle weight z — {1 + idpgT,{p))~^ are finite as T and therefore do not contribute to 
the leading order. 



The total frequency-momentum k is the same for all vertices entering Eq. (2.43). It is therefore 
useful to define the functions 

rA,fe^(fci,fc4) = rA(fci,fc-fci,A:-fc4), 

3^^^(ki,k4) = 5(ki,k-ki,k-k4). (2.44) 



and rewrite Eq. (2.43) as 



r^Tiki^k,) = -.gr^(ki,k4) + ^ ^gr^(ki,p)D^^,(p)5r^(P,k4) - • • • (2.45) 



In analogy with the discussion at the end of Section 2.2.1 we emphasize that after the change of 
variables p ^ ^,9, the dependence of g^'^^ on the frequencies fcoi,fco4 and on the band energies 
^1, ^4 is not relevant to leading logarithmic order. One can replace r^^'^(fci, and g^'~'^ {ki, ^4) 
in Eq. ( 2.43| ) by its value at fcoi = fco4 = Ci ~ ?4 = and write 



rfTiSi,04) = -.9^^(01,04) + d0gS^'{e,,9)Bl^,{e)gE^^{eM - • • • (2.46) 



where we have introduced the angle-resolved p-p bubble 

KkW = 4E [d^.H0,ODAAp)\p=Pim, (2.47) 

with the Jacobian J(0, ^) coming from the change of integration variables. 
For k — T = Q with the energy-cutoff A one obtains 




"^de ) ^ (2.48) 
^ \J-W JK J i 

and thus 

from which the logarithmic behavior 

Sa(6') = J{e, 0) log ^ + "finite" (2.50) 

can be deduced. The limit limA^o >^(^7 A) can be safely taken, provided that no van Hove singularity 
is present. As shown in Section A can be replaced by Min{A, T, fco}, i.e. the cutoff A can be 
interpreted as the temperature or a finite frequency to leading logarithmic precision.^ 



See for example Section 5.7 of Nozieres' book [ Noz6^ for a mmprehensi ye exply i riation iigornii s mathematical 



FST9 



statements about the moving Fermi surface have been presented recently in |FST96, FST98, 

■^The dependence on the momentum k is less simple and depends on the form of the Fermi surface. For example 
if the ladder has fiat portions, then the p-p bubble diverges for finite values of k parallel to the fiat Fermi surface. 
At the same time, the flat portions introduce divergences in the p-h channel. 
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Introducing the measure d/i(6') — d9J{6,0), Eq. (2.46) can be further simphfied to 

rr^(0i, 04) = "5^'"^(0i, 04) + log ^ J d^,{9) 9^^^{9,,0)g^^HO, 9,) - ■ ■ ■ (2.51) 

Note that the integral over the Fermi surface J dfi{6) equals the density of states at the Fermi level 
u{0). 

Eq. (2.51) is like a geometric series of functions g^^^{6,9'). In order to sum the series ( |2.51| ), 
we consider g^''^^ as the kernel of an operator acting in the space of functions f{6) as 

f{9)^!^ J dfi{9')g^^''{9,9')fie'). (2.52) 

This operator can be diagonalized. This means that g^'~^^ has the spectral representation given by 

5^^5(0,0') = ^A„/„(0)/„(0'), (2.53) 

n 

where A„ are the eigenvalues of the operator and fn{9) the corresponding eigenfunctions. The latter 
satisfy the orthogonality relation 

' dfi{9)f„{9)fmi9)=S„m. (2.54) 

It is important to note that the eigenvalues An not only depend on the interaction g, but also 
on the properties of the band structure. For example in the simple case of a constant coupling 
gBCS(^0^ 0') = U (th e Hubbard model), the only non-zero eigenvalue is given by A = v{0)U. 
The series ( ^.51 ) can now be summed easily 

rr^(0,0') = ^/„(0)/„(0')(-A„ + A^log^-A^V^ + ...) 
„ -A„ -log — 

Clearly, negative eigenvalues produce a singularity in the vertex as A is lowered to the energy scale 

Ac^We^. (2.56) 

The divergence of the vertex F^*-^"^ at A = Ac is usually interpreted as a signature of the onset 
of superconductivity. First it has to do with p-p pairs and it requires some attractive channel in 
the interaction. This suggests that the instability is associated with the formation of Cooper pairs. 
Moreover formula (2.56) for the critical energy scale corresponds exactly to the BCS formula for 
the critical temperature. To make the connection to BCS-theory, consider the gap equation 

= E.9^'''(k. k')^ tanh(^). (2.57) 



Close to the critical temperature, the mean-field dispersion i?k' = \/ -I- is replaced by |^k|- 
Furthermore one can neglect the dependence of A and g on the band energies and ^k' in the 
weak coupling limit. Using Eq. (2.34), this leads to 

A{9) = - y" d^' g^^^iO, 0')A{9')Bt^ (0'), (2.58) 
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where the angle-resolved p-p bubble has been defined in Eq. (2.47). To logarithmic precision, 
BtAO') = J{0,0)logiW/Tc). One finally obtains 

Ai0) = - iog{w/T,) J d^l{0') 5^^^(0, e')A{e'). (2.59) 

This is precisely the eigenvalue equation for the operator g^'^^ introduced above. We conclude that 
the eigenfunctions f{0) correspond up to a normalizing factor to the gap function (s-wave, p-wave, 
etc.) and find the correspondence — A^^(O) = log(Vl^/Tc) for the eigenvalues. 



2.4 Superconductivity from repulsive interactions 

The theory presented above predicts a superconducting in stabili ty, if one of the eigenvalues of g^'^^ 
is negative. However, according to Kohn and Luttinger |KL6£], superconductivity occurs for any 
short-range interactions, even if it is purely repulsive. 



The ladder series (2.51) satisfies the following differential equation 



AaArf^^( 



'3)= /dMe)rr^(0i,« 



^BCSf 



0- 



(2.60) 



This equation describes how the vertex changes as we vary some energy scale, i.e. the temperature, 
the frequency or the artificially introduced cutoff. Quite generally, the map Fa FA-dA is called 
a renormalization group (RG) transformation and equations such as Eq. (2.60) are called RG 
equations. 

A good strategy for the calculation of the vertex at low energy is a combination of naive per- 
turbation theory and ladder approximation. First use simple (second order) perturbation theory 
to calculate Fao . Then integrate Eq. ( 2.60| ) from Ao down to lower energies. The energy scale Aq 
must be chosen not too small, such that perturbation theory works, i.e. gB^^ ^ 1. On the other 



hand it must be small enough, such that the ladder approximation is justified, i.e. B^^ <C B 



ph 



Eq. ( 2.60 ) can be solved using the same strategy as for the su mmat ion of the series ( 2.51 ), 
namely by w riting F^'-^'^ in a spectral representation similar to Eq. ( 2.53 ). Because of the special 
form of Eq. ( 2.60[ ) only the eigenvalues, but not the eigenfunctions, depend on the energy scale A, 

9')= V7n(A)/„(^?)/„(^?')- (2.61) 



1. e. 



^BCSi 



This ansatz solves Eq. ( 2.60 ), provided that the eigenvalues 7n(A) satisfy the RG equation 



AaA7n(A)-7f.(A), 



and thus are given by 



7n(A) 



7r7'(Ao) 



log 



(2.62) 



(2.63) 



The vertex calculated by this method does not correspond to the ladder series ( pl3| ) but to 
a similar series, where g is replaced by — F^^p'^ and 15a is replaced by Z3a ^ ^Ao- Therefore the 



eigenvalues A„ in Section 2.3 are replaced by — 7„(Ao), and log ^ is replaced by log ■ 



I have calculated the superconducting instabilities for an extended Hubbard model away from 
half-filling. The electron hopping is restricted to nearest neighbors, i.e. — —2t{cos k^+cos ky)— fi. 
Fig. 2.3 shows the Fermi surface for a typical electron density. The variable 9 is defined as the 



radial angle. The interaction contains a nearest neighbor term V in addition to the Hubbard U 
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Figure 2.3: The Fermi surface of the nearest-neighbor tight binding band for /i = —0.8 (i.e. n « 0.7). 



(See ( 2.13 ) and ( 2.14 )). The same calculation was done before by Zanchi and Schulz | ZS96 for the 
repulsive Hubbard model. 

The study of superconducting instabilities involves two steps. First calculate the vertex r^^'^(6'i, 6*3) 
by second order perturbation theory. The calculations presented here have been done for Aq = 
Min{— /1/4, + 4)/4}. To make a numerical treatment feasible a discretization of the 9 variable 
into 24 patches was introduced. 



The second step consists in writing F^^'^ in the form ( 2.53| ) by diagonalizing the corresponding 



operator. This is done by choosing an orthonormal basis bn{9). If the basis is suitably chosen, the 
infinite matrix 



MO) 



(2.64) 



The 



can be cut off at some finite value of the indices n, m and diagonalized by standard algorithms, 
biggest positive eigenvalue of F„.„j gives the leading superconducting instability. 

I have used the harmonic functions cos{n9) and sin{n9) divided by the square root of J{9, 0) as 
a basis. Due to the symmetry of the square lattice, the matrix has a block-diagonal form, where 



each block corresponds to an irreducible representation of the point group D4. Table 2.1 lists the 
five irreducible representations and the corresponding basis functions. 



Irreducible representation 


Basis functions 


Ai 
A2 
Bi 
B2 

E 


cos 4m(? 

sin 4:m9 

cos (4m -1- 2)9 

sin (4m -|- 2)9 
( cos (2m + 1)9 
I sin (2m + 1)6* 



Table 2.1: Table of the irreducible representations of the D4 point group and the according basis. 
In order to obtain an orthogonal system, the functions are divided by -y/ J{d, 0), but this doesn't 
change their symmetry properties. 



All the calculations have been done for a repulsive on-site interaction U ~ Q.5t. In the case of 
an attractive < 0, the appearance of superconductivity is not surprising. In fact, the first order 
vertex F = —g has positive eigenvalues in this case. The results can be seen in Fig. 2.4. They 



show triplet p-wave superconductivity for n < 0.65 and singlet dx2_y2-wave superconductivity for 
0.65 < n (close to half filling, additional instabilities arise in the p-h channel). The second order 
contributions to T^^^{di, 6*3) give only minor changes to the present result, although they introduce 
small positive eigenvalues in the symmetry blocks B2, Ai and A2. 
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Figure 2.4: The positive eigenvalues of T^'~'^ as a function of the electron density for an attractive 
nearest-neighbor interaction V — —0.05t and U — Q.bt, as obtained by a first order calculation. The 
only positive eigenvalues are in the E and Bi symmetry block. The corresponding eigenfunctions 
are p-wave like (for E) and (i^2_j,2-wave like (for Bi). 



In contrast, in the pure Hubbard model {V = 0, U > 0), positive eigenvalues of Tf^^^ {61,93) 
appear only via second order terms. Since the sum of three diagrams PH2 is zero in the Hubbard 
model and PP is merely a constant (i.e. independent of the angles 6,6'), positive eigenvalues can 
only be generated by t he dia gram PHI. Kohn and Luttinger have investigated the same diagram in 
their historical paper [KL65|. The difference is that we are dealing with an anisotropic 2D system 



as opposed to the isotropic 3D case considered in |KL65 
The results in Fig 



2.5 for the Hubbard model show dxy-wave symmetry for n < 0.55 and 



ci3,2_ 2-wave for higher densities. Note that the eigenvalues of Fig. 2.5 created by second order 



terms are much smaller than those in Fig 2.4 with an attractive V, leading to a superconductivity 
at much lower energy scales. The eigenvalues of the Bi block (dxy) and the B2 block {dx2-y2) are 
degenerate in the limit of low electron density, reflecting the approximate rotational symmetry at 
low filling. 

The appearance of superconductivity by second-order corrections to the vertex has a physical 
interpretation in terms of an effective attractive interaction, which is mediated by spin fluctuations. 
In fact, the vertex Fa can be interpreted as an effective coupling function gjy = —Fa- This point 
will be made more precise later in Section 3.4.2. The contribution of PHI to gf^'^'^{di, 63) equals 



U^BP'\ke, + ke,) = U\o{^0^ + k^s), 



(2.65) 



where Xo(q) is the spin susceptibility of the non- interacting systein|^ The momentum dependence of 
the effective interaction (and thus the attraction) comes therefore from the spin susceptibility Xoil)- 
One can thus say that the effective attraction is created by the exchange of spin fluctuations. The 
situation is analogous to the conventional superconductivity of metals, where an effective attraction 
is created through the exchange of phonons. However with the important difference that the same 
electrons which feel the effective interaction are also responsible for the spin fluctuations. 

The idea of spin-fluctuation-induced superconductivity has led to semi-phenomenological theo- 
ries of materials with strong magnetic correlations such as the cuprates which are close to antifer- 
romagnetism |BSW89|, pSSl |MTU90|, |MP92|, |KM99|, but more recently also for Sr2RuOi which is 



nearly ferromagnetic [ ML99[ 



^See Section 



or standard textbooks. 
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Figure 2.5: Second order calculation for the pure Hubbard model V = and U — O.bt (Kohn- 
Luttinger superconductivity). The biggest eigenvalue of each symmetry block is shown. The corre- 
sponding eigenfunctions are dj^^-wave like for B2 and dx2_y2--waYe like for Bi. 
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Figure 2.6: Same as Fig. |2.5| for repulsive nearest-neighbor interaction V — 0.05t and U — 0.5t. 
The leading eigenvalues from the Ai and A2 blocks correspond to highly oscillating eigenfunctions 
with eight nodes along the Fermi surface. 



2.5. DIVERGENCES IN THE PARTICLE-HOLE CHANNEL: NESTING 
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If a positive nearest-neighbor term > is added to the Hubbard model, all second order 
contributions PP, PHI and PH2 to Pf^^ (6*1, 6*3) are non zero and angle-dependent. As it is shown 
in Fig. |2.6| , the V^-term has a destructive effect on superconductivity. The effect is strongest in 
the d^2_y2-wa,Ye channel (corresponding to the Ai symmetry). While the dxy regime for n < 0.5 
is not affected much by the nearest-neighbor repulsion, the d^2_y2-wa,ve regime has completely 
disappeared. The leading instabilities arise in two highly oscillating exotic channels instead. 

These results indicate strongly, that there is indeed superconductivity in the repulsive Hubbard 
model, at sufficiently low temperature. Furthermore they illustrate how the RG equation Eq. 
( 2.60| ) can be successfully used to investigate the weak-coupling instabilities of 2D lattice electrons. 



Simple perturbation theo ry is used to calculate the vertex at an energy scale Aq. They serve as 
initial conditions for Eq. ( 2.6C| ). Superconducting instabilitie s man ifest themselves as poles in the 



vertex function, or equivalently as diverging solutions of Eq. ( ^.60 ) 



2.5 Divergences in the particle-hole channel: Nesting 



In Sections 2.3 and p.4| , it was assumed that the divergences in the perturbation series of the vertex 
are exclusively generated by the p-p bubble. 

The divergence of BPP{k) at k = depends on the parity relation 

Cp = e-p. (2.66) 

Since parity is a general symmetry, this divergence is generally present. 

In contrast, divergences in the p-h channel only arise in special situations. The p-h diagrams 
are proportional to the p-h bubble 

B^^k) = J-J2CA{p)CA{p + k) 

p 

where k = — ki in PHI and k = — ki in PH2. The p-h bubble is not diverging at fc = 0, 
because of the numerator. In fact it is easy to see, that 

hm i3^"'(0,k) = J d^i^iO^ - KO) (2.68) 

and SP''(fco,0) = 0. 

Consider now the special case of the half-filled nearest-neighbor tight-binding band. In order to 
be more general, one can allow for an anisotropy in the hopping parameter, i.e. two different values 
—tx, —ty for hopping between nearest neighbors in the x- and y-directions. The electron dispersion 

= —2tx cos — 2ty COS Py Ilas the special property 

Cp = -?P+Q (2-69) 



for the nesting vector Q = (tt, tt). The Fermi surface is shown in Fig. 2.7. The immediate 
consequence is the following relation between the p-p and the p-h bubbles 

BP''{k)^-BPP{Q-k), (2.70) 



where Q = (0, Q). Therefore, everything that was said before about the p-p bubble at momentum 
A; = applies as well for the p-h bubble at momentum k = Q. 
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Figure 2.7: The Fermi surface of the nearest-neighbor tx-ty-hand. It is perfectly nested by the 
vector Q = (tt, tt) 



The condition (2.69) is in general too restrictive, i.e. it is not a necessary condition for diver- 
gences in the p-h channel. We now formulate more general conditions for divergences both in the 
p-h and p-p channels. They depend only on the geometry of the Fermi surface and are valid for 
every periodic lattice in any dimension. It is then useful to remember that the Fermi surface divides 
the Brillouin zone in two parts, the Fermi sea with the occupied single particle states (^k < 0) and 
the remaining part with the empty states (fk > 0). 

• B^^ is diverging at a fixed momentum k, if: 

1. The Fermi surface FS has some finite overlap with its own translation k + FS. 

2. The overlap is such that the occupied side of FS is covering the empty side of k -I- FS 
and vice versa. 

This property is called p-h nesting with a nesting vector k. For example every bipartite lattice 
(i.e. a lattice with electrons hopping only from one sub-lattice A to another sub-lattice B) 
leads to a p-h nested Fermi surface at half filling^. 

• B^P is diverging at a fixed momentum k, if: 

1. The Fermi surface FS has some finite overlap with k — FS 

2. The overlap is such that the occupied (empty) side of FS is covering the occupied (empty) 
side of k - FS. 

This property is called p-p nesting with a nesting vector k. For example, every Fermi surface 
with spatial inversion symmetry is p-p nested by k = 0. 

These are necessary conditions if the density of states at the Fermi level is finite, i.e. in the 
absence of van Hove singularities. 

The situation is particularly rich if the Fermi surface has flat portions. In this case there is 
usually a continuum of nesting vectors k satisfying the above conditions. For example the square 
Fermi surface of the nearest-neighbor tight-binding model is both p-h and p-p nested by every 
vector k || (1, ±1), which is parallel to one of the edges of the square. 



'The reason for this is p-h symmetry as explained in Section 



3.6.1 



Chapter 3 

Renormalization group formalism 



In this chapter I derive the RG equation for the vertex function to the lowest (one-loop) order. As 
opposed to some more modern formulations [ pha94 ^SOOj , |HMOO| , [HonOOt , approach chosen here is 
quite close to the parquet summation technique [EIGN69|, where the aim is a consistent summation 
of the leading divergent terms of the perturbation series (see the expansion (2.35)). We have done 
this in sections 2.3 and 2.4 for the case where the divergences exclusively arise in the p-p bubble. 
In the following it is generalized to the case where divergences appear in the p-h channels as well. 

I feel that the less modern route chosen here makes the meaning of the one-loop approximation 
more transparent. In particular it gives a natural guideline for deciding which terms are to be 
included in the calculation and which ones are to be neglected. 

The chapter starts with the derivation of several exact Bethe-Salpeter equations (Section [3.11 ). 
These are used in Section B.2 to derive the one-loop RG equation for the vertex function. In Section 
3.3, the RG flow of the vertex is related to various generalized susceptibilities associated with pairing, 
density waves and flux phases. In Section 3.4, I make the connection to Wilson's effective action 
and compare different one-loop RG approaches to 2D fermions known in the literature. In Section 
3.5 , I reconsider a general (non-nested) Fermi surface and argue that in this case, the one -loo p RG 
equation is not better than the ladder approximation of Chapter Finally, in Sections 3_^ and 3/7 I 
consider two specific Fermi surfaces, first the case of perfect nesting in an anisotropic tight-binding 
model and second a non-nested case at van Hove filling. The nearest-neighbor tight-binding model, 
where perfect nesting and a van Hove singularity occur simultaneously, will be addressed later in 
Chapter ||. 



3.1 Bethe-Salpeter equations 



If the divergences arise in both p-p and p-h bubbles, the ladder approximation used in section pT3| is 
not sufficient. In fact, in such cases the leading term of the perturbative expansion ( 2.31 ) is given 
by the so-called parquet diagrams. These diagrams are obtained by retaining the five one- loop (or 
second order) diagrams shown in Fig. 2T and by replacing any bare vertex by one of the one-loop 
diagrams and continuing this process to any order. An example is given in Appendix 

To keep track of all the logarithmic divergences, it is useful to define the set of two-particle 
reducible diagrams. The diagrams contributing to the IPI vertex T\(ki, k2, k^) (with an infrared 
cutoff A) can be two-particle reducible in three possible "channels" . A diagram is called reducible 
in the p-p channel if it has the structure shown in Fig. 3.1 a), i.e. if it can be divided in two 



disconnected pieces by cutting two particle lines. Similarly, the diagrams of the form shown in Fig. 
jlb) and c) are called reducible in the channels p-h 1 and p-h 2, respectively. 
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:: 




c) 



Figure 3.1: The structure of diagrams, which are reducible in the p-p channel (a), the p-h 1 channel 
(b) and the p-h 2 channel (c). The grey circles represent any sub-diagram. 



Let us denote by 

R^{ki, k2, fca) the set of reducible diagrams in the p-p channel, with an infrared cutoff A, 



^a'^(^1i ^2: fcs) the set of reducible diagrams in the p-h 1 channel, 



i?^ (fci, fc2, fca) the set of reducible diagrams in the p-h 2 channel, 
If^{ki, fc2, fcs) The set of two particle irreducible diagrams. 

It is rather simple to check that a given diagram is either two-particle irreducible or reducible 
in only one of the three possible channels p-p, p-h 1 and p-h 2. Hence 



rA = /A 



Rf 



R 



phi 



R 



ph2 



(3.1) 



It is also useful to define /' 
o = pp,phl or ph2. 



A 



R^, the set of irreducible graphs in each of the three channels 



Clearly, R'^ is given by an infinite series similar to the ladder summation Eq. (2.43), where — g 



is replaced by /^'^ and the single particle propagators Ca are replaced by full Green's functions Ga- 
We write this series formally as 



nPP _ jPP r-iPP jPP 
^A — ^A ^A ^A 



iPPDPPirDPPiPP + . . 



(3.2) 



The first-order term is missing in Eq. ( |3.2| ) , since it corresponds to an irreducible contribution. The 
series Eq. (3^) satisfies the integral equation 

(3.3) 
(3.4) 



and since T = PP + Rpp 



nPP 



IT DPP T A. 



This is the Bethe-Salpeter equation^ With all the functional dependences included it reads 



RPP{ki,k2,k3) 



p 



/^^(fci, fc2, k - p)DP^j^{p) Tk{p-, k-p^ks), 



(3.5) 



^See for example [AGD61, Noz64|. The Beth e-Salp eter equation was originally introduced to calculate two- particle 



bound-states within quantum electrodynamics [3B51 



3. 1 . BETHE-SALPETER EQUATIONS 
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Figure 3.3: The general structure of a reducible diagram in the p-h 2 channel. The wavy lines drawn 
horizontally stand for —I^^^, but the wavy lines drawn vertically stand for —I^^. The electron 
lines stand for full propagators Ga. 



First, RP is given by a series over 



where fc = fci + fc2 and D^^j,{p) = G\{p)GA{k -p). 

There are similar equations for the two other channels 
the p-h ladders shown in Fig. 3.2, where the wavy lines f 
denote the full propagators Ga. In complete analogy with the p-p case one finds the Bethe-Salpeter 
equation for the p-h 1 channel 



If^^ and the electron lines 



R^^^{ki,k2,k3) = "572 



-pjl ^a''^ (fci ' P + 91 , h)Df^g, (p) Ta (p, k2,p 



-qi), 



(3.6) 



where qi — k^ — ki and 



(3.7) 



The p-h 2 channel is more involved. In fact, a general i?P -diagram has the structure shown 



in Fig. 3.3, where the wavy lines drawn horizontally stand for —I^"^, but the wavy lines drawn 



vertically correspond to — . The electron lines stand for full propagators Ga. Note that the 
diagram must have at least one horizontal wavy line (in the opposite case it is a p-h ladder and not 
reducible in the p-h 2 channel). 



The simplest examples of p-h 2-reducible diagrams are shown in Fig. 2.1. The first PII2 diagram 
of this figure is given by 



(3.8) 
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where (72 = ^ k2. We will write this specific convolution of the functions 
formal product 

rj jph2 T^ph jph2 

Note that this notation differs from the formal product introduced before in the p-p case 
With this notation one can write 

j^ph2 2 jph2 jjph jph2 _j_ jphl^jjph jph2 



ph2 1 T^vh 

^ and D'f^ as a 



(3.9) 



ph2jjph(^j^jphi^ + "higher order terms," (3.10) 

where the overall A-index has been omitted for the simplicity of notation and XF{ki,k2,k3) = 
F{k2, ki, k^) for any function of three energy-momenta. The p-h ladder series can also be written 
using the same notation: 

XRP^^ = {XP^^)DP^{XP^^) + "higher order terms." (3.11) 

The following definitions turn out to be very useful. 



R" = 2RP''^ - XRP^^, 

= r + R'' = {2-x)r 



(3.12) 



and 



R' = -XRP''\ 

P = -XP''^ (3.13) 
r" = + i?" = -XT. 

The superscripts c and s refer to charge and spin, respectively. Note that F^/'' — F^^ ztF^-^, so these 
vertex functions ent er na turally in the calculation of charge- and spin- response functions. 
Eqs. (|to1) and ( ^Tll) can now be written in the simple form 

j^s^c ^ _p,cj^php,c ^ aj^jgi^gj. oj-der terms." (3.14) 

The exact expressions for i?"'^ are given by the infinite series 

j^s.c _ js.c £jph j-s,c js,c jjph js,c jjph js,c _ ^g-j 



Proof of Eq. ( |3.15| ) 

We first write i?^'''^ as an infinite series Rp^^ = Yl^=i^r!^^^ where R^^ is the set of diagrams 
which are exactly n-times reducible in the p-h 2 channel. A similar decomposition is done for R^^^ 
The recursion relations between i?fj'^\ , and R^^^R^^"^ are shown graphically in Figs. 3^ 

and |3J. The relation for R^^"^ (Fig. |j) can be understood as follows. A general i?^ , j^-diagram 
as shown in Fig. |3.3| can either start on the left with a vertical or a horizontal wavy line. If it 



starts with a vertical line, then it is given by the second diagram of Fig. 3.4. If it starts with a 
horizontal line, there are two cases. If there are no other horizontal lines except for the starting 
one, the diagram is given by the third term of Fig. 3.4, If there is more than one horizontal line in 



the diagram, it is given by the first diagram in Fig. 3.3 
Analytically the equations depicted in Figs 



3.4 and 3.5 read 



<+i = -2 jph^DP^'RP''^ + XP'^^DP'^RP^^ + P'^^DP'^XRP^^ 



XRt 



phi 
■n+1 



XP^^DP^Rl^^, 



or, with the definitions (|I|) and ( ^A3|) 
From this, it is easy to deduce Eq. ( 3.15| ) 



(3.16) 
(3.17) 

(3.18) 

□ 



3.2. RG FLOW OF THE VERTEX 
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Figure 3.4: The recursion relation for R^^i. Wavy lines drawn horizontally stand for 1^^^ and 
wavy lines drawn vertically stand for /p''^. The grey rectangles with index n stand for Rfl^'^ {R^^) 
if they are drawn horizontally (vertically). 
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Figure 3.5: The recursion relation for -Rf^'V^j. The conventions are as in Fig. 3.4 



In analogy to the p-p case, we find now the Bethe-Salpeter equations for F* and F'^ 

3.2 RG flow of the vertex 



The Bethe-Salpeter equations ( |3.5[ ) and ( [3.19 ), three integral equations for the three unknown 



functions Rpp , R'^ and R^ , are in general difficult to solve. 

We will address the more modest task of calculating the leading term in the perturbative ex- 



pansion (2.35). This will be accomplished via the flow equation, i. e. a differential equation for 
9a F where we keep only the leading terms in A ^ 0. 

Within this approximation, the two-particle irreducible vertex / is given by the bare interaction 



-g. Thus by Eq. (3.1) 



t{ki,k2, h) = RPP{h,k2,k3) + RP''\ki,k2, fca) + R'"'Hki, k2,k3), (3.20) 
where the dot means a partial derivative with respect to A. 



Consider the first term Rpp. The derivative of the Bethe-Salpeter equation (3^) has three 
contributions 

RPP — jPP £)PPY + JPP£)PPY + pPjjPPY. (3 21) 

Although IPP is by no means negligible by itself (since it contains terms which are reducible in the 
p-h 1 and p-h 2 channels), the first contribution, IppDppT, can be shown to be of subleading order 
in A. I will not prove this here, but an example is discussed in detail in Appendix The last term 
is written as /pp^pp/pp -i- JppdppRpp and IppdppIpp is neglected for the same reason as IppDppT. 
Therefore 

RPP = IPPDPPT + pP DPP RPP. (3.22) 
This equation can be iterated to give 

RPP = jPPjjPPY + JPP DPP IPP DPPT + JPP DPP IPP DPP IPP DPPT + ■ ■ ■ 



28 



CHAPTER 3. RENORMALIZATION GROUP FORMALISM 



fph2 



Figure 3.6: Negligible terms in the flow equation for OaT. The grey circles stand for any sub- 
diagram. 



TDPPT, 



or written out with the full functional dependencies 

R^^{ki,k2,k3) = ^ r A (fci , /C2 , - p) D^A^k ip)'^A{P,k~p,k3), 



(3.23) 



(3.24) 



where fc = fci + A:2 is the total frequency-momentum. 

The same kind of differential equations are obtained in an analogous way for the p-h channel^ 
Using Eqs. ( p^ and ( p^ , we obtain 



^s,c _js,cj-^phY'S,c js,c j~^ph ^s,c 

js,c j-^phY^s,c _|_ js,c j~^ph js,c j-^phY'SjC 

js.c j-^phY'S,c ,c j^phj^s ,c 

The only point to be verified is that the terms /s.c£)p/ips,c js.cj^phjs.c^ .^^j^jj,^^ 
out in Eq. (|^), are in fact of subleading order in A. Writing out /"'^ in terms of 
one can verify that all the neglected terms are of the form shown graphically in Fig. 
be shown to be negligible in analogy with the example of Appendix 
Eq. ( p8| ), expressed in terms of Rp^^ and Rp''^ leads to 



or written out: 
R^'^{ki,k2,k3) 

R^'^{ki,k2,k3) 



Rp^^ = X ((xr)i)f''(xr)) 

j^ph2 ^ _2 TDP'^T + iXT)DP''T + TDp''{XT), 

X! ^A,gi ip) Fa (p, , P + gi )rA (fci , _p + gi , fca) 
Jf2 E ^ti2^P) [-2rA(A:i,p,p + g2)FA(p + 92, ^2, A^a) 



(3.25) 
(3.26) 
(3.27) 
(3.28) 

have been left 
jP^i and iP^"^, 
They can 



(3.29) 
(3.30) 

(3.31) 
(3.32) 



-|-Fa(p, ki,p + (72)Fa(p + 92, fe, fcs) + FA(fci,p,p -I- (72)FA(fc2,P + 92, fcs)] , 



where qi = k^ — ki, q2 = k^ — k2 are the direct and exchanged transfered momenta and D^(p, q) 
Ga(p)Ga(p + 9). 



^Remember that the formal products of functions ninan som nthi riF different in Eq. (3.2^ ) than in Eq. (3.2J). The 
formal product in t he p -p ch anne l is defined by Eqs. ( p.4[ ) and (3.5), whereas the formal product in the p-ncnannels 
is defined by Eqs. (Esh and (Eol). 
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Eqs. ( 3.20 ), ( 3.24 ), ( 3.31 ) and ( 3.32 ) are the one- loop RG equations we were looking for. They 
describe the behavior of the vertex under a differential change of the energy scale A. Because of 
their central role in this thesis, I will put everything together once more: 



RTih,k2,h) 



R''^'^{ki,k2,k3) 



One-loop RG equation 

r(fci, ^2, fcs) = RPP{ki,k2,k3) + 7?f"(fci, fe, fca) + i?'"''(fci, ^2, ka) 

■j^ X! ^Zk(p) ^A{ki,k2, k ~ p)Ta{p, k-p, fcg) 
p 

-Jf^ X! ^a(p, k2,p + <?i)rA(fci ,p + qi,k3) 

X! -^^92 (P) [- 2r A (fci , p, p + 92 )rA (p + g2 , A:2 , fcs ) 



(3.33) 



+Tiy{p, ki,p + q2)TA{p + q2,k2,k3) + TA{ki,p,p + q2)TA{k2,p + q2, k^)] , 



where k = ki + k2 is the total frequency-momentum and qi — k^ ^ ki, q2 — k^ — k2 are the 
direct and exchanged transfered frequency-momenta, respectively. The p-p and p-h propagators 
are defined as D^f^,,{p) = GA{p)GA{k-p) and D^^^^{p) = Ga{p)Ga{p + q). 



This equation has a very simple diagrammatic interpretation. The three terms RPP, RP^^ and 
Jiph2 correspond to the one-loop diagrams shown in Fig. 2.1, where the wavy lines now represent 
full vertices Fa and the product of two single particle propagators has to be derived with respect to 



A. Since D^A^^^ip) — Ga(p)G'a('? ± p) + Ga{p)Ga{<1 ± p), each diagram can be viewed as the sum 
of two terms, where one of the two lines represents the propagator Ga{p) and the other is Ga{p), 
the propagator restricted to the energy scale A. 

Although we have defined A as a sharp infrared cutoff in the band energy, the RG equation 
is also correct for a cutoff in the frequency or for smooth cutoffs. It is straightforward to use the 
finite temperature instead of a cutoff to regularize the theory. The RG equation for Ot^t is then 
obtained by replacing in Eq. (3.33) 



1 



^b^ip.q)... J2 ^ dr (^D''{p,q)\p^^27rT{2n+l) 



(3.34) 



where o = pp,ph. 

Different regularizations of the theory give the same result to leading logarithmic order. For 
instance, the vertex at zero temperature and finite infrared cutoff A is, within logarithmic precision, 
equal to the vertex at temperature T = A without cutoff. The following calculations will be done 
at zero temperature and a finite cutoff, but at the end the parameter A can be interpreted as the 
temperature. 

Since the one-loop RG equation is only correct to leading order in A, only the leading terms 



are to be taken seriously. Taking into account subleading terms in Eq. (3.33) is inconsistent and 



arbitrary, since these subleading terms are comparable to terms that have been neglected in Eq. 

(Hi)- 

This remark about subleading terms concerns in particular the self-energy corrections to the 



propagator Ga- As already pointed out after Eq. (2.43), the self-energy changes the shape and lo- 
cation of the Fermi surface. This effect is by no means negligible in perturbation theory. However, in 
some special situations, that will be addressed later, the Fermi surface is fixed by exact symmetries, 
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so the self-energy cannot affect it. Other effects, such as the renormahzation of the Fermi velocity 
or the reduction of the quasi-particle weight are of subleading order, as will be checked later from 
case to case. So these effects cannot be addressed consistently within the one-loop approximation. 
For this reason, I will replace the dressed electron propagators Ga in Eq. (3.33) by the bare ones 



Eq. (3.33) differs from the one-loop RG equation derived by Salmhofer and Honerkamp [HonOO 



HSFROl, 3H01| only in such self-energy terms, which should be neglected for consistency reasons. 



In fact, if the cutoff is introduced by some function multiplying the bare propagator 



Ca(p) 



then the full propagator is given by 

Ga{p) 

The derivative with respect to A gives 



Oa{p) 
ipo - ' 

Ca{p) 
1-Ca(p)Sa(p)' 
Oa{p) 

ipa - S.P - (^a{p)^a{p) 



Ga{p) 



dA{p){WO - Cp) 



9ib)SA(p) 



(*Po - - 0a(p)Sa(p))' (*Po - ^p - 0a(p)Sa(p))2 



(3.35) 

(3.36) 
(3.37) 

(3.38) 



The first term of the right hand-side is the single-scale propagator of Ref. [ HSFROl |. The one- loop 
RG equation of Salmhofer and Honerkamp is obtained from Eq. (3.33) by omitting the second term 
ofEq. (E3i). 



3.3 Linear response 

To study the instabilities of the electron gas, we investigate the linear response of the system to a 
weak external perturbation. 



3.3.1 Generalized susceptibilities for superconductivity 

Let the system be coupled to a time dependent, real valued field A(r, i), via an additional term in 
the Hamiltonian 

H..t., t^-Yl ^(r' t) (^A(r) + A^(r)) , (3.39) 

r 

where 

^A(r)=E/^('"-^')^r'iCr'T, (3.40) 

r' 

is the locally defined superconducting order parameter. The index A stands for the internal wave 
function — r') of the Cooper pair. This term can also be written in Fourier space as 

k 

where 

A^(k) = ^ e''''-A^(r) = ^ /^(p) Cp^Ck-pt (3.42) 

r p 
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and the functions A(k) and /yi(p) are the Fourier transforms^ of A(r) and /^(r). 

One wants to test the effect of this term on the expectation value (A^(r))t, where the internal 
wave function /s, can in general be different from /a. The coupling of the system to the field A 
induces in general a non-zero value of (A^(r)). Physically, this is a strongly idealized version of the 
superconducting proximity effect. The linear response is given by the retarded response function 
or susceptibility XretBAi^ — r', i — t'), via 

(As(r)),= jdt'Y,X^Z'BA{^-v\t-t')Kv',t') (3.43) 

or, after Fourier transformation in space and time, 

(As(k)). = xfc?lA(k,^) A(k,c.). (3.44) 



The global superconducting order parameter is given by A(k — 0), but later in Section 4.5, I will 
also consider the more exotic case A(k = (tt, tt). 

A more convenient quantity is the Matsubara response function given by 



xiS^(k,^) = ;^^''dre- (A5(k,r)Ayk) 



(3.45) 



where the Bose-Matsubara frequency v is restricted to even multiples of 27r//3 and the imaginary 
time dependence of operators is given by Eq. ( ^.IS]) . The retarded response function is obtained 
by analytic continuation 

x{^)"^=^Xr.t{u:), (3.46) 

where 5 is an infinitesimal quantity, which mimics dissipation. 
From Eq. (3.45) and using Eq. ( 2.24| ), wc calculate 



xlT{k) = rdre-^^/B(p)/A(p') (cpx(r)ck_pT(r)ctp.TC^,.i) (3.47) 

" p.p' 



1 



/b(p)/a(p') G^\p,k~p,k~p'), (3.48) 



(/3L2)2 ^ 



/b(p) [PlH,,,, + Dl^p) rf ^^(p,p')] cr(p') /a(p'), (3.49) 



where k = [v, k), T^'-'^ {p^p') = T[p, k — p, k — p') and /^(p') is the complex conjugate of /a(p')- 

Let us perform the change of variables p — )■ ^^pj^p, where ^ is the band energy and is a 
suitable (angular) variable. In the weak coupling limit the low energy behavior is determined by 
the degrees of freedom close to the Fermi surface. Therefore the dependence of the form factors on 



^ is irrelevant (the argument is analogous to that in Eqs. (2.38) and (2.39)). The most elementary 



superconducting susceptibility is a function of two angular variables, which is obtained by putting 



fsip) - 5{9p - e) and /^(p') = 5(0p' - 9') in Eq. (3.46) 



X^^^(fc,0,0') = /de jA£,'J{^,6)J{^',e')^Y.G^^^Pie,i)^k-pf^e,i).k-Pi^e',i')) (3.50) 

Pa-.p'o 



''Throughout this thesis I use the conventions F{r) = e*'"'F(r) and -F(k) = J^^. e ''"'_F(k) for any quantity 

F(r). The only exceptions are the creation, and annihilation operators. There, I use a pre-factor for both 
transformations. Note that -Ft(k) denotes the Hermitian conjugate of -F(k) rather than the usual Fourier transform 
of Ft(r). The Fourier transform in (real) time is defined as F{t) = / |f e-'"'F(a;) and F{uj) = J di e''^*F{t). 
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where J(^, 6) is the Jacobian, such that L~'^ J2p ~^ I^^ I ^^•^(^7 

These susceptibihty allows for an analysis of the superconducting instability without any prej- 
udices on the form factor /(p). The natural form factor is obtained by writing x(^)^') ™ the 
diagonalized form 

X(f?,^') =Exn/«W/«(^'), (3.51) 

n 

where Xn and fn are, respectively, eigenvalues and eigenfunctions of the operator 

f{e) fde'x{0,o')f{0'). (3.52) 



If Xn diverges at an energy scale Ac, this indicates a (spin-, charge- or superconducting) instability 
with the form factor /„. 

To simplify the calculations, one often chooses by hand a certain form factor /(p) and calculates 
the susceptibility for /^(p) = /b(p) = /(p)• 
3.3.2 Density waves and flux phases 



In the same way as we coupled the system to a particle-particle pair creating operator Eq. (3.40), 
one can perturb the system with a term that creates particle-hole pairs. We again allow for two 
different form factors /a and fs- 



A 



A/B[ 



(r) = ^ E '-fA/B{r' - r)cJ,^Cr„ (3.53) 



r .cr 



or in Fourier space 

^A/B{^) = 2 X! «'^/^/s(P)Cp<TCp+q<x> (3-54) 

where either = cr for the spin susceptibility or So- = 1 for the charge susceptibility x'^- 

What is the physical interpretation of the order parameter AB(q)? Writing this operator in real 
space gives 

^b{ci) = IY. ^^^"""'fsW M'.^ra- (3.55) 

Its physical nature depends on the orbital form factor fs, the spin structure Su and also on the 
momentum q. Different cases will be discussed in this thesis. 

• Spin- and charge-density waves 

The simple choice /b(p) — 1 yields a charge- or spin-density wave if = 1 or cr, respectively. 
Most prominent is the checkerboard case, where q = (tTjTt). In Section 4.5, we will also 
consider the limit q — s- 0, where we obtain the particle number and the total spin operators, 
respectively. The corresponding susceptibilities are the charge compressibility and the usual 
uniform spin susceptibility. 

• Flux phases 

If /s(p) is a function with d2.2_j,2-symmetry, q = (tt, tt) and So- = 1 (c), then the nearest- 
neighbor-terms in Eq. (3.55) yield circular charge (spin) currents flowing around the plaque- 
ttes of the square lattice with alternating directions (see Fig. 3.7). These four charge and 
spin instabilities have been discussed a long time ago in the context of the excitonic insulator 
HR68|. 
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Figure 3.7: The pattern of charge (spin) currents along the bonds of the square lattice in a charge 
(spin) flux-phase 



We call t he phase with circulating charge currents the charge flux-phase (CF) |]AM8g , Kot88 , 
AZHA88 |, it is sometimes also called d-density wave, charge-current wave or orbital antifer- 
romag netism. The ch arge flux- phase (CF), closely related to the concept of the chiral spin 
liquid |WWZ8S|, |LZ8g||, stiU pla ys a prominent role in the strong coupling SU{2) theory of the 
t — J model I WL9(j , LNNW98| . Recently it was proposed to be the competing order param- 
eter to d-wave superconductivity and responsible for the pseudo gap phase of the cuprates 



CLMNOl] 



The phase with circulating spin currents is called the spin flux-phase (SF). Other names 
encountered in the literature are "spin current wave" or "spin nematic state" (because it is 
a state with broken rotational symmetry and unbroken time reversal symmetry). The low- 
temperature thermodynamics of both the charge and spin flux-phases have been investigated 
in mean field approximation in Refs. | NV8E , NJK9l| . 



• Fermi surface deformations 

I will discuss the case q 



also with a d„ 



-v 



-wave form factor. In the 



In Section 4.5, 

case So- = 1, the term i/cxt modifies the dispersion relation et, i.e. it deforms the Fermi 
surface. It makes the electrons move preferably in the x-direction, than in the ^/-direction. 



The nearest-neighbor terms of Eq. (3.55) then measures, the difference — Ty, where (Ty) 
is the kinetic- or hopping energy of the x (y) bonds. In the case Sa- — cr, the Fermi surface 
deformation is in the opposite direction for up- or down-spin electrons. 



The generalized charge- and spin susceptibilities are given by 



Xiiliq) = ^^''dre-^(A5(q,r)(Ali(q) + A^(-q) 
° p,p' 



(3.56) 



(3.57) 



where q — {v, q). We consider in this thesis, either q — > (tt, tt) (Section [I.2| ) or q ^ (Section [4.51) . 
In both cases, 2q 0, and A^/5(q) can be made Hermitian if we choose /(p + q) = /(p). Within 
this choice, the spin- and charge susceptibilities can be written in perfect analogy with Eq. ( ^.49| ) 

xbKi) = Tgl^J2Mp)7I(v^{G^Hp',p + q,p)TG^Hp',p + q,p))- (3.58) 



34 



CHAPTER 3. RENORMALIZATION GROUP FORMALISM 



(/3L2)2 



J2 fsip) \-PL^Sp,p, + DP\p) rfip,p')\ Dl\p') /a(p'), (3.59) 



p,p' 



where Ti{p,p') = -XT{p' ,p + q,p) and V^Jp.p') = {2 - X)T{p' ,p + q,p). 



3.3.3 Flow equations for the susceptibilities 

The RG formahsm can be used to calculate the susceptibilities to leading order in the infrared cutoff 
A. For this, rewrite Eqs. (Ooh and (|3.59D as 



Xba,a(9) 



(3L 



]-,Y.Z^%{p)DZ,(p)U{v) 



where the effective field vertices Z are defined as 



p' 



z-JlM = fBip) - ^Y.fBip')Dt(p)<>'^p)■ 



{3m) 

(3.61) 

(3.62) 
(3.63) 



p' 



The vertex function T^^'^ in Eq. ( |3.62| ) can be expressed by an infinite series in the p-p irreducible 
part IPP. Formally 



Z 



BCS ^ ^ f^j^pp^BCS 

= fB + /bOPPIPp + fj^DPPpPDPPpP 
^ !b + Z^^^DPPpP 



(3.64) 



To get the flow equation, we take the derivative with respect to A of Eq. ( |3.64 ). It gives three 
terms 

Z^^^ = Z^^^DPPpP + Z^^^DPPpP + Z^'^^DPPjPP. (3.65) 



The last term of Eq. ( p.65| ) is neglected for the same reason as in Eqs. (|3.21| ). The remaining 
equation is iterated to give 



^BCS _ ^BCS £)ppjpp _|_ ^BCS £)PP jpp jjpp jpp ^ ^BCSjjppjppjjpppp jjpp jpp . . . 
^ ^BCS ^pp^BCS 



(3.66) 



This is the RG equation for the field- vertex Zg'^'^ . The form factor enters as an initial condition 

zEZM^fBip)- 

It is now easy to obtain a RG equation for x^'^^ as well. From Eqs. ( 3.601 ) and ( 3.66| ) 



■BCS 



zBCS 7 



Z§'-^DPPfA + Z^-u"jA 

zBCS tSpp rBCS r>pp j: , ryBCS A 



ZBObj^pp j.BL.i,j~)Ppj^ _^ Z^'-'=DPPfA 



^BCS r 
^BCSfjPPZBCS 



DPP {fA 



^BCS 



DPPfA 



(3.67) 



In the last equation, we have introduced the function ^^'^■^^.(p), which is the same as Zg'~j^^{p)^ 
except that the initial condition fB{p) is changed into /a(p). 



3.4. RELATION TO THE WILSONIAN APPROACH 

Written out, the RG equations for the pairing susceptibihty are 
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(3.68) 
(3.69) 



where o = A or B, and rf^^{p',p) = T^ip', k~p' p). 



Once the vertex function r-^'-^'^ is known from Eq. ( 3.33| ), one can solve the Hnear equation 
( ^.69| ) for the initial condition Zf^^j^{p) — foip) and finally one can integrate Eq. (3.68) to obtain 
the susceptibility. 

All the steps can be repeated for the charge- and spin susceptibilities, just by replacing the 
quantities z^^'^ ,T^^^ , Ipp by ^^/c ps/c js/c propagator Dpp by -Dp''. 







(3.70) 


KLip) 


- E Zl[l,{p')DtM) Tt.ip'.p), 
p' 


(3.71) 


where o = yl or i?, q{p' tP) = 


-XTA{p',p + q,p) and ri^^ip',p) = (2-X)rA(p',pH 


i-q,p)- 



3.4 Relation to the Wilsonian approach 

A key ingredient to the Wilsonian RG is the idea to replace the given problem by a different one 
with less degrees of freedom but with the same low energy behavior. To achieve this, the effect 
of the eliminated degrees of freedom is incorporated in a renormalization of the parameters of the 
effective low-energy theory. 

This strategy has been successfully followed using Brillouin-Wigner perturbation theory in the 
strong coupling limit of various many-body problems (see for example |Eme79, Aue94, Faz99, 
ABCO(| ). The best known example is the Hubbard model at half- filling which is represented, 
in the limit of strong coupling, by a Heisenberg spin Hamiltonian in the limit of strong coupling. 

In the weak coupling limit, the application of the Brillouin-Wigner formalism to compute an 
effective Hamiltonian is less evident and has not been followed to my knowledge. A tractable imple- 
mentation of the same idea (i.e. elimination of high energy degrees of freedom and renormalization 
of parameters of the effective theory) is given more easily in the functional integral representation, 
which is the content of the following section. 

It should nevertheless be mentioned that the idea of effective Hamiltonians on a reduced Hilbert 
space led to a most po werful n umerical tool for one-dimensional systems: the density matrix renor- 
malization group (see [ NW99|] for a n introd uction). An alternative route to effective Hamilton ians 
was presented recently b y Wegner |Weg94| and applied to the 2D Hubbard model | GKW01 | and 
other correlated systems | Keh99 , KchOl , EIUOl| . 
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3.4.1 Functional integral formulation 

The functional integral formulation of quantum many-particle systems is presented in detail in 
N088|. Here we only mention some of the results. 

In the functional integral formulation, the annihilation- and creation operators c,.^ and cj^. are 
replaced by anti-commuting Grassmann fields 4'{t,t,<t) and '4>{T,T,a}- They depend on the "imaginary 
time" variable r G [0, /?] in addition to the labels of the single-particle states and they satisfy 
anti-periodic boundary conditions in the r variable. 

The Hamilton operator is transformed into the action 



sm = - 



dr ^ir,r,a){dr " ^i)^ir,r,a) + H[t, , (3.72) 



where H[T,'ip] is obtained from the Hamiltonian by the replacements Cra — > V'(r,r,cr) s-^d cj^, 

V'(r,r,<T)- 

The partition function is given by the functional integral 



Z = J P^e*''^', (3.73) 
where is a short-hand notation for J^^ dipx'i-^'x and the Green's functions are 

G"(a;i,...,x„|a;2„,...,x„+i) = ^ / e^^'^'l V,, • • • V,„i^.„+i • • • i^.,„ (3.74) 



= i'^Px^- ■ ■'^x„^x„+i- ■ ■^X2J, (3.75) 

where we have introduced the notation (. . .) for averaging over Grassmann monomials. 
It is useful to perform a Fourier transformation 

^.fc = (/3L2)-i/2 /%re^^«^ Ve'*"-V(.,r,.), (3.76) 



where k = (A:o,k) contains the Matsubara frequency fco as in section 2.1.2 



The Fourier transformed one- and two-particle Green's functions are 

G{k) = -i^jka^jka) (3.77) 

G"^'^ (fci, A:2, fca) = f3L'^ {i>kio-ipk2a'i>k3C7'i^ki+k2-k3a) (3.78) 
After Fourier transformation, the action reads 

S[^P] = Som - W[^], (3.79) 

with 

SoW = J2 Miko - a)V'<rfc, (3.80) 

where <^k = ek — /i and 

W^M^^^ X! 5(ki,k2,k3)^V^fc,^Vife2<T'V'fc3T''0fei+fe2-fc3<T. (3.81) 

ki,k2^k3 a^a' 

The functional integration measure is also readily expressed in terms of the Fourier transformed 
fields = Hfcff dipkad-ipka- 
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It is convenient to introduce the partition function with source term 

= /d/ic[7/']e-'^['^l+(*'''^)+('^'''), (3.82) 



where we used the short-hand notation (x, ip) := X^trfc Xa-k'ipak and the normahzed Gaussian measure 
is defined by 

The connected part of the correlation functions are obtained as functional derivatives [N088| 



{ipl ■ ■ -Ipnlpn+l ■ ■■1p2n)c = 



6^'' log Z[r]] 



Srj2n ■ ■ ■ Srjn+iSrjn ■ ■ ■&n 
where we have written tpi instead of ipkiUi- 



(3.84) 

r;=0 



3.4.2 Low energy effective action and relation to one-particle irreducible 
vertices 



In the spirit of Section 2.2.2 , the bare propagator is now endowed with an infrared cutoff A on the 
band energy CA(fc) = 9(1^1 - A)C(fc). 

One now defines the effective interaction 

WaM = - log / dA^CAMe-^^I^+^l (3.85) 



which depends on a Grassmann field x (iiot to be mistaken with the susceptibilities, for which 
I have chosen the same symbol). Note that the integration with respect to cI^CaH'] 1^ perfectly 
defined, although is not. This can be seen most easily in the expansion of yVA[x] in terms of 
Feynman diagrams. The evaluation of these diagrams involves only Ca and never C^^. Whenever 
CJ^^ appears in an intermediate step of a calculation (see below), it may be regularized by replacing 
the zero in the Heavyside function by an infinitesimal number. 

Wa[x] lias a twofold interpretation. On the one hand, we can restrict the field x to the low 
energy degrees of freedom ii^ka = 0(A — l^kDV'feo-- The object 

Sf[ij<] = (^<,C-V<) - Wa[V'<] (3.86) 

corresponds then to Wilson's effective action, which describes the system in terms of ■0^ only. In 
fact, for observables (or Green's functions) which depend only on the low energy fields, one shows 

(0[V'<]) = ^ Jvib< Jvi^> e'^^^'''^'''^ 0[^<] (3.87) 

and 

Z= fv^< f V^> e^['/'"^^>l = / I?V< e^A*['^'l. (3.89) 



On the other hand, Wa is the generating functional of amputated connected correlation functions 
with infrared cutoff A because of the identity 

logZAW = -(??, Ca??) - Wa[Ca^], (3.90) 

where Za is given by Eq. ( 3.73| ), with C replaced by Ca. 
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Proof 



To prove Eq. ( 3.9Cl| ), we write the Gaussian measure (Eq. ( 3.83 )) explicitly 

then, only in the numerator, we perform a shift of the integration variables ijj ip ~ x 



(3.91) 



-WaM 



Eq. ( ^.90 ) follows by putting x = Ca?? and taking the logarithm. 



□ 



As a consequence of Eq. ( 3.9C ), the quadratic part of Wa is related to the self-energy Sa by 

— WaW 



SXakSXak 



SA(fc) 



(3.92) 



(3.93) 



1-Ca(A:)Sa(A:)' 

where we have used the following identity for the full electron propagator GA{k) = —{ipcrk'^'ak)A = 
CA(fc)(l - CA(fc)i;A(fc))~^ Therefore in the case |^k| < A the right-hand side of Eq. ( 3.93D simply 
becomes T,\{k). 

Similarly the quartic part of Wa is related to the one particle irreducible vertex Fa. In fact, 
differentiating Eq. ( |3.90D we find 

4 

= -(V'(TfciV'cr'fc2'/'o-'fe3'0crfe4)c,A C^"^(fci) (3.94) 



SXakiSXa'k3SX<T'k2^<Tki 



x=o 



i=l 



ih,...,ki) 



f3Vl\tli'^-CAih)^A{k^ 



(3.95) 



In the last line we have used Eqs. ( 2.24 ) and ( p.7i 



The quartic part of Wa is of the same form as Eq. ( 3.81 ) with an effective coupling function 
QAiki, k2, k^) that now depends on the frequencies as well as the momenta. Taking functional 
derivatives of Eq. ( |3.81 ), we find for |^kj < A, 



(1 - S,„,X)gA{ki, k2,k3) = -rX"'(fci, k2, fcg 



(3.96) 



and thus 



gA{ki,k2,kz) = -rA(fci,A:2,^3)- 



(3.97) 



(7a is equal, up to the sign, to a connected amputated correlation function if all j^k; | > A and 
to the IPI vertex in the opposite case, ^a is therefore not continuous at j^kj = A. A formal and 
non-perturbative proof of these relations was given by Morris [Mor94] for a bosonic field theory. The 
derivation given above is perturbative, but a generalization of the non-perturbative proof of Morris 
to fermions appears to be straightforward. Morris has also shown that Sa and Fa are continuous 
at iCkil = A, in contrast to ^a- 
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3.4.3 Other one-loop RG equations 

The effective interaction satisfies the following exact RG equation for 9a Wa = Wa 



^ oXakSxak 



(7,k 



ak 



(3.98) 



This equation was first derived by Polchinski in t he cont e xt of a scalar field theory [ Pol84 |. A 
proof is given in Appendix ^. Zanchi and Schulz |^an96 , ZSOC | proposed to develop Wa up to 
order six in the fermionic variables and to neglect terms of higher order. Terms of order six are 
not present in the original interaction but they are produced by the RG procedure. Their effect 
is then to renormalize the effective coupling function g^. The result is a closed one- loop equation 
for the coupling functio n g\{ ki, . . . , k^) where |^ki| < A. It is identical, within the correspondence 
gA — — Ta, to our Eq. ( |3.33 ) with one difference. In the RG equation of Zanchi and Schulz, the 
vertices on the right hand-side are not evaluated at the scale A, but at a higher scale A, which is 
given by the band energy of the single-particle propagators, i.e. A = Max{|^p|, |^k-p|} in the p-p 
term of Eq. ( 3.33| ), A = Max{|^p|, |^p+qj} in the p-h 1 term and A — Max{|^p|, ICp+qsl} in the p-h 
2 term. Their equation is thus non-local in A. It is a flow equation with memory, i.e. the flow at 
scale A doesn't only depend on the vertex function Fa but also on the history of the flow. 



Since this is not very convenient, it was proposed |IIMOO| to develop Eq. (3.98) into Wick ordered 
polynomials of the fermionic variables instead of monomials as it was done above. Wick ordering 
with respect to the low energy propagator Da = C — Ca results in the same one-loop equation as 
above but now all the couplings are evaluated at the actual RG variable A and — d [Ca(p)Ca((z)] /dA 
has to be replaced by d [I?a(p)£'a('z)] /dA, i.e., one propagator is at the energy A and the energy of 



the second propagator is now restricted to be smaller than A. This is different from our Eq. (3.33), 
where the second propagator is restricted to higher energies. 

We have seen above that the coupling function of the low-energy effective theory equals, up to a 
sign, the IPI vertex Fa. This correspondence can not be generalized to higher order vertices, i.e. it 
would be completely wrong to say that the sixth order term of Wa is related to the IPI part of the 
three particle Green's function and so forth. In addition, the correspondence Eq. ( [3.97 ) relies on 
the choice of the cutoff which ensures that Ca{p) = 0, for |^p| < A. The correspondence no longer 
holds for an alternative scheme, where for example the finite temperature is used to regularize the 
theory instead of the infrared cutoff. 

The general IPI vertices are obtained by performing a Legendre transformation on the functional 
Wa- Wetterich has presented a renormalization group scheme for bosonic field theories, working 
with this Legendre transformed quantity rather that with the effective action defined in Eq. (3.85) 



|Wet93, 


TW94| 


HSFROl 


, |SH01| 



The idea was implemented recently for the many-fermion problem in [HonOO 



(see end of section |3.2D 



The resulting RG equation is identical to Eq. ( 3.33 ) apart from the self-energy 



The RG equations of the three groups [ ^SOOj , |HMOO| , |HSFROl| differ in the treatment of the 
diagrams which give not a leading order contribution in the limit A — s- 0. For example, the (leading 
order) p-p diagram at total momentum fc = features two internal propagators with exactly 
opposite momenta. They have exactly the same band energy. So the non-locality of the Zanchi- 
Schulz equation is not present in this diagram. For the same reason, the energy-constraint of the 
Wick-ordered scheme can not have any effect in this case. So the three RG schemes treat such 
diagrams identically. The same is true for the p-h diagrams in the case of perfect nesting. However, 



as I have argued in section 3.2, only the leading order terms of these RG equations make sense. The 
subleading terms are of the same order as others that have been neglected. In fact, a consistent 
treatment of subleading terms requires going beyond the one-loop approximation. 

The equivalence of the different RG equations to leading order can also be understood in the 
following way. By successively integrating the one-loop RG equation ( p. 33 ) and expressing the 



40 



CHAPTER 3. RENORMALIZATION GROUP FORMALISM 




Figure 3.8: Example of a two loop parquet diagram. 



result in terms of Faq « —g, one obtains the full series of parquet diagrams. However the structure 
ofEq. ( |333| ) introduces a constraint on the energies of internal lines. For example the parquet 
diagram of Fig. is generated by Eq. (3.32) with the following constraint on the propagators 
1,2,3 and 4 

Min{|aU6l}<Min{|e3U^4|}. (3.99) 

Higher order parquet diagrams are generated with similar constraints, i.e. with a certain "ordering" 
of the band energies, when one goes from the inner loops to the exterior loops of a given parquet 
diagram. It can be checked by introducing this constraint into Eq. ( p3.5| ) of Appendix ^ that the 
constraint docs not change the value of the diagram to leading logarithmic order in A. 

The different RG equations ZSOO , HMOO , HSFROlj all generate the whole series of parq uet 
diagrams, but the constraints are different. In the Wick ordered scheme, the constraint (|3.99| ) is 
changed into 

Max{|Ci|,|6|} <Max{|e3|,|C4|}. 
The RG equation of Zanchi and Schulz introduces the most restrictive constraint 

Max{|ei|,|6|}<Min{|e3Ue4|}, 



i.e. both propagators of the inner loop are higher in energy than both propagators of the exterior 
loop. All these constraints are irrelevant for the leading logarithmic order in A. I conclude that the 
RG equation Eq. (|33|) and those of Refs. HZSOOj , |HMOO| , |HSFROl| are all equivalent to leading 
order in A. 



3.5 The case of a general Fermi surface revised 

In the following, we identify the leading contributions to the one-loop RG equation in the limit of 
small energies. The main part of the work on the half-filled lattice with nearest-neighbor hopping 
will be treated in the next chapter. Here I shortly reconsider the case of a general non-nested Fermi 



surface which has been treated before in Sections 2.3 and 2.4 



3.5.1 Low-energy scattering processes 

First one identifies the possible scattering processes which connect four momenta on the Fermi sur- 
face because only they will be relevant to leading order. These processes have to satisfy momentum 
conservation, i.e. we are looking for solutions of the equation 



ki -t- k2 = ka + k4. 



(3.100) 
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where ki, . . . , k4 are momenta on the Fermi surface (^ki = 0). Note that this equation is understood 
modulo a reciprocal lattice vector 27rZ^. 

One obvious solution is that incoming and outgoing momenta are equal. The corresponding 
scattering processes are called forward and exchange scattering, respectively with vertices and 
T'^ defined by r-''(k,k') = Xr<=(k,k') = r(k,k',k'). To simplify notation I omit the subscript A 
in Fa from now on. A second general solution corresponds to scattering of particle pairs with zero 
total momentum. The corresponding vertex is F^'^'^(k, k') = F(k, — k, — k'). 

If the Fermi surface is big enough, there are additional Umklapp processes, where the momentum 
conservation is violated by an integer multiple of 27r in one or even in both directions. A simple 
geometric criterion for the existence of low-energy Umklapp processes for the case of a closed and 



convex Fermi surface is as follows [KR97]. Let FV = {p;Cp < 0} be the Fermi volume. Given two 
momenta k and k' on the Fermi surface, then the center of mass ^(k + k') is in the Fermi volume, 
due to convexity. Inversely, for every point p of the Fermi volume there exists at least one pair of 
Fermi momenta k and k' such that ^(k + k') — p. Low-energy Umklapp processes exist, if the 
Fermi volume has some non-zero intersection with its own translation G/2 + FV, where G/2 is half 
of a reciprocal lattice vector (i.e. (7r,0), (0,7r) or (7r,7r)). Given two Fermi momenta ki,k2 such 
that the center of mass |(ki -|- k2) is in the intersection FV n (G/2 -I- FV), then there exists a pair 
of Fermi momenta k3,k4, such that ki -t- k2 = ks + k4 + G. Umklapp processes are continuous 
two-parameter families like forward-, exchange- or BCS-scattering, not restricted to certain spots 
on the Fermi surface, nor are they related to some "Umklapp surface", as claimed in |HSFR01|^. 



3.5.2 Analysis of the bubbles 

In order to identify the leading contributions to Eq. ( |3.33 ), we neglect for the moment the angular 
dependence of the vertex. The three different contributions p-p, p-h 1 and p-h 2 arc proportional 
to 

1 >p d[CA(p)CA(fcTp)] 

/3L2 ^ dA 



(3.101) 



where k = ki + k2 in the p-p contribution, k = k^ — ki in the p-h 1 contribution and k = k^ ~ k2 
in the p-h 2 contribution. 

We consider the thermodynamic limit and zero temperature and therefore replace by 
/ (2tt)'^+'^ calculations. Taking explicitly the derivative with respect to A and integrating over 

the frequency po (for fco = 0) we find 



_fp_ 

(2^)2 



e(|gq|-A)e(±gpgq) 



(3.102) 



where q = k ^ p. 

Within the Wick ordered scheme |HMOO|, the first step function in Eq. (3.102) would be replaced 
by 6(A — l^ql), since the second propagator in Eq. (3.101) is restricted to be in the low energy 



part. I have verified t hat th is alternative scheme would not change the final results, confirming the 



conclusion of Section 3.4.3 



''In |HSFR01|, the dispersion relation = —2t{cospx + cospy) + 4t' cos cos Py — fi oi a generalized tight- 
binding model was considered. In the parameter range < t' < 0.5 and fi > — 4t', the Fermi surface consists of 
four disconnected arcs in the Brillouin zone BZ =] — 7r,7r]x] — 7r,7r]. If we choose on the other hand the region 
BZ' = [0, 27r[x [0, 27r[ to represent the momenta p £ R'^/27rZ^, then the Fermi surface is closed and convex and the 
abovementioned criterion applies. But attention: the Brillouin zone BZ' is equivalent to the Brillouin zone BZ only 
modulo 2ttZ^. Thus some Umklapp processes in BZ' correspond to non-Umklapp processes in BZ. But the lattice 
model does not really care about which processes are Umklapp and which arc not. BZ and BZ' are mathematically 
equivalent 
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For general values of k one finds that both bubbles i?^^'^''(A, k) are proportional to log A. For 
the special value |k| — 2, one finds B^^ ^ B'p^ ~ A~^/^. The strongest divergence comes from 
the p-p bubble at k = 0, namely Bpp{A, 0) ~ A^^. Correspondingly, the dominant contribution in 



the low-energy regime A ^ comes from the p-p channel. As it was argued before in Section 2.3, 
the dependence of F^*-^"^ on the frequencies ko, k'^ and the band energies ^k,'Ck' is irrelevant so that 
F^'^'^ can be treated as a function of two (angular) variables instead of six, a result that is already 



known from the more standard scaling analysis used in [Pol92, 3ha94, CFS96]. We therefore write 
pBCS^]^^ k') instead of T^'~'^{k, k'), because the frequency dependence has been neglected. 



Does it make any formal sense to take into account the non-leading p-h terms of Eq. (3.33) 
in the absence of nesting? One has to compare the non leading term s ~ F^ logA with other 



contributions, that have been neglecte d in t he derivation of Eq. ( |3.33 ). According to a careful 



analysis by Sal mhofer and Honerkamp |^H01 |, the first neglected terms are or the order ~ F^ log^A 



(see Eq. 108 of | SH01 |). The contributions to the exact RG flow can be classified into three classes: 
1. leading terms of the one-loop equation, 2. subleading terms of the one-loop equation and 3. 
terms which are neglected in the one-loop equation. Schematically, 

Fa — "leading terms" + "subleading terms" + "neglected terms" . (3.103) 
~riA-i -rJiogA -riiog^A 

The RG flow is split into three different energy regimes. 

1. High energies, in which A is not small. In this regime, the subleading terms of the one -loop 
equation dominate the neglected two-loop terms, provided Fa logA is small. Thus Eq. ( 3.33|) 



can in principle be used to to calculate Fa, starting with a cutoff equal to the bandwidth, 
where the vertex is given by the bare coupling. In this regime, there is no small parameter 



to further simplify the functional RG equation (3.33), so using it is technically very difficult. 
But in the high energy regime (glogA <ti 1) naive perturbation theory provides a controlled 
and much more feasible method of calculation than the RG. 

Small energies, where A is small such that (AlogA)~^ 3> Fa logA > 1 . Na ive perturbation 
theory breaks down in this regime, but the leading terms of equation ( 3.33| ) are superior to 



the neglected terms even if F is not necessarily small. Thi s is w here the one-loop RG equation 
is most useful. Note, that the subleading terms of Eq. ( [3.33| ) are no longer superior to the 
neglected terms in this regime. They have to be neglected to be consistent. 

3. The critical regime, where A is close to the critical energy scale Ac at which Fa is diverg- 
ing. If A is too close to Ac, the neglected terms are no longer negligible and the one-loop 
approximation is no longer accurate. Note however that the weaker the initial interaction is, 
the more A can approach Ac before the one-loop approximation breaks down. 

It is worthwhile to discuss the behavior of the bubbles for small but finite values of k. To be 
specific, I consider the example of very low filling, where the single-electron spectrum is approxi- 
mately parabolic Cp = ~ 1 (^U energies are given in units of the Fermi energy and all momenta 
in units of kp). The energy shell |^p| = A consists of two circles with radius vT±A. For the p-p 
bubble one finds in the limit |k|, A <^ 1 



-1 2arctan 



if Ikl < A, 



BPP(A,k) = <^ }LA±^ II-' /3^^Q4) 

2^V|k^-A^| \ log l'-l+vf^ if|k|>A. 



If we renormalize a vertex with a small total momentum k = ki -I- k2 the p-p contribution to its 
flow is approximately independent of k as long as A ^ |k|. In this case we replace FA(fci, fe, fca) 
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by r^^^{ki, — fca) even if the total momentum is not exactly zero. This replacement is no longer 
justified when A is of order |k|. The flow then depends strongly on k and cannot be controlled. 
Nevertheless, there is no danger because a few renormalization steps later, if A <C |k|, the flow is 
suppressed and the coupling under consideration no longer contributes. 

We can do the same type of analysis for the p-h bubble with a small momentum transfer, 
which renormalizes couplings that are close to the forward- or exchange scattering (^^^(k, k') := 
Xgl{k,k') := .gA(k,k',k',k). We obtain 

f if |k| < A, 

B-(A.« = | ^,<„.M±4!3 if|k|>A. 

It gives a big k-dependent contribution if k is of order A. But this flow is again suppressed if A is 
further reduced. 

The results presented above for the isotropic case should remain valid as long as the Fermi 
surface is both far away from van Hove singularities and not nested. The presence of Umklapp 
scattering does not change the result, since the contribution of Umklapp processes to the flow is 
not of leading order. 

It is not excluded that Umklapp processes influence the flow in an important way before one 



enters the low-energy regime |EISFR01|. But these are most likely non perturbative effects and I 



doubt that they can reliably be accounted for within a one-loop approach. 



3.6 Nesting 



I now consider the case of the half-fllled — ty model, presented in Section 2.5 



3.6.1 Particle-hole symmetry 



It is interesting t o not e that the nesting property ( 2.69 ) and thus also the relation of the p-p and 
p-h bubbles Eq. ( 2.70 ) are consequences of the particle- hole symmetry of the model. 

For a general single-band model on the square lattice, where the hopping is restricted to nearest 
neighbors, the one-body Hamiltonian Hq is invariant under the particle- hole transformation 



c. 



k+QcT- 



(3.106) 



The interaction given by Eq. (2.4) does in general not satisfy this symmetry. However the term 

Hi = Hi- H' (3.107) 

where 

^ " ^ " (3.108) 



^' = 1^2(72 E(2 - P, P) J n.k, 

(T,k \ p / 



is particle-hole symmetric, provided 

5(ki , k2, kg) = .g(ki -f Q, k2 + Q, kg Q). 



(3.109) 



For example, the nearest-neighbor interaction Eq. ( ^.14 ) satisfles the condition ( 3.109[ ), except 
for the Coulomb assisted hopping term K. li K = 0, the term H' is proportional to the particle 
number and can be absorbed in the chemical potential. 
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Figure 3.10: Examples of scattering processes T^{k,k') for different choices of k and k'. The 
"center of mass" of the two incoming and outgoing particles which are are fixed to (±7r/2, ±7r/2), 
are indicated by points. 



The transformation ( ^.106 ) interchanges occupied and empty single-particle states. The particle 
hole symmetry therefore relates the system at particle density n to the system at density 2 — n. At 
half filling (n = 1), the particle hole symmetry implies the following exact identities 

G(fco,k) = -G(-fco,k + Q). (3.110) 

and (using time reversal invariance, spin rotation invariance and parity) 

T{h,k2,k3) ^TiQ ~ h,Q - k2,Q ~ ks). (3.111) 

As a consequence of Eq. ( ^.110 ) , self energy corrections might change the details of the Ferm i 



surface, but not t he n esting property. In fact the dispersion Cp = + S(0,p) satisfies Eq. ( 2.69 ) 



and the relation ( 2.70| ) remains true if one replaces the bare propagators C in the definition of the 



bubbles by exact Green's functions G. 
3.6.2 Low-energy scattering processes 

While the generic forward, exchange and BCS scattering processes are connected to zero momentum 
transfer and zero total momentum, the peculiar (particle-hole symmetric) geometry of the present 
Fermi surface allows for low-energy scattering processes associated with the momentum Q = (tt, tt). 
The processes with direct or exchanged momentum transfer Q are r'^(k, k') = Xr^(k, k') = 
r(k, Q + k', k') and those with a total momentum Q are r''(k, k') — r(k, Q — k, Q — k'). These are 
low-energy scattering processes for k and k' chosen freely on the Fermi surface, since ^Q±k ~ ^S,k- 



Some examples are shown in Figs. |3.9| and 3.1C 



This list of the possible low energy processes is complete but the classification into F-'', F*^, r^'-^'^. 



(defined in Section 3.5.1) and F , F^, F'' is not unique. In fact, two two-parameter families intersect 
in a one-parameter family as T^^'^i'k, k + Q)^ F''(k, Q - k), F^<=''5'(k, Q - k) = F^(k, Q - k), etc. 
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3.6.3 RG flow in the p-p and p-h channels 



As a consequence of the particle- hole symmetry, the p-p and p-h bubbles are related by Eq. ( 2.70 ). 
Thus in addition to T^'~'^ the vertices V and have to be renormalized. 

However the renormalization of a general BCS vertex is the same as without nesting. For 
example we have 



BCS 



where 



r 



(k,k') 



R 



phi 



R 

r/ph2 



ph2 



BP^{K, k + k') and R^^^ ~ ^^''(A, k - k'). We may neglect 



BCS 



RP^^ and RP^'^ for small values of the cutoff, except when k± k' = Q -I- 0(A). Similarly, for generic 
k and k', r^(k, k^) i s renormalized exclusively by the p-h 1 term and r''(k, k') only by p-h 2. The 

and r'^ = 2T'^ - have a 

(3.112) 
(3.113) 



flow equations ( |3.33 ) written in terms of the functions T 
simple form. They read 



r^^^(k,k') i?fp(k,-k,-k') 



f"/"(k,k') 



7?^/^(k, k'-|-Q,k') 



The notation V"/" is the same as in Section |3.3| , except that the transfered momentum is understood 
to be Q. The vertices F'*/'^ are thus relevant for spin- and charge susceptibilities of momentum Q. 
Because of the similarity of Eq. ( 3.28|) with Eq. ( 3.23 ), one obtains three closed RG equations 



(k,k') 



BCS 



r/^(k,k') = 



I3L^ 
1 



(k,p)i?r.o(p)r^^^(p,k') 



p 

p 



r/^(k,p)l)P^(O Q)(p)r/^(p,k'). 



(3.114) 
(3.115) 



Note that we have systematically replaced r(fc, k') by r(k, k'), be cause the d ependence on the 
frequencies is irrelevant. Due to particle-hole symmetry, Eqs. ( |3.114 ) and ( 3.115| ) have the identical 

form 

f-(k,k') = ^^r(k,p)i)^^„(p)r-(p,k'), (3.116) 

where o stands for BCS, s or c. After integrating over the frequency pq, one finally obtains 

-2Ar(k,k') - -l^^dCpl - A)r-(k,p)r-(p,k'). (3.117) 



In spite of the apparent decoupling of r^'^'^, and F'^, there is in principle some coupling by 
hybrid vertices such as F(p, — p, — p -I- Q). However the weight of such terms in the sum on the 
right hand-side of Eq. ( 3.117 ) is negligible. For example, the "mixing" vertices F-^ '-^^(k, p), where 
p = Q ± k + 0(A), give merely a 0(A) contribution to the right hand side of Eq. (3.117) (see also 



Sha94|). The situation is more difficult if the Fermi surface contains a van Hove singularity, i.e. in 
the case tx = ty considered in Chapter ^. The non-uniform density of states increases the weight of 



certain points in the p-integration of Eq. (3.117) and thus weakens the phase space argument used 
here. 

We conclude that if ^ ty, the charge, spin and superconducting instabilities do not influence 
each other in leading order. Moreover the charge and spin instabilities can be studied in exactly 



the same manner as it was done for superconductivity in Section 2.4 



For the repulsive Hubbard interaction, the dominant instability will be a spin-density wave, 
since F^^ = — F^^ « [/ is positive, whereas F^^ and F^'^'^ are negative. One expects thus a Slater 
insulator with a spin-density wave at weak coupling and, as U is increased, a crossover to the 
antiferromagnetic Mott insulator. 
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Figure 3.11: The Fermi surface at van Hove filling: t' = 0.3 (solid line) and t' = (dashed line). 



3.7 Van Hove singularities 

We now consider the case where the Fermi surface passes through a van Hove singularity, i.e. 
a saddle point of the dispersion relation. The density of states is logarithmically diverging at 
the saddle points and leads to a more singular behavior of the p-p bubble at zero momentum. 
Differentiating Eq. ( 2.41 ) with respect to A, one gets 



^PP(A,0) 



KA) 



-A) 



2A 



(3.118) 



Assuming a logarithmic behavior of the density of states, it follows that i3''P(A, 0) ^ —A ^logA 
and thus Bpp(A, 0) - log^ A. 

To be specific, consider the dispersion relation of a generalized tight-binding model: — 
— 2(cosfca; + cosfcy) + 4t'(cos/c2, cos/c^ + 1). The unit of energy is given by the hopping amplitude 
between nearest neighbors. A finite electron hopping < < 1/2 between next nearest neighbors 
has been included and the chemical potential is fine-tuned such that the Fermi surface contains the 
saddle points at Pi = (±7r,0) and P2 = (0,±7r) (see Fig. fj.ll ). 

If t' ^ 0, the leading (log^ A) terms arise only in the p-p channel at zero momentum. In addition, 
there are divergences in i3PP(Q), BP^{Ci) and in BPP{q) for small q, but these divergences are only 
~ log A and are thus not of the leading logarithmic order. If one restricts the RG flow to the 
leading terms, only the BCS coupling function r-^'-^'^(k, k') is renormalized and the only possible 
instabilities are of the superconducting type. The flow of r^'^'^(k, k') is similar to the one explained 
in detail in the next section, but without competing charge and spin instabilities. 

For repulsive interactions, where the BCS flow is towards weak coupling, the approximation is 
clearly not sufficient and non-leading terms might play an important role. Many efforts to include 
non-leading terms in the RG flow have been made [ |LMP87| , |Dzy96| , |FRS98| , |IKK01| , |IK01| , [HSOlbt , 
but a consistent treatment has not yet been attained. In the following, I mention some of the 
problems which would have to be solved to obtain a consistent treatment of the first sub-leading 
terms (i.e. terms of the order A~^). 

Let us first study the behavior of the p-p bubble at small but finite |k|. It is instructive to focus 
on the contributions of a small patch surrounding a saddle point, say the region 



Pi = {p ; VT^2F|p, 



2f\py\<2p} 



The parameter p is small enough so that can be replaced by its limiting quadratic form close 
to the saddle point Pi = (tt, 0). I have computed the bubbles Bp^(A,k) and i3p^(A,k) defined 



by Eq. (3.102), with the summation restricted to p e Pi- The values of i?p^(A,k) and i3p^(A,k) 
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depend sensitively on CPi+k- Both bubbles are negligible if |CPi+k| ^ A, but for |^Pi+k| ^ A (and 
A, < p^) we get 



-1 



log 



4Ap2 



(27r) VI - 4i'2 A " (A + pfc+)(A + p/c_)" 



and 



S?.'(A,k) 



(27r) Vl - 4t'2 A 



e(/9fc+ - A) 



pk+- A 
pk+ 



+ e(pfc_ - A) 



- A 



(3.119) 



(3.120) 



where fc± := |Vl - 2i' ± y/1 + 2t' ky\. Note that k+ fc_ = |^Pi+k|. 

The behavior is most interesting in the regime k^k^ <C A <C pk±. The p-h bubble, zero for 
sufficiently small |k|, gives a contribution ^ A^^ there. Thus for a given small momentum transfer 
k, the p-h contributions are considerable over many RG iterations, in contrast to the case without 
van Hove singularities. Even more strikingly, the first subleading (i.e. A~^) part of the p-p 
diagram is very sensitive with respect to both the size and the direction of the total momentum 
k in this regime. Thus for a given small momentum k, the p-p contribution has a complicated 



behavior on k over many RG iterations, in contrast to the case of Sections 3.5.2 without van Hove 
singularities. 

We conclude that the RG equation to order A^^ for the vertex rA(ki,k2,k3) depends in a 
delicate and non-trivial manner on the momenta. 

In contrast to the case of the circular Fermi surface, self-energy corrections are not negligible in 
the present situation. In fact, they are expected to change the shape of the Fermi surface. 

For the corrections to the single particle dispersion and the qua si particle weight, one finds in 



second order perturbation theory (see for example Eq. 7 of [ Dzy96 



VpEA(p) 



3' log' A 
g^log^A • VpCp. 



(3.121) 



These corrections are negligible within the leading-order approximation, where we assume g log A ^ 
1. An attractive coupling diverges at this scale. However, in the case of repulsive interactions the 
RG flow can be followed to smaller energies such that glogA ~ 1. At these energy scales, self 
energy corrections have to be taken seriously. In other words, a consistent calculation to the order 
A~^ requires to follow simultaneously the RG flow of the self energy and of the vertex. 

Due to the additional logarithm in the p-p bubble the theory is not renormalizable in the usual 
sense of field theory, i. e. it is not possible to send the bare momentum cutoff to infinity while 
keeping some physical correlation functions finite (even after the introduction of counter-terms in the 
microscopic Hamiltonian and of wave- function renormalization) . Gonzalez, Guinea and Vozmediano 
| GGV96| , |GGV97b| , pGV99| , |GGVOC| ] proposed a field theoretical RG scheme where the coupling 
constants for forward and exchange scattering are renormalized only by the p-h diagrams. In this 
approach the p-p channel is treated separately in connection with a renormalized chemical potential. 
It is however not clear from our approach that the RG equations will not mix p-p and p-h diagrams, 
if all contributions ^ A^^ are taken into account. 



CHAPTER 3. RENORMALIZATION GROUP FORMALISM 



Chapter 4 



Instabilities of the half-filled lattice 
with nearest-neighbor hopping 

4.1 Motivation 

This chapter is dedicated to the nearest-neighbor tight-binding model (i.e. — ty — 1 and t' — 0) 



at half filling where the Fermi surface is a square (see Fig. 3.11). This special case has common 
features with the two cases considered before. There is exact nesting associated to particle-hole 
symmetry and the Fermi level is at a van Hove singularity. This case is more difficult, than the 
anisotropic model of Section R^, b ut it leads to nontrivial results already to leading (one-loop) 
order, unlike the case of Sectioii [3.7| without nesting. 

As a motivation, I would like to relate the subsequent analysis to the partially numerical in- 



vestigation of the Hubbard model by Zanchi and Schulz |ZSOC]. According to them, close to half 
filling the flow RG separates into two regimes. For not too low energies A > /i, where /i is the 
chemical potential, the RG flow is not yet susceptible of the deviation from half-flUing (see Fig. 



4.1). Consequently, the RG flow is complicated by the simultaneous appearance of divergences in 
the p-p and p-h channels and by the presence of van Hove singularities. Zanchi and Schulz call 
this the parquet regime, since the interplay of fluctuations of different kinds is the essence of the 
parquet method. 

For lower energies A < /i, the distance of the Fermi level from the van Hove singularity as well 



as the absence of nesting becomes relevant (Fig. 4.2). There, the RG flow is entirely in the p-p 
channel as in Chapter | This is called the BCS regime. 

Zanchi and Schulz have found that the critical energy scale Ac, where the RG flow diverges, 
depends strongly on the electron density (i.e. on the chemical potential /i). It is largest at half 
filling and becomes ex pone ntially small upon doping. One can now draw the schematic phase 



diagram shown in Fig. O. The instabilities at Ac are qualitatively different in the two regimes. 
In the BCS regime, the instability is entirely in the p-p channel and the diverging quantity is the 
vertex function r^'-''^, i.e. the scattering of particle pairs with a vanishing total momentum. As 
a consequence, the superconducting susceptibilities diverge, but not the density-wave or flux-phase 
susceptibilities which depend on nesting. The interpretation of this instability as the onset of 
superconductivity is therefore quite obvious, although there is no proof. 

The situation is much less clear if the instability takes place in the parquet regime. There the 
numerical investigations find many diverging vertices of different kinds and simultaneously diverging 
susceptibilities of the SDW- and of the superconducting type, although the SDW susceptibilities 
are clearly dominating at half filling. How should one interpret this kind of instability? 
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-3-2-10 1 2 3 



Figure 4.1: The Fermi surface of the nearest-neighbor tight-binding model shghtly below half- filling. 
The grey region corresponds to high energies, |^k| > A. The case shown here is in the parquet regime 
A > /X, where RG equations are not yet susceptible for the doping away from half-filling. 




Figure 4.2: Same as Fig. 4.1, but for a lower cutoff A < /i. This is the BCS regime, where the 
saddle points have already been integrated out and the curvature of the Fermi surface becomes 
relevant. 
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I consider the half-filled case as an important key to a better understanding of the instability 
in the parquet regime because it allows us to take the limit of weak interactions and thus very low 
energies A without leaving ultimately the parquet regime. 

The remaining part of the chapter is organized as follows. In Section the renormalized 



couplings are classified according to both the location of momenta with respect to the van Hove 
points and the channels characterizing the different instabilities. It is argued that to leading order 
there is no mixing between superconducting, charge and spin instabilities except from momenta 
very close to the van Hove points. A simple way of disentangling this s pecia l behavior at the 
van Hove points and the generic behavior elsewhere is presented in Section 4.2 and contrasted to 
an earlier approach where the momentum dependence was altogether neglected. The asymptotic 
behavior of the RG flow allows to draw a phase diagram including superconductivity, density waves 
and flux phases, depending on the values of the bare couplings. This phase diagram agrees with 
symmetry considerations linking the various order parameters, as shown in Section 4.4. The content 
of of Sections 4.2-4.4 was published in [ BBD02( |. Finally in Section 4.5, I investigate the behavior 
of some additional physical quantities in the vicinity of each instability. These quantities are the 
uniform spin susceptibility, the charge compressibility, but also more exotic susceptibilities which 
measure the tendency towards ry-pairing (pairing of Cooper pairs with momentum {tt,tt)) or the 
tendency towards a deformation of the Fermi surface (Pomeranchuk instability). 



4.2 RG flow of general vertices and susceptibilities 

4.2.1 Low energy scattering processes 

I flrst identify the possible scattering processes which connect four momenta on the Fermi surface 
and satisfy momentum conservation modulo a reciprocal lattice vector. In addition to the six classes 



of processes discussed in Section 3.6, there are now further processes, which are related to the fact 
that the Fermi surface consists of straight lines. If three points ki,k2,k3 are chosen freely on two 
parallel sides of the square, the resulting k4 = ki + k2 — ks lies automatically on the Fermi surface. 



giving rise to a three-parameter family r" (ki, k2, ka). Some examples are shown in Fig. 4.4 



As in Section 3.6, the list of the possible low energy processes is complete but the classiflcation 
into r^, F'^, T^'-^^/T'^, F^, F'' and F^ is not unique. For example if k and k' belong to the same 
pair of parallel sides of the square, the two-parameter families r-^(k, k'), F'^(k, k'), F-^'-^'^(k, k'). 



52 



CHAPTER 4. INSTABILITIES OF THE HALF-FILLED LATTICE. 



rii(ki,k2,k3) 




Figure 4.4: Examples of processes described by r". The momenta ki, k2 and ks can be chosen 
freely on two parallel sides of the square Fermi surface. 



r''(k, k'), r^(k, k') and r''(k, k') belong to the larger three parameter family r". Furthermore, 
two two-parameter families intersect in a one- parameter family as T^^'^ik, k + Q) = r''(k, Q - k), 
pBCS^j^^ Q — k) = r^(k, Q — k), etc., where Q — (tt, tt). Finally three two-parameter families can 
intersect in scatterings between the two saddle points as r^'^'^(Pi, P2) = r''(Pi, P2) = r^(Pi, P2). 
Where Pi — (tt, 0) and P2 = (0, tt) are the two saddle points. 



4.2.2 Self-energy effects 



As it has been explained in Section 3.6.1 , the nearest-neighbor tight-binding model is particle- 
hole symmetric. We further assume that the interaction also respects this symmetry. The self 
energy then satisfies the exact relation I]a(po,p) ~ ^Sa(— Po,P + Q)- Imposing in addition the 
point symmetries of the square lattice, one can conclude that the self-energy vanishes on the Fermi 
surface {po = Cp = 0)- This means that the Fermi surface is not modified by self-energy effects. 

For the corrections to the si ngle particle dispersion and the quasi particle weight, one can derive 
relations similar to Eq. (3.121), namely (in second order perturbation theory) 



dp„T,A{p) 



VpI]A(p) 



< const. • 5^1og^A (4.1) 

< const. • 5^1og^A • Vpfp. (4.2) 

The calculations are given in Appendix In t he absence of p erfect nesting (van Hove filling, but 
t' ^ 0), it is known that the bounds Eq . (|4.2| ) and Eq. ( [l.l| ) overestimate the self-energy terms 
by one power of the logarithm (see Eq. ( 3.121 )). I do not know whether this is also true if t' = 0. 



In any case the bound given above is sufficient to exclude a leading order correction ~ g log A. I 
conclude that self-energy corrections are negligible within our approximation, where glog^ A ~ 1. 

In spite of the fact that self-energy effects are negligible within the leading order (or one-loop) 
approximation, it cannot be excluded that they influence the physics considerably. In fact, due 
to properties of the self-energy, the half-filled 2D Hubbard model is not a Fermi liquid above the 
critical temperature for the SDW |Riv02|, in contrast to the 2D jellium model which is a Fermi 
liquid above the superconducting critical temperature |DROO]. But quite generally, a consistent 
treatment of self-energy effects requires to go beyond the one-loop approximation, i.e. it would be 
inconsistent to include self-energy effects at our level of approximation. 



4.2.3 RG flow of general vertices 



We now investigate the one-loop corrections to the coupling constants by an analysis of the bubbles 
(3.102). As a consequence of particle-hole symmetry, the p-p and p-h bubbles are related by Eq. 
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Figure 4.5: Plot of the p-p bubble for k = (k, k), parallel to the square Fermi surface and A = 0.01. 
The small mark on the vertical axes indicates the maximal value —Bpp{A, 0,0). 



Both the p-p and p-h bubbles are of order whenever k = n(7r, tt) + k(1, ±1) (n G Z, k G M). 
The reason is that the Fermi surface consists of straight lines. If k is parallel to such a line, it satisfies 
the conditions of Section both for p-p and p-h nesting. It follows that the three parameter vertex 
funct ion is renormalized by contributions of order A^^ from every term Rpp , RP'^^ , and RP'^^ in 
Eq. (HH). 

In Fig. 4.5 we show a plot of i3PP(A,k) for k = {k,k). The analytic expression is given in 
Appendix |e . Here I only mention the asymptotic result for small k 



BPP{A,k,k) 



1 



log 



16 



27r2A ^A + 4:K 



(4.3) 



By restricting ourselves to the logarithmically dominant terms of order A~^ log A we can go a 
long way using an analytical approach, as will be shown now. We consider the two-parameter vertex 
functions r^'~^^ , F^ and F'* and use the same argumentation as in Section 3.6.3. For generic k and 
k', the vertex F^'-^'^(k, k') is renormalized only by the p-p term, r^(k, k') only by the p-h 1 term 
and F'^(k, k') only by p-h 2. 

As a consequence, one obtains the RG equations for the fimctions r^'^'^, F" = — F^ and F"^ = 
2F'^ - F^ as in Section Isjll 



2A— F*(k,k') 



(4.4) 



where o stands for BCS, s or c. 



Eq. (4.4) is valid for momenta k and k', such that k±k' is not too close to Q. It was argued in 
Section 3.6.3, that the contribution of these special momentum configurations to the right hand-side 
of the RG equation is negligible. We will see shortly that this is no longer true here because of the 
non-uniform density of states in momentum space. 



A large contribution to the integration over p in Eq. (4.4) comes from a small neighborhood 
of the saddle-points, due to the diverging density of states. To identify the logarithmically di- 
verging contribution to Eq. (4.4) we consider the patches Pi = {p ; \px — 7r| + \py \ < 2p} and 
P2 = {p ; I + |Py — Ti" I < 2p} of a sizc p <C 7r/2 around the two van Hove points and separate the 
integral into X^pePi +SpeP2 +Spes Z-P1-P2' where B.Z. is the whole Brillouin zone (see 
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Figure 4.6: The Brillouin zone with the energy shell |^p| = A. It is separated into two patches Pi 
and P2 and the remaining part. 



Fig. 4.6). We compute the weight of the patch and of the remaining part of the Brillouin zone, 
assuming that r*(k, p) is of the same order of magnitude for every value of p. Comparing the 
values 



pePi 



V 



5(|^p|-A) = 4log(4)' 



(4.5) 
(4.6) 



- y 

peB.z.-Pi-P2 

we conclude that the patch contribution dominates the remaining part if A ^ p'*. Under the 
hypothesis that the functions r^(k, k') are slowly varying, it is then consistent to replace Eq. (4.4) 

by 

'-r"(k,k')= Yl n(k)rt(k'), (4.7) 

i=l,2 

r*(k) for p e P, and 



d_ 



where we have set r^(k, p) 



1 



Z = S-(A,0) = — W(^) 



(4.8) 



as the new RG parameter. B^p{K, 0) is the value of the p-p bubble, where the loop momentum is 
restricted to the saddle point patch. 

Correspondingly, setting r^(k') w for k' e Pj in ( p7| ) leads to 



F^(k)= rnk)n 



(4.9) 



i=l,2 



Eqs. (4.7) and ( [4.9D are only correct under the hypothesis that r^(k, k') and F^(k) are slowly 
varying functions. This hypothesis is certainly justified in the beginning of the RG flow, when 
F « —g. It will be shown now that, starting from this hypothesis, we obtain an RG flow where the 
abovementioned functions remain well behaved. The hypothesis is thus consistent with the final 
result. 

We are now left with the problem of renormalizing the vertices F^^ , which describe scattering of 
particles near the saddle points. Eq. (4.4) cannot be used here because it is not valid for k±k' w Q, 



i.e. if k e Pi and k' e P2- The RG flow of the vertices F*- will be discussed in detail in Section 4.3 
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Once the evolution of the F*- as a function of / is known, one can integrate Eqs. (p 



Introducing r^(k) = r*(k) ± r^(k) we get 



r^(k) = r^(k)| 



exp 



and 



r-(k,k') = r-(k,k')b=,, + i /'d/ [n(k)r^(k') + r-_(k)r--(k')] • 



and @. 
(4.10) 

(4.11) 



We will see that the RG equations of as a function of I yield diverging solutions at a finite 
value I = lc- Near this critical value they behave asymptotically like 



F? 



+ oih-ir, 



(4.12) 



where the constant F^ can be determined from the RG equation and a > 



-I. Eqs. (4.10) and 



(4.11) then give 



F^(k)«A±.F-±(k)L^,^ . {i,-iy^±+o{h-i) 



(4.13) 



and 



F-(k,k') « F-(k,k')U,„ + E r:(k)F:(k')U,„ (/c - 1)'-'^^^ + 0{k - 



\Min{0,2-2r*+a} 



(4-14) 

where Fj. F^iiF^g ^-^d A±, B± are positive constants. The vertex function F*(k, k') is diverging 
if f* or ft > 1/2. 

The functions F^(k) have s-wave symmetry, i.e. they respect all the point symmetries of the 
square lattice. On the other hand we see that F* (k) is of the d^2_y2-wa,ve type [Tt{kx, ky) = 
^-{kx, —ky) = —Tt{ky,kx)]. Thus the diverging part of the vertex function has s or dx2^y2-'wave 
symmetry. 

In the preceding calculation we have distinguished strictly between points far from the saddle 
points and those close to the saddle points. The scale which distinguishes between "far" and "close" 
is the patch size p, which was introduced by hand. 

The behavior of the overall vertex function F^(k, k') near the criti cal p oint depends on the 
biggest positive value of the constants f j.. It will be shown in Section 4.3.2| that Max{f j.} = 1. 
In this situation the function F^(k, k') diverges everywhere with the same power {I — Ic)^^ and 
F*(k, k') remains a smooth function upon renormalization even if k or k' (or both) approach the 
saddle points. This justifies a posteriori the estimation of the p-p, p-h 1 and p-h 2 terms of Eq. 
( [3.33| ) by the bubbles ( 3.101 ) as well as the approximations of Eqs. (17) and (O), namely that 



F*(k, k') approaches continuously F*(k) as k' Fi and F*(k) approaches F*- as k — Fj. If instead 
we would find < MaxjFj.} < 1, the scattering of particles near the van Hove points would diverge 
more rapidly than that of particles at the remaining Fermi surface. This would mean that different 
regions in the Brillouin zone behave differently. The region around the saddle points would become 
strongly interacting while the remaining Fermi surface would remain weakly interacting. 



4.2.4 RG flow of the susceptibilities 

In order to calculate the susceptibilities we consider Eqs. ( |3.68 ) - (3.71). We only consider the lead- 
ing logarithmic terms, i.e. the pairing susceptibilities (Eq. ( |3.68 ) and (3.69)) at a total frequency- 
momentum fc = and the susceptibilities for density waves and flux phases (Eq. ( |3.70| ) and ( |3.7lD ) 
at a transfered frequency- momentum q = (0, Q). 
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In the low energy regime A ^ 0, we assume that the frequency dependence of Z is not important 
and replace Z{p) by Z{p). We can then exphcitly perform the frequency integral and obtain, because 
of the exact nesting, identical equations for the charge, spin and pairing susceptibilities 

2Ax" = -^5]<5(|ep|-A)Z^(p)Zl(p), 
p 

2AZ^/5(p) = --1^5(|Cp,|-A)Zl/5(p')r*(p',p), (4.15) 
p' 



where o = c,s, BCS. The index A in the symbols x, Z and T have been omitted to simplify the 
notation. 

We treat Eq. ( 4.15| ) in the same way as Eq. (4.4), i.e. we take the leading contribution from 



the patches Pi,P2 around the two saddle points and assume Z^^^(p) « Z^^^ ^ for p e P;. It 
is sufficient to consider complex conjugate form factors /a(p) = /s(p) = /(p) and thus Z\^ = 
Z'^ . = . The " more general" angle-dependent susceptibilities ( |3.50D are in fact not more general 
in this caseQ We get 



d 



2 

\ ^ \ ryO\2 



d/^ - 2.1^..^ 

J^Zt = E^;r^„ (4.16) 



which can be written as 



= ZlTl, (4.17) 

where Z^ := Z{ ± Z^ and =: x\ + X-- 

If Fj. is diverging asymptotically like rj./(Zc — I) and if Fj. 1/2, the corresponding suscep- 
tibility diverges with a critical exponent 1 — 2r : 

- {Ic-iy-'^'^ - (A- A,)i-2n. (4.18) 

On the other hand Fj. < 1/2 leads to a finite value of the susceptibility. 

We thus naturally identify six possible instabilities. In each of the charge, spin or pairing 
sectors the form factor can be either even or odd under the exchange of the two van Hove points (i.e. 
Zi = Z2 or Zi = —Z2Y In the pairing sector, x+*^'^ clearly corresponds to s-wave superconductivity 
(sSC) and X-*^"^ to (ia;2_j,2-wave superconductivity (dSC). The charge- and spin-density waves 
(CDW and SDW) are related to the susceptibilities x+ and x+, respectively. The two remaining 
susceptibilities x^l^* correspond to a form factor with d2.2_j,2-wave symmetry in the charge and spin 
sectors. They describe the tendency towards the formation of charge or spin flux-phases (CF and 
SF, see Section^. 

These six instabilities of a system with two van Hove singularities have been discussed long 



ago by H. J. Schulz |3ch8£]. Here it was shown that they appear naturally in the leading order 
renormalization group. 

^The reason is as follows. The angle-dependent susceptibility x(d, d') is replaced by a matrix Xij, where i,j G {1, 2} 
run over the saddle points. The eigenvectors of this two by two matrix (i.e. the natural form factors) are completely 
determined by the symmetry with respect to the exchange of the two saddle points. One eigenvector is odd, and the 
other is even with respect to the permutation of the two saddle points. 



4.3. FLOW OF THE VERTICES BETWEEN SADDLE POINTS. 
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Figure 4.7: Classification of tlie vertex function into Fi, . . . ,F4. 

4.3 Flow of the vertices between saddle points. 

We now return to tlie renornialization of vertices for scattering processes botli witliin and between 
saddle point patches in the nearest-neighbor tight-binding model (t' = 0). The most simple approach 
is to treat the saddle points in close analogy to the Fermi points of a one-dimensional system 
where there are just four types of scattering processes, one restricted to the region of a single 
Fermi point, the other three involving both right and left movers (forward, backward and Umklapp 
scattering) . This one-dimensional scenario with four effective coupling constants gi, . . . , (and thus 
four vertices given by F^ — —gi) implicitly assumes that the more detailed wave vector dependence 
in this region is irrelevant. While this appears to be true for one-dimensional Fermi systems, where 
going away from a Fermi point means leaving the Fermi surface, we will argue that in the present 
case the functional dependence of the vertex is relevant in the neighborhood of the saddle points. 



4.3.1 "One-dimensional" solution 



In early contributions to this subject |Sch87, Dzy87], it was assumed that the vertex function 
r(ki, . . . , k4) takes only four different values according to how the momenta ki, k2, ka and k4 = 
ki + k2 — ka are distributed over the two patches, namely 



r(ki,k2,k3) 



Fi ; ki, kg e Pi and ka, k4 e P2 
Fa ; ki, k4 e Pi and ka, kg e P2 
F3 ; ki, k2 e Pi and kg, k4 e P2 
F4 ;ki,...,k4 e Pi 



(4.19) 



or symmetry-related configurations of the external momenta. The four corresponding scattering 
processes are illustrated Fig. 4.7. The parameters F^. 
instabilities (o = s, c, BCS), are readily expressed in terms of the vertices Fi 



F^j^ ± F*2, which control the various 



SDW/SF 
CDW/CF 
sSC/dSC 



= -r2 T Tg 

= 2Fi - F2 ± Fg 

= F4±Fg. 



(4.20) 



Fi has to be renormalized by the diagrams p-h 2, because the direct momentum transfer kg — k2 
is close to the nesting vector Q, but the other contributions coming from p-h 1 and p-p are negligible. 
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Similarly T2 is renormalized only by p-h 1 and r4 by p-p. The remaining vertex on the other 
hand gets leading contributions from all three channels, p-p, p-h 1 and p-h 2. 

The RG equation is obtained by locating the external momenta ki, . . . , k4 in Eq. (3.33) exactly 
at the saddle points Pi and P2 and restricting the sum over p to the two patches Pi,P2- The result 
is 



1^. 



2ri(ri 



2r3(ri-2r2-t-r4 



- 4- 



(4.21) 



For most initial conditions the numerical solution of these equations diverge asymptotically like 
in Eq. ( 4.12 ) with coefficients Ti satisfying 



fa 

f4 



2fi(fi-f2) 



f2 
^ 3 



2r3(ri-2r2 + r4 



(4.22) 



Eq. ( 4.22 ) has many solutions, but the ones which are relevant for the divergences of Eq. ( 4.21 ) 
are Fi = 0, r2 = — r4 = — 1/6 and Fs = ±-\/5/6, depending on whether the initial value of F3 is 
positive or negative (no te th at F3 c anno t change its sign). The special feature of these two solutions 
is that in view of Eqs. ( 4.18 ) and ( 4.20 ) three out of the six dominant susceptibilities are diverging 
with the same critical exponent. Namely 



XSDW ' 
XCDW 



XdSC ■ 
' XsSC 



XCF 
• XSF 



(A - Ae)-^ 
(A-A,)-^ 



if 
if 



Ta < 0, 
Ta >0, 



(4.23) 



where 7 = {-sfb — 2)/3 « 0.08. Note, that Fa < corresponds to a repulsive coupling ga > 0. The 
divergence of the susceptibilities is thus weak compared to the m ean fi eld behavior x^ [T ~ Tc)^^. 

However, for some initial conditions the solutions of Eq. (4.21) are not diverging but flow 
towards the trivial fixed point Fi = F2 = Fa = F4 = 0. This means that the RG flow of the vertices 
between the two saddle point patches does not develop an instability. In this case, the flow can be 
followed down to small energies such that F log A ~ 1 and then, the restriction to the saddle-point 
patches is no longer valid. It means that the low energy behavior can be controlled in this case by 
subleading terms such as the part of the Fermi surface which is far away from the saddle points, 
self-energy effects or others. 



4.3.2 Towards functional renormalization 

The hypothesis of constant vertices Fi, . . . , F4 neglects the fact that in reality all these parameters 
are functions of incoming and outgoing momenta ki, . . . , k4, three of which can move freely on the 
Fermi surface within the saddle point patches. Unfortunately, a true functional renormalization 
is presently beyond the reach of an analytical approach. I therefore search for an appropriate 
simplified characterization of the momentum dependence of the vertex. 

The p-p bubble, which rcnormalizes the vertex varies most rapidly near k = 0, as it can be seen 
in Fig. 4.5. It is thus natural to admit that the vertex F(ki, k2, ka) depends strongly on the total 
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Figure 4.8: Schematic view of the space of three momenta ki,k2 and ka moving along the Fermi 
surface, where the function Fa is defined. The three planes correspond to the special relations 
ki + k2 = 0, ka — ki = Q and ka — k2 — Q. The point of intersection of the three planes 
corresponds to the configuration, where the momenta are exactly located at the saddle points. In 
the approximation proposed here, the vertex takes a different constant value (Fg, F| or F;f'-^'^) on 
each of the three planes and yet another constant value Fa in the remaining phase space. 



momentum ki + k2 in the vicinity of ki + k2 =0 and similarly on the two transfered momenta 
ka — ki and ka — k2, if one of them is close to Q. A simple way of modeling this behavior is the 
following. The vertex is assumed approximately constant except in the vicinity of one of the three 
special configurations (ki + k2 = 0, ka — ki = Q or ka — k2 = Q). 

We therefore mimic the true momentum dependence by introducing, in addition to the constants 
Fi, . . . , F4 representing general values of the momenta (in the patches), other vertices corresponding 
to specific combinations of momenta. For example, we allow Fa(ki,k2,k3) to take four different 
values: 



F3(ki,k2,ka) 



Ta 

pBCS 

^ a 



if 
if 
if 
if 



|ka 
ki 
ka- 
ka 



k2 

k2 = 
ki = 
k2 = 



-Ql, 


Q 
Q 



Ika-ki-QI, |ki+k2| >0(a/A) 



(4.24) 



The approximation is illustrated schematically in Fig. 4.8 



Similarly, Fi(ki, k2, ka) takes the values Fj or Fi, F2(ki, k2, ka) — F2 or F2 and F4(ki, k2, ka) 



^BCS 



or F4. We thus separate the vertices Fg 



BCS 



rf-^^ with zero total momentum, Tf, 



Ff with 



a direct momentum transfer equal to Q and Ff,F| with an exchanged momentum transfer of Q 
from the general ones (Fi, . . . , F4), where none of these special relations among the in- and outgoing 
momenta applies^ In the transition domain (where for example < |ki + k2| < \/A) the function 

is unknown, but, as I will argue shortly, its knowledge is not essential. 

The special vertices are related to the vertices Fj. introduced in Section 4.2.3| by 



SDW/SF : 


Fi 


_ ps _(_ ps 
— ^ 11 =■= J- 12 


T^X -r- T^X 

— 2 + a 


CDW/CF : 




_ pc _(_ pc 
— ^ 11 ^ 12 


= 2Ff - F| ± (2Ff - Ff ) 


sSC/dSC : 


pBCS 


_ j^BCS 1 j^BCS 
— ^ 11 =■= ^ 12 


_ j^BCS 1 Y-BCS 



(4.25) 



^Other special configurations of the momenta will be considered later in Section 4.5 
different vertices and the notations is given in Appendix ^ 



An overview of all the 
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k' 

PP: p 



-k' k' 
-p PHI: p 
-k k 



-k 

f p-k-k' 
-k' 



Figure 4.9: Two diagrams involved in Eq. (4.26). k and k' are typical vectors belonging to the 
patch Pi and P2, respectively; p is the integration variable 



Since r^'-''^,r^ and get the strongest contributions from the p-p, p-h 1 and p-h 2 channels, 
respectively, we say that r^'^'^ is resonant in the p-p channel, in the p-h 1 channel and F'^ in 
the p-h 2 channel. Furthermore the term with the largest contribution, say Rp^ for a BCS vertex, 
again only includes BCS vertices and does not mix with non-BCS processes 

For example, let k g Pi and k' G P2 such that r(k, — k, — k') 
equation for this process: 

-^ri^cs ^j^PP^j^phi^j^ph2^ .4 26) 

dA ^ ' 



Tf^^ and consider the RG 



The dominant p-p diagram, shown in Fig. 4.9, involves the vertices F'^'-^'^ (k, p) and r-^'-^'^(p, k'). 
Within our approximation they are replaced by the constants Tf-^^ or F^'-^'^, respectively. The p-p 
contribution to Eq. (4.26) becomes 



BCS-rBCS 



RPP = 2BPP{A,0) Ff'-^F 



(4.27) 



where B^{A, k) is the p-p bubble restricted to a saddle point patch of size p, given by Eq. ( 3.119 ) 
for t' = 0. 

By contrast, the non-resonant diagram p-h 1 (also shown in Fig. 4.9) involves vertices like 
F(k, p — k — k', — k') which for almost every value of p do not satisfy one of the special relations 
ki + k2 = 0(\/A), ks - k2 = Q + 0{^/X) or ka - ki = Q + 0(VA). This diagram therefore includes 
only the general vertices F2 and F3 and no special vertices p^C'-S'^p'* or F^. Thus 

RP^^ ^ 2Bf{A, k + k') F2F3, 



(4.28) 



where (A,k) = -B^^(A,k - Q). It follows that the RG flow depends on the ratio (A,k 
k')/i3p^(A, 0). In order to obtain a closed set of equations we replace this ratio by a constant. 
Let us first define 

.(A,k + kO:^^^^%|±^, (4.29) 



BPP{A,Q) 



which varies in principle between and 1. Its value is 1 if k -|- k' = Q, i.e. if the process to be 
renormalizcd is at the same time a F^'^'^ and a F^. Such processes exist of course, but they will 
not influence the RG equations in a relevant way. According to Eq. (4.5) a region in the Brillouin 
zone of width ^ -s/A or smaller can be safely ignored within logarithmic precision. We thus assume 
that k -I- k' — Q ~ -s/A or bigger. The biggest values are obtained if k -|- k' is parallel to the Fermi 
surface k + k' — Q (k, k). For k > \/A we get from Eq. (3.119) 



a(A,k + k') < 



log 



log 



A 



(4.30) 



In the following we replace q;(A, k -I- k') by a constant a < 1/2. 
are treated in the same way, i.e. they include general vertices Fi, 
reduced with respect to the resonant ones by a factor a < 1/2. 



All the nearly resonant diagrams 
. . . F4 only and their amplitude is 
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Our approximation scheme leads to a set of RG equations for the special vertices Vf, , T^'~'^ , 

r^, rf and rf 



d_ 

di 
d_ 



J- 1 



2) 

2 /™^2 



2r^(r^ 



(r^)^ - (r^)^ 



2pBCSpi3CS 2ar3(ri - 2r2) 



-2rfr^ + 2ar3(ri-r2 + r4) 



(4.31) 



^BCS 



dl 



2(2r?r 
fr 



2ar3(r4-r2) 



The general vertices Fi, . . . , r4 ar e resonant in none of the three channels. The RG flow of 
Fi, . . . , r4 is thus given by Eqs. ( 4.21 ), with the right hand side multiplied by a. 

Eqs. ( |4.3l| ) can be rewritten in terms of the vertices which are associated with the dominant 
instabilities (see Eq. (4.25)) 



di^^ 
di^^ 

rBCS 



dl 



(ri)'T2ar3(ri-r2 + r4) 
(r±)'T2ar3(ri + r2-r4) 
(pBcs)2^2ar3(2r2-ri). 



(4.32) 



For a = 0, this is a set of six independent equations, one for each instability. In fact, if the non- 
resonant diagrams are completely neglected, the RG becomes equivalent to the summation of ladder 
diagrams. Even for < a < 1 the special vertices associated with the different instabilities still do 
not influence each other, but each RG equation has a source term coming from the general vertices 
Fi, . .^JU- For a = 1 and initial conditions Ff = Fi, Ff = F2, Ff = Ff = Ff^^ = F3, Ff = F4, 



Eq. (4.21) is perfectly recovered (since these conditions are then conserved by the RG flow). 



4.3.3 Results 



One can search for asymptotic solutions of the form F(/) = f • (/^ — l)^^ of Eq. ( [1.32| ) by solving the 
resulting algebraic equations for the F. We first consider the possibility of diverging general vertices 
Fi, . . . ,F4. In this case it follows from our analysis of Eq. (4.21) that the asymptotic behavior of 
the general vertices is given by fi = 0, f2 = -f4 ^ -l/(6a) and f3 = ±^/5/(6a), depending on 
the sign of F3. By inserting this behavior into Eq. (4.32) it is easily seen that a real solution for 
Fj_,Fj_ and Y^^° requires a > v/80/81 w 0.994. But as we argued above, the appropriate values 
of a are < 1/2. 

It follows that for acceptable values of a a special vertex can only diverge if Fi = . . . = r4 = 0. 
This means that this special vertex diverges at a higher energy scale than the general vertices. 

To illustrate the behavior of the RG equations for different values of the parameter a, I have 
plotted some numerical solutions of Eq. (4.31) in Fig. 4.1C. The initial conditions have been chosen 
such that all the vertices are initially equal to —1, i.e. the bare value of the repulsive Hubbard 
model at U = t = 1. The RG flow is quite similar for a = 1/2 and for a = (Fig. 4.10| b) and c)). 
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a) 



b) 



0) 




pBCS _ -nBCS 
^3 — ^ 4 

Fi = • • • = r. 



pa: pd pa: pBC5 



and 



pBCS 



as a 



Figure 4.10: Three plots of the ten vertices Fi, Fs, Ff , i 2, ^ 3, j 
function of the RG variable I, for different values of the parameter a. All the vertices are ini- 
tially chosen equal to —1. The six non diverging vertices of Figure b) are, from top to bottom, 
F4,Fi,FfC"5,r3, andFs. 



rBCS 
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Both RG flows are governed by the divergence of the and vertices, while other vertices stay 
smalj^. I conclude that the RG flow depends little on the parameter a in the range < a < 1/2. 

I n con trast the case a = 1 (Fig. 4.1C| a)), which is equivalent to the "one-dimensional" approach 
E g. ([4.2l| ), shows a completely different flow. If one w ould choose a close to 1 but still smaller than 
■y/80/81, the flow would first resemble the one of Fig. 4.10| a) and the dominance of the F^ and F'^ 
vertices would start only close to the critical scale Ic. But for acceptable values < a < 1/2 the 
flow differs from Fig. 4.1Cl| a) already long before the critical scale. 

The most striking difference to the "one-dimensional" solution is that some of the six special 
vertices can diverge, while the others remain flnite. This occurs here because the mixing of the flow 
for these ve rtice s has been neglected on the basis of a phase space argument (see the discussion 
before Eq. ( 4.28 )). This argument is certainly valid as long as the vertex function is slowly varying, 
but it may be quest ioned close to the instability, where the vertex function gets peaked. It is argued 
however in Section 4.3.5 that the non leading vertices can nevertheless stay finite in the case of an 
instability. 

One important feature however is shared by Figs. 4.10 a) and b), in contrast to Fig. 4.10 c). 
Namely the RG flow for a > leads to a positive value of the quantity F^*-^"^ — Tf'-'^ . This would 
lead to d-wave superco nductivity if there was no competing SDW. The main conclusion of Schulz's 
original paper | Sch87 that the Hubbard model with an electron density slightly lower than one 
becomes a d-wave superconductor remains valid in our more elaborate sche mq^. 

The divergence of the leading vertices is characterized by Fj. = 1. Eq. (4.1S) then implies the 
asymptotic behavior x ~ (A — Ac)~^, corresponding to a mean field exponent. Fig. 4.11 shows the 
RG flow of the susceptibilities xsdw, Xdsc and xcf in the case of Hubbard- like initial conditions 
(all the vertices equal to —1). It has been obtained by solving Eq. (4.17). The SDW susceptibility 
dominates the two others in both cases a = 1 and a = 1/2. But whereas all three susceptibilities 
diverge with the same asymptotic behavior ~ {Ic — lY^^^^i/^ in the "one-dimensional" solution 
a = 1, the non-dominant ones saturate in the more realistic case a = 1/2. The difference between 
a = 1 and a = 1/2 is less pronounced for the susceptibilities than for the vertices, but it is clearly 
visible. For a better comparison, see Fig. f4.12| where the ratio between the non-leading and leading 
susceptibilities is shown as a function of l/lc- 

Finally, in view of Eq. (4.14) the coupling function F*(k, k') diverges as (A — Ac)^^ everywhere 
on the Fermi surface and remains a smooth function upon renormalization, as anticipated in Section 



In conclusion, the RG scheme presented in Section 4.3.2 interpolates between the "one dimen- 
sional" solution and the generalized ladder approximation (or generalized RPA), and shares features 
of both of them. Like in the "one dimensional" solution, our RG flow generates enhanced supercon- 
ducting fluctuations in addition to the leading SDW and strongly suggests a SDW-dSC transition 
upon doping in the repulsive Hubbard model. But like the RPA, our RG flow leads to only one 
diverging susceptibility and to the mcan-flcld critical exponent x ^ {T — T^)^^. 



4.3.4 Phase diagram 

For a curved Fermi surface away from the van Hove singularity, the dominant instability is su- 
perconductivity, as discussed in Chapter ||. For the square Fermi surface the flow equations ( 4.32 ) 



show that there are several possible instabilities, s- and d-wave superconductivity, charge- and spin- 



^The case a = features the exact degeneracies Ti = r2 = Ta = rf<^^ = Tf^^, = Tf and Tf = r|. 
This is an artifact of both, the special initial conditions and of a = 0. 

■^The argument for the occurrence of superconductivity was as follows. In a doped system with an electron density 
slightly lower than one, the RG flow will first be similar to the half-filled case and Fj"^^ — F^'-^^ develops a positive 
value. Finally for low enough energies, the competing density-wave and flux-phase instabilities are suppressed because 
of the absence of nesting and d-wave superconductivity can develop unhamperedly. 
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Figure 4.11: RG flow of the susceptibilities xsdw, XdSC and xcf- Two of them, XdSC and xcf, 
are exactly degenerate for the initial conditions at hand. This is due to a special symmetry as 
explained in Section 4.4. The location of Ic is indicated by a vertical line. 
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Figure 4.12: Plot of the ratio between the non-leading and leading susceptibilities as a function of 
l/lc. It shows a qualitatively different behavior for a — 1 and a — 1/2. 
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Charge Flux 
CF 




4V/U 



Figure 4.13: The phase diagram for [/ > and a = 1/2. In the dotted region aU the vertices 
flow to zero. The size of the dotted region depends on the parameter a, the other features are 
a-independent 



density waves, charge and spin flux-phases. The values of the initial vertices will determine which 
of these instabilities, if any, occurs flrst, i.e. at the largest energy scale. 

The initial conditions for the flow equations are given by the vertices Fi, . . . at a cutoff 
^0 1 when the flow enters the asymptotic regime. They are calculated by naive perturbation 
theory, i.e. to lowest order Fa — —5, where g is the microscopic coupling function. This choice 
is a good approximation if the interaction is small enough such that the vertices vary little before 
entering the asymptotic regime. 



To be specific I consider the interaction given by Eq. (2.13) consisting of on-site and nearest- 
neighbor terms. It yields the following starting values of the vertices 



Fi = F? = -[/ + 4y + J + 

= -U + AV - 3J - 4W, 



F2 = Ff ^ -C/-4y - J + 4PF, (4.33) 



BCS 



= Ff = -[/ - 4y + 3 J - 4W. 



The Coulomb-assisted hopping term K does not appear in Eq. ( 4.33| ), because it vanishes exactly 
at the Fermi surface of the half-filled nearest-neighbor tight-binding band. 

I have solved Eq. (4.32) numerically for various initial conditions and obtained the phase diagram 
as a function of U, V and J. The result is shown in Figs. 4.13 and 4.14 for, respectively, positive and 
negative values of U. W has been put to zero. The predicted phase for the repulsive Hubbard model 
{U > V = J = 0) is a spin-density wave (SDW) as expected and quite strong nearest neighbor 
terms are needed to establish a fiux phase or (only for attractive V) a d-w&ve superconductor. An 
unexpected feature of Fig. 4.13 is that the SDW can be destabilized by positive values of J. 



In contrast to the "old" RG equations (4.21) our more elaborate scheme Eq. (4.31) produces 
only one diverging susceptibility while the others remain finite. Only at the phase boundaries, 
where the two neighboring phases are degenerate, both susceptibilities diverge. 



In a certain parameter range (the dotted region in Figs 4.13 and 4.14) all the vertices flow to 
zero. In this case the behavior is not necessarily dominated by the saddle points and subleading 
contributions to the RG flow such as self energy corrections can play an important role. 

It is interesting to note that most features of the phase diagram are independent of the parameter 
a. In fact, only the size of the dotted regime depends on it, but neither of the transition lines 
between two neighboring phases. The reason is that these transition lines are entirely defined by 



some symmetries, a point which will be further explained in Section 4.4 
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J/IUI 




Figure 4.14: The phase diagram for U <0 and a — 1/2. 



Since the phase diagram is independent of a, we can choose a — 0. There it is clear, that the 
leading instability is given by the biggest positive value of the initial vertices F , namely 



sSC : 


irBCS 

2^ + 


dSC : 


IpBCS 
2^ - 


SDW : 


2^ + 


SF : 


2^ - 


CDW : 


1 "pc 
2-L + 


CF : 


1 pc 
2-^ + 



= -U -AW 

= -AV + 3J 

= U + 2J 

= AV - J -AW 

^-U + W 

= 4F + 3 J + AW. 



(4.34) 



A numerical integration of Eq. (4.32) is only required to check whether the vertices diverge or flow 
to zero. 

Let me conclude this section with a remark about superconductivity in a slightly doped system, 
where the density-wave and flux-phase instabilities are suppressed. As we have seen in Chapter ^ 
a nearest-neighbor repulsion 1/ > is efficient to suppress d-wave superconductivity. This is also 
evident in Eq. ( [4. 34] ). However, it is easy to see from Eqs. ( [4.32| ) and ( [4.2l| ), that t'^^^ > 0, if 
initially F3 < and 2F2 — Fi < 0. This means that the c?-wave superconducting fluctuations are 
enhanced by the RG flow even in the presence of nearest-neighbor repulsions AV < U (assuming 
J = W — and U > 0). I have checked that for V < 0.18U, the enhancement is sufficient to 
turn the initially negative value of T^'~'^' into a positive one before the SDW instability takes place. 
In this parameter regime, d-wave superconductivity is the most favorable instability for the doped 
system. 



4.3.5 A consistency test 



The approximation scheme presented in section 4.3.2 leads to a complete decoupling of the p-p and 
p-h channels for a — 0. In this case, the solution of the RG equation is given by the summation of 
the p-p or p-h ladders Figs. 



2.2, 3.2 and 3.3 



We have seen that, although Eq. (4.31) for a > goes beyond the RPA, its solutions are 
asymptotically the same as for a = 0. The reason is that in our approximate treatment of the non- 
resonant diagrams only the general vertices Fi , . . . , F4 intervene and not the (diverging) special 
vertices Ff , Ff etc. This approximation is certainly justified in the beginning of the RG flow, i.e. 
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Figure 4.15: Plot of the vertex r(A, k, k) normalized by its maximal value at k = 0, for the choice 
of parameters Aq = O.lp^, 5 = — 1 and different values of A. The vertex depends relatively weakly 
on K for A far above Ac, whereas the peak at k = gets more and more pronounced and narrow 
as Ac is approached. Fig. 1.15 is in nice accordance with our approximation, where the vertex 
function is replaced by a constant for A/2 < k < p. The thick points indicate k = ^ K/2. 



as long as the vertex does not vary much as a function of the momenta. However, in our solution 
we obtain diverging special vertices and finite general vertices, i.e. the vertex functions get a strong 
peak close to ki + k2 = 0, ka — k2 = Q or ks — ki = Q, depending on the nature of the instability. 
Using the schematic picture of Fig. the vertex becomes peaked at one (or several) of the three 
planes. It is not a priori clear whether our approximation is justified, once the vertex is strongly 
peaked. 



For example the (non-resonant) contribution of the diagram p-h 1 of Fig. l.£ to the renormal- 
ization of V^'-^^ involves an integration over the vertex function r(k, p — k — k', — k', p), where p 
is the integration variable moving along the one-dimensional energy shell |^p| = A. It is an integral 
over a one dimensional curve in the space of momenta (ki,k2,k3). For special values of p, the 
curve crosses one of the planes specified by ki -I- k2 = 0, ks — k2 = Q or ks — ki = Q. For example 
for p = k' the vertex is equal to F^'^'^. The question is whether one can neglect the contribution 
close to this point even when T^'-'^ diverges. 

In order to answer this question we have to know the value of the vertex for small but finite 
total momentum. Within the ladder approximation one obtains 

where q is the (small) total momentum^ and g^'^^ = ff(Pi,q— Pi , q — Pi)ibff(Pi,q — Pi,q — P2) 



Eq. ( 1.35 ) can be obtained either by solving the RG equation ( 1.3S ) for a — , wh ere the RG 
variable / is replaced by _Bp''(A, q), or by explicitly summing the ladder series ( 2.46 ) where the 
integral over the angle 9 is replaced by a sum over the two saddle points. 

A similar expression is obtained for F" and F'^ as functions of the deviation of the momentum 
transfer from Q = (tt, tt). The couplings _|_ depend weakly on q and can thus be approximated 
by their value at q = 0. I will therefore study the function 

^(^'^) = TT7^fei)' (4-^^) 

For negative values of 5, the vertex diverges at a critical scale Ag. For A > Ac, F(A,q) has a 



maximum at q = 0, which diverges for A — > K^. Fig. 1.15 illustrates the behavior of the vertex for 



""q = p — k' in the example of section 4.3.2 
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Figure 4.16: Numerical evaluation of Iging and I^eg for Lq = O.lp^ and g = —1 as a function of 
log A^Tp- (Taking smaller values of Aq or go reduces the weight of Ising with respect to Ireg-) 



q = (k, k) parallel to the Fermi surface where, from Eq. (3.119), the p-p bubble is given by 



dA- 



1 ^ 4p2 
■log^ 



A (27r)2A A + 2pK 



(4.37) 



We estimate the contribution of the peak in r(A, q) to the non-resonant diagram by integrating 
this function over a curve in q-space. The biggest contribution is obtained from q = {k, k) (parallel 
to the Fermi surface). Within this worst case scenario, we compare the contribution from the 
vicinity of the singular peak 

/s^ng(A) - / dKr(A,K,K) (4.38) 

with the regular remaining contribution 



Ireg (A) ^ I dKT{A,K,K) 
'VA 



(4.39) 



The non uniform density of states is not tak en int o account for the moment. The integrals can be 

„ is suppressed with respect to Ireg j because 



calculated numerically and are shown in Fig. 4.16. / 
the integration range is much smaller. Finally for A 



A„ 



is diverging like log (A — Aj,) (note 



the positive slope in Fig. 4.16) whereas Ireg saturates. But the divergence of Iging manifests itself 



only exponentially close to Ac and is in any case negligible compared to the resonant diagrams 
~ (A — Ac) ^. I conclude that the approximation presented in section 4.3.2, namely evaluating the 
non-resonant diagrams with the constant general values Fi, . . . ,r4, is consistent. 

In view of the preceding analysis, the RG equation of a non- leading vertex is affected by the 



singularity through a logarithmically diverging term, i.e. log(A — Ac). Since J^J" dAlog(A 



Ac) < oo, this confirms the result stated above in Section 4.3.2 that non leading vertices remain 
finite at Ac. The same applies to non-leading susceptibilities!^] 

®The result that Ising diverges as log(A — Ac) can also be obtained analytically by noticing that the vertex function 
behaves close to the peak like 

r(A,K, k) : 



A - Ac + 



Mil 



where a and b are positive constants. This behavior is an artifact of our zero-temperature calculation. If the same 
calculation would be done for A = but a finite temperature T, the vertex would be analytic in q and behave like 



rT(A = 0,q): 



( T I 
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Note that this behavior can only be obtained because the momenta are aUowed to move con- 
tinuously on the Fermi surface. If instead we would discretize the Fermi surface and replace the 
continuous vertex function by a finite set of constants, a divergence of one vertex Fc at a scale Ac 
would imply the divergence, at the same scale Ac, of all these vertices that have Fc appearing in 
the right hand side of the RG equation. 

Recently, the discretized RG has been studied in detail for a simpler model without van Hove 
singularities and without Umklapp scattering [DVdAD02|. It was found that there is a factor 
between the biggest non-dominant vertices and the dominant ones, where N is the number of 
patches. This is consistent with our result in the continuous case {N — > oo), that non-dominant 
vertices stay finite at Ac while the dominant vertices diverge. Similar results were also found in the 



large- hmit of half-filled A^-leg ladders [ ]Led01 |. 

Taking into account the non uniform density of states could increase the contribution of the 
peak, if K = corresponds to the maximal density of states. But this implies k' — (0, tt) in the 
example of Fig. 4.£ and is typically not the case. In order to test the really worst case, I have 



evaluated the following integrals 



I sing (A) 



F(A, K, k) 



(4.40) 



and 



IregW 



VA \/A + K- 



: F(A, K, k) 



(4.41) 



The denominator in the integrand mimics the situation where the path of integration |fp| = A gets 
close to the va n Ho ve singularity at the same time as q = p — P.; approaches zero. The result is 
shown in Fig. 4.17 , The singular part is no longer small compared to the regular part, but the 
diver gence of Ising is of course still not stronger than ^ log (A — Ac). As I mentioned before, Fig. 



4.16 corresponds to a more typical situation then Fig. 4.17 



where c is another positive constant. But in a finite temperature calculation, the right-hand side of the RG equation 
is not given by one-dimensional integrals over energy shells |^p| = A, but over a region \t;p\ < T with a finite width. 
As a consequence, the contribution of the peak in ry(q) is given by a two-dimensional integral 

Is^ng = f d^qFrCq) '^"^'= > ~ log{T - Tc). 



It shows the same logarithmic behavior close to the instability as in the zero temperature case with a cutoff. 
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4.4 Special symmetries 



I will now discuss the phase diagram at half filling (Fig. 4.13 and 4.14) in terms of special symme- 



tries on the lines separating two different phases. For that purpose I consider the p-h symmetric 



Hamiltonian of Eq. (3.107) 



The symmetry group of the noninteracting Hamiltonian Hq is extremely large. From any func- 



tion dk we can build an operator = l/2X]o-k'^kC 



(Tk^CTk 



that commutes with Hq and thus 



generates a continuous group of symmetry transformations exjp(iaNd). We find that Nd commutes 
with the complete Hamiltonian H if and only if 



(dki + rfka - d]^3 - dki) 5(ki,k2,k3) = V ki, . . .k4, 

where here and in the following k4 = ki + k2 — ka. 

Let the order parameter for superconductivity with a form factor / be given by 



A 



(4.42) 



(4.43) 



The symmetry Nd relates superconducting order parameters with form factors / and / • c? by 
[Nd,A.f] = — 2A(j./, where d ■ f means the point-wise multiplication of the two functions. This 
relation can be iterated to [Nd,Ad.f] — —2Ad2.f. If d is chosen as a sign function, such that 
d^(k) = 1, then the two order parameters A^^ and A^.^^ are related by the symmetry transformation 

(4.44) 



Ad.f = exp {iTTNd/4:)iA / exp (-i7r7V<j/4). 



A proof is most easily given by solving a second-order differential equation for the operator 0(0) = 
exp {i(j)Nd)Af exp {-ic^ Nd). 

Clearly, Eq. (4.44) transforms an s-wave superconductor into a c?— wave superconductor and 
vice- versa. As a consequence the susc eptibilities for s- and d-wave superconductivity must be 
exactly equal, provided condition ( 4.42 ) holds. The symmetry Nd relates in a similar way spin- 



or charge-density waves to the corresponding flux phases. It might therefore control the transition 
hues SDW/spin flux-phase, CDW/charge flux-phase and sSC/dSC. 



Similarly the operators Sd — 1/2 X], 



cr,(T' ,k 



'^Irk'^a-a' Ca'k relate the spin-density wave to the 



charge flux-phase and the charge-density wave to the spin flux-phase. They commute with H if the 
following two conditions hold 



(dki - rfka + dk?. - dkjgi^i, . . . , k4) V ki, . . . k4 

(dki + rfka - dk3 - dkj (1 - X)g{'ki, . . . ,k4) = V ki, . . .k4. 



(4.45) 



Another symmetry introduced by Lieb | Lie89| and then further investigated by Yang and Zhang 
YZ90, Zha90| is generated by the pseudo spin operator rjs = X^k '^k C|Q_kCxk, where the function 
Sk satisfies SQ_k = Sk- It turns a s-wave superconductor into a charge-density wave and a d-wave 
superconductor into a charge flux-phase (and vice versa), rys commutes with Hq because of the 
exact nesting (^q-k + ^k — 0) and it commutes with the full Hamiltonian provided 



EpSp.9(P>P>k,k) = SkEp(2-_^)5(k,P,P,k) Vk_ 

Ski(l - -'^)5(ki, . . . ,k4) - Sk3g(k3,k2,k4,ki) + Sk45(k4, k2, kg, ki) = 



(4.46) 



V ki, . . . k4, where k := Q — k. 



Finally Zhang [Zha97] considered the operators H^ — 1/2 X^cr ct' k CaQ-k ("tt^)^^, Ccr'k con- 



necting a spin-density wave to a d-wave superconductor and a spin flux-phase to a s-wave super- 
conductor. The symmetry condition is of the same form as ( |4.46| ) but with Ski replaced by a 
function dki that satisfies dq-k = — dk. 
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It is in general difficult to satisfy the conditions ( |4.42| ) to (4.46). For example tliey do not hold 
for the U — V — J — W interaction Eq. (2.13) (the Coulomb assisted hopping K is supposed zero, 
because it breaks the p-h symmetry). The only exception is the pseudo spin symmetry 77^ which is 
exact for V = W — and Sk = 1. 

However the restriction of the model to the two saddle point patches has more chance of being 
symmetric. We take Sk = 1 everywhere whereas d(k) = 1 for k G Pi and d(k) = — 1 for k S P2- 
For this simple choice the symmetry generators Nd, Sd, rja and Hd together with the total spin- 
and charge operators form a so(6) ® so(2) Lie algebra. The commutation relations of the symmetry 
generators and the relevant order parameters are listed in reference [MV98|. 

We now assume that the coupling function .g(ki, . . . , k4) takes only four different values gi, ... 174 
as indicated in Fig. [1.7| . The symmetry conditions are then (?3 = for Nd, 51 = for Sd, 52+54 = 251 
for "qs and 52 + 54 = for Hd- These hyper-planes in our four dimensional coupling space define 
exactly the transition planes of the phase diagram (shown in Figs. 4.13 and 4.14 for 51, ... ,54 
parametrized by U, V and J) . 

We have thus shown that the transition planes of the phase diagram are fixed by exact symme- 
tries of t he ffi , . . . ,54- model. This is a strong indication that the phase diagram shown in Figs. 
4.13 and 4.14 is the correct one at sufficiently weak coupling. Such a determination of an exact 



phase diagram by simple symmetry considerations was also possible for a one-dimensional system 
i BFSOj . 

Let us n ow loo k at the sy mmetr ies of the effective interaction 5A(ki,k2,k3) = — rA(ki, k2, ka) 
(see Section ^.4.2 ). In Section 4.3.2 , we have parametrized this function in terms of ten instead of 
four constants. The symmetry conditions for the effective interaction read 



Ti = T'i = T 



BCS 



2Ti - Ta - r4 = 2rf - 

r2 + r4 



Fi 

= F? + Ff : 



^ 1 



^ 3 



= F3 = 

3 - - 

Ff - 



for 
for 
for 
for 



Nd, 
Sd, 

lis, 



(4.47) 



These symmetries are respected by our approximate RG equation (4.31) for every value of the 
parameter a. In fact it is easy to show that if one of these conditions is satisfied at the initial scale 
Iq it remains to be so at any scale This is not completely trivial since an arbitrary approximation 
scheme might violate the symmetries of the model. 



4.5 Uniform susceptibilities, Fermi surface deformations and 
?7-pairing 

4.5.1 More conventional and generalized susceptibilities 

So far we have identified six leading instabilities of the half-filled nearest-neighbor tight-binding 
band: spin- and charge-density waves (SDW and CDW), spin and charge flux- phases (SF and CF) 
and finally s-wave and d-wave superconductivity (sSC and dSC). Each of these instabilities is now 
further characterized by looking at some additional physical quantities. The framework is linear 
response theory, which was presented in Section In the preceding sections we have focused on 
the spin- and charge susceptibilities at momentum Q = (tt, tt) and pairing susceptibilities at zero 
total momentum. The quantities addressed here are obtained from them by a momentum shift of 
Q, i.e. we consider spin and charge susceptibilities at small momenta q — > and pairing with a 
total momentum of the pairs close to Q. 
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The uniform spin susceptibility xs is defined by 



Xs= limx^(go = 0,q), (4.48) 



where x'^il) is the spin susceptibiUty of Section 3.3 with a trivial form factor /a(p) = /-b(p) = 1 



The small momentum limit needs to be specified here, because the p-h bubble BP'''{q) has no well 
defined limit q ^ 0. From Eqs. ( 2.67| ) and (2.68), one obtains 



limBP\qo^O,q)^^^-u{T), (4.49) 

q^O 



whereas BP'^{qo, q = 0) = 0. It is important that there is no infrared cutoff A in Eq. (4.49), because 
otherwise, the result would be zero rather than —v{T). Therefore, in the remaining part of this 
section, we assume a finite temperature T but A = 0. 

An enhancement of the uniform spin susceptibility indicates ferromagnetic tendencies, and a 
suppression of this quantity signals the absence of low energy spin excitations (spin gap). 

The uniform charge susceptibility is the analogous quantity in the charge sector. The external 
field then couples to the total particle number N and corresponds to a change in the chemical 
potential. The charge susceptibility is thus defined as 

1 d{N) 

XN = lim Y^(<Zo = 0,q) = 7^^, (4.50) 
where the partial derivative is taken at constant temperature and system volume. This quantity 



is often called charge compressibility. The relation between Eq. (4.50) and the compressibility 
is explained in Appendix ^. A suppressed charge compressibility at low temperature signals a 
charge gap and is symptomatic of insulating systems (both band- or Mott insulators) . A diverging 
compressibility would indicate phase separation (see Appendix ^ . 

We also investigate the linear response with respect to the operators Nd and Sd, which have been 
introduced in the preceding section. These operators are obtained by introducing a d^2_y2-wave 
form factor into the particle number and total spin operators, respectively, in a similar way as one 
obtains the charge and spin flux-phases by introducing such a form factor into the corresponding 
density waves. The susceptibilities XNd ^^'^ XSd have the same definition as xn and xS: but with 
/a(p) = /b(p) = cos Pa; — COS Py (oT another smooth function with the same c?-wave symmetry). 
They have the following interpretation. 

If xnj is diverging, the system is likely to establish a non-zero value of (Nd) ~ (A^i) — {N2), 
where Ni is the occupation number of the saddle point patch Pi. If one patch, say Pi, is more 
occupied than the other, it means that the electrons are moving more easily in the x-direction than 
in the y-direction. I interpret this as a spontaneous deformation of the square Fermi surface towards 
an open Fermi surface as shown in Fig. EJ, breaking the discrete D4 lattice symmetry. Such an 



instability (called Pomeranchuk instability) has been predicted recently in |HM01]. This instability 
is likely to move the Fermi level away from the van Hove singularity and it destroys the fiatness of 
the Fermi surface, but it does not destroy perfect nesting. In fact, as it was shown in Section |3.6| , 
perfect p-h nesting will persist unless the particle-hole symmetry is spontaneously broken. Since the 
flat Fermi surface and the van Hove singularity both increase the mixing of p-p and p-h diagram^ 
a Pomeranchuk instability wou ld increase the asymptotic decoupling of the p-p and p-h channels 



(the main result of Section LS. 



^In fact, a flat Fermi surface has an infinity of nesting vectors. This increases the value of non-ladder parquet 
diagrams. The non-uniform density of states at the van Hove singularity focuses on special scattering processes, like 
Fs, which are simultaneously resonant in more than one channel. For a uniform density of states, these processes 
would be negligible. 
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A diverging susceptibility xSd would indicate that there are two different Fermi surfaces for spin 
up and spin down electrons, respectively, with opposite deformations relative to the square Fermi 
surface of the bare system. 

From the particle-hole symmetry of the model, the pairing susceptibilities with a momentum 
close to Q can be treated in the same way as the p-h susceptibilities at a small momentum q ^ 0. 
We therefore investigate the linear response with respect to the operators 77^ (singlet s-wave Cooper 
pairing with a total momentum Q) and 11^ (triplet c?-wave Cooper pairing with a total momentum 
Q). These operators have been introduced in the preceding section. The susceptibilities x-qs ^^'^ 
correspond to limt^q X'^'^'^(0, k), with s-wave and d-wave form factors, respectively. 

Unfortunately, the present investigation has to stay rather qualitative. The reason is that we 
are addressing quantities, which do not feature a leading logarithmic divergence in the perturbation 
theory. For this reason it is not straightforward to calculate them within a one-loop theory. The 
one-loop RG is nevertheless able to give a qualitative estimate of the six susceptibilities XS: Xn, 
XSd J XNdi Xvs ^^'^ XHd in the vicinity of each of the six instabilities SDW, CDW, SF, CF, sSC and 
dSC. In the following section, we only present and discuss the results. The calculations are shown 
in Appendix 



4.5.2 Results 



The results which are established in Appendix ^ (Eqs. ( |G.24| )-( |a29| ), ( |G.38| ), ( |G.42D and ( |G.43D ) 
are summarized in Table 4.1. It turns out that the result does not depend on the form factor (s- 
or d-wave) of the instability, i.e. it is the same for SDW as for SF, etc. 





SDW/SF 


sSC/dSC 


CDW/CF 


XN 





7 





Xs 


? 








XNd 


00 


00 


00 


XSd 








00 


Xvs 







00 


XHd 











Table 4.1: Behavior of various susceptibilities in the vicinity of spin, superconducting or charge 
instabilities, according to the results of Appendix ^. "0" means that the susceptibility is suppressed 
and "00" indicates a divergence. In the two occasions, where "?" appears, there is a competition 
between two terms of opposite sign, which can not be decided within our approach. The sign "-" 
means that from our approach no dramatic behavior is expected at the instability. 



In the cases where a susceptibility is suppressed, I obtain a vanishing susceptibility at the critical 
energy scale, corresponding to a behavior x I log(r — Tc)!"^ (see Appendix^. This can of course 
not be true, since no susceptibility is strictly zero at finite temperature. The real behavior is more 
likdy of the form x{T) ~ exp(— Tc/T). This is probably how one should interpret the RG result 
X — ^ " > 0, namely that x is suppressed at an energy scale given by Tc. 

Quite remarkably, the susceptibility XNd diverges for any of the six instabilities. Whereas the 
suppression of susceptibilities above Tc is very slow (only logarithmic) , XNd diverges even at a higher 
energy scale than T^ according to the results of Appendix |^. This indicates that the square Fermi 
surface is generally unstable with respect to a deformation towards the open Fermi surface of Fig. 
p.7| (Pomeranchuk instability |IIM01 |). As it was argued above, this deformation does not imply 
the destruction of perfect nesting, i.e. it does not hinder the nesting related instabilities SDW,CF, 
etc. In view of Section |3.6| , I suspect that the deformed Fermi surface will increase the decoupling 
of the different p-p and p-h channels, which was anticipated in Section 4.3. 
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Apart from that, we see that the spin instabihties SDW and SF are characterized by a suppressed 
charge compressibihty xn, indicating the opening of a charge gap. This is expected from the 
mean field description of the SDW and SF phases. But it is not trivial that a weak coupling RG 
theory, where no quasi-particle gap is implemented, can reproduce this result. For the uniform 
spin susceptibility, there is a competition between two effects, which presumably cancel each other. 
In fact, the SDW ground state allows for low energy spin excitations and there is no reason for a 
vanishing or diverging spin susceptibility. The numerical works [ HMOC | and | HSFROl | both obtained 
a finite spin susceptibility at the SDW transition. 

In the case of s-wave or d-wave superconductivity, there is a suppressed spin susceptibility xs- 
Within BCS theory, this is due to the fact that the electrons get paired into singlet Cooper pairs. To 
polarize the spins, one has to break these pairs and this needs energy (spin gap). It is remarkable 
that a RG theory of the (unpaired) normal state already yields the suppression of xs- For the 
charge compressibility, I get a competition between two effects. Although the present calculation is 
not precise enough to give a conclusive answer, I suspect that the compressibility is neither zero nor 
divergi n g at the superconducting transition. In the numerical calculations away from half filling 
HMOO , HSFROl I, a strong upturn of the charge compressibility was observed, but only very close 
to the superconducting transition. 

The charge instabilities CDW and CF have vanishing susceptibilities xn and xs- This indicates 
the opening of a charge- and spin gap. More unexpected is the divergence of the susceptibilities xSa 
and Xri- T his suggests that the formation of a charge fiux-phase (or d-density wave as proposed in 
CLMNOl]) could be accompanied by even more exotic phenomena such as ry-pairing. The diverging 
XSd indicates that there could be two Fermi surfaces, one for spin up and one for spin down electrons. 
These points certainly require further investigations before definite conclusions can be drawn. 



Chapter 5 

Conclusion 



The present thesis provides a systematic analysis of the possible instabilities of weakly interacting 
electrons on a 2D square lattice. This is done using a one-loop Wilsonian RG approach, i e. to the 
leading logarithmic order in the energy-scale variable A. 

In the general case where the Fermi surface is neither nested nor at a van Hove singularity, 
the one-loop RG equation becomes equivalent to the ladder approximation and leads to Kohn- 
Luttinger superconductivity at very low energies. Other instabilities do not occur in the weak 
coupling limit. In particular, the inclusion of particle-hole contributions and of certain self-energy 
effects is not consistent within a one-loop calculation. The superconducting instabilities of the 
extended Hubbard model, including a nearest-neighbor interaction, have been studied in Chapter 

i . 

If the Fermi surface passes through a van Hove singularity without being nested, the one-loop 
approximation is probably not sufficient, even for weak coupling. An appropriate calculation for this 
case has still to be invented. It has to include self-energy effects as well as higher-order (two-loop) 
diagrams. Furthermore it may be necessary to use explicitly a finite temperature to regularize the 
theory instead of the band cutoff A, since the two a re equi valent only within leading logarithmic 
precision. A step in this direction has been made in | IIS01a |. 

Most results have been established for the half-filled nearest-neighbor tight-binding model with 
its square Fermi surface. Besides s- and c?-wave superconductivity I have identified commensurate 
density waves and flux phases in both the charge and spin sectors as the dominant instabilities. 
The transition lines of the phase diagram as a function of some parameters of the interaction are 
flxed by exact symmetries and therefore robust for various approximation schemes. 

In addition, the results of Section 4.5 show that the square Fermi surface of the half-filled 
nearest-neighbor tight-binding model is generally unstable with respect to a deformation towards 
an open Fermi surface. This result is in line with the partly numerical work of Halboth and Metzner 
HM01|. It is argued that such a deformation, although it breaks the D4 symmetry of the square 
lattice, cannot destroy perfect nesting (unless the particle-hole symmetry is also spontaneously 
broken). The deformed Fermi surface would then be well described by the anisotropic half-filled 
nearest-neighbor t^-ty model. The deformed Fermi surface will be shifted away from the saddle 
points (the latter are pinned to (0, tt) and (tt, 0) by symmetry) and it is no more a straight line. 
Both features, curvature of the Fermi surface and detuning from the van Hove singularity, contribute 
to an enhanced decoupling of spin, charge and superconducting instabilities. If the deformation is 
weak, such that the saddle point regions still contain a major part of the density of states, the 
instabilities and also the phase diagram will be exactly the same as for the square Fermi surface. 

In the case of a charge-density wave or a charge flux-phase (or d-density wave |CLMN01| ), Ihave 
detected additional diverging susceptibilities, which signal a tendency towards ?7-pairing (pairing 
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into Cooper pairs with momentum (tTjTt)) and spin-dependent Fermi-surface deformations. Further 
studies would be required to give a definite interpretation of these features. 

On a technical level, I have found that the dominant RG flow of the two-particle vertex function 
in the limit of small energies is controlled by scattering processes with momenta close to van Hove 
points. Nevertheless, the low energy effective action contains relevant couplings between electrons 
everywhere near the Fermi surface. The couplings far away from the saddle points diverge at the 
same critical energy and with the same power '--^ (A — Ac)~^ as the couplings at the saddle points. 
From this point of view, there is no sign of a scenario with strong effective couplings at the saddle 
points and weak couplings on the remaining Fermi surface. 

It was shown that the RG equations for different superconducting and density-wave instabilities, 
although strongly coupled at the initial stage, become decoupled in the asymptotic limit of small 
energies. Thus the asymptotic result turns out to be similar to that of a generalized random phase 
approximation (RPA). The decoupling arises because the effective coupling function is strongly 
enhanced only for special configurations of the external momenta, i.e. in a small region of k-space. 
In fact according to our analysis, the diverging part of the effective interaction restricted to the two 
patches Pi and P2 (formulated in terms of creation and annihilation operators) writes 
2 

J2 J2 •^A(q) cjk4-k+qCi-k'+qCTk' 

^ 2 

+ 2 4 H /A(q)5I^^k4'k' + Q+qC<T'k'C^k+Q + q (5.1) 

= l kePi.k'ePj.q CTcr' 



1 ^ 

2^12 9il /A(q)5I^^k4'k' + Q+qC<T'k+Q+qCak' 



i,j = l kGPi,k'GPj,q era' 



where /A(q) is a strongly peaked function such that |q| < VA, but k and k' arc within patches of size 
p. (The effective A-dependent coupling constants are related to the vertices introduced in Section 
4.2.3| by g^j = — r*j-) The strongly peaked nature of /A(q) includes that the effective interaction 



gets long ranged in real space. If instead we would insist on a local interaction, /A(q) would b e 
a constant and Eq. ( ^?l| ) is simplified to the model considered by H. J. Schulz | Sch87 , 3ch89| . 



The non-locality emerges here from the unbiased treatment of the most diverging terms in the RG 
equations. 

In a one-dimensional system one never generates non local effective interactions. The reason is 
that the low energy excitations of a one-dimensional electron gas are constrained to two privileged 
points in k-space: the Fermi points. This special geometry allows for a strong mixing between 
the various interaction channels, because reducing the arguments of the gA(fci, A;2, fca) function 
to the Fermi points fixes at once the total momentum and momentum transfer. Consequently, 
superconducting and density wave instabilities remain coupled in one dimension. It is however not 
unusual to encounter non-locality in the vertex function rA(fci, ^2, ^3) where CO^i) > A is allowed. 

The saddle points of the two-dimensional dispersion have a different status than the Fermi points 
of the one-dimensional electron gas. They are privileged only due to the diverging density of states, 
but the low energy excitations exist on the whole Fermi surface. As a consequence, the effective 
interaction (5a(Aij ^2, ^3) = — rA(Ai: ^2, ^3)1 where ^(k^) < A) depends on k points which can 
move continuously on the Fermi surface. The high sensitivity of the RG equations on the external 
momenta leads to a non local effective interaction at very low energies. 

In spite of the asymptotic decoupling, our RG flow develops enhanced correlations in the rf-wave 
superconducting sector, although they are dominated by the SDW at half filling. This supports 



the by now well established result of Zanchi and Schulz |ZSOO| that d-wave superconductivity will 



occur in a slightly doped system, where the density wave fluctuations are suppressed. 
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Our analysis has nevertheless shown that a discretization of the Fermi surface in terms of 
a finite number of patches can enhance artificially the coupling between the different scattering 
channels in the low energy regime. In fact, in the numerical studies ^S0(\ , |HMOO| , |HSFROl| different 
susceptibilities are found to diverge at a single energy scale. By contrast, the present results show 
that only the susceptibility of the dominant instability diverges. This behavior has been referred 
to as the "moving pole solution" by the Russian school |DK72, GD74|. 

The decoupling of the RG equations admittedly has only been established to leading logarithmic 
order in the energy cutoff A. Subleading terms of the one-loop RG equation which would couple the 
competing channels, are neglected. It is argued however that a consistent treatment of subleading 
terms requires to go beyond the one-loop approximation. 

In order to estimate how small the bare interaction must be, we recall that our approach requires 
the patch around the saddle point to be small compared to the size of the Brillouin zone (p <^ 7r/2). 
On the other hand, the density of states of this patchQ has to be big as compared to that of the 
remaining part of the Brillouin zone (i.e. log4/p^ <C log4p^/A). The bare interaction g ^ U must 
be small enough such that 

2C/i?,>,0)^^log^^^l. 

(The factor 2 appears because there are two saddle point patches in the Brillouin zone.) If for 
each of the two "<C" signs above, a factor of 10 is introduced, this amounts to U ^ 0.02 i. A less 
stringent factor of 3 for each of the two inequalities would correspond to U ^ 0.6 i. 

To conclude, we have seen that the weak coupling analysis of the 2D Hubbard model reveals two 
features which are common to the high-Tf, cuprates. Namely the SDW phase at half filling, which 
is continuously transformed into the antiferromagnetic Mott insulator if the coupling is increased, 
and the existence of a c?-wave superconducting phase in the doped regime. Unfortunately we have 
gained no better understanding neither of the under-doped regime nor of the unusual "normal- 
state" properties at optimal doping. In fact, strong- or intermediate coupling seems to be a key 
ingredient in a theory for these materials. 

It is nevertheless possible that a two-loop calculation would change the picture. The self-energy 
of the half-filled Hubbard-modcl docs not satisfy the requirements of Fermi liquid theory above the 
SDW-critical temperature |Riv02], which indicates that self-energy corrections arc important. But 
a consistent two-loop calculation for a 2D Fermi system is a very difficult task. 



iSee Eqs. (15) and <^^. 
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Appendix A 

Time reversal invariance 



In quantum mechanics, the anti-unitary time reversal operator T maps wave functions to their 
complex conjugate. In a single-band model, it is reasonable to assume that the Wannier wave 
functions are real, because otherwise there would be two degenerate wave functions per lattice site. 
The Wannier states are therefore symmetric under time reversal (i. e. T|r) — |r)). Including the 
spin degree of freedom, we find T\ra) — a\r — a), where the number 1 (—1) is assigned to t (i)- 
This implies 

T c^T-' = a 4_,, T Cr,r-i = a c,_, (A.l) 
and finally after Fourier transformation 

^cL^"' = ^c^-k-.- Tc^aT-^ = ^C-^-a- (A.2) 

Suppose that the system is symmetric under time reversal, i.e. that H commutes with T. In 
order to investigate the consequence of this symmetry for thermal averages, recall that {Til}\T(j)) = 
{T^\T\4>) = {il}\4i) ■ Furthermore, if \n) is an orthonormal basis of the Fock space, then is T\n). 
One then shows 



(0> = ^Y.^n\e-^'"'''''0\n) 

n 



1 

Z ^ 



{TO^T^^). (A.3) 



This equation, applied to the two particle Green's function, leads to Eq. 2.27 . Note that the 
hermitian conjugate of c(r) is c\—t). 
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Appendix B 

Parquet diagrams: An example 



As an example of a typical parquet diagram, I consider the diagram shown in Fig. [B.1[ For 
simplicity, I have chosen the total momentum and the external frequencies fco and fcg to be zero and 
I assume a constant coupling function g. The value of the diagram (at T 0) is then 



I3L 



Po 



-y 



«?)o->-|Cpl 



(B.l) 



The sum is over all momenta p allowed by the infrared cutoff, i.e. satisfying |^p| > A. The 
value of the internal p-h bubble is 



r 2 



1 ^, e(eq)~e(|^p| + eq+p+k) 

^'V Cq-l^pl-Cq+P + k ' 



(B.2) 



where the sum over q is restricted to momenta such that |^q| , |^q+p+k| > A due to the infrared cutoff. 
The p-h bubble gives a leading order contribution, if p + k is a nesting vector. For simplicity, and in 
order to consider a case where the p-h contribution is most relevant, I assume that ^q+p+k = — ^q 
for a non negligible set of p-vectors. Within this assumption. 



ipo^-|Cp| 



A 



2e + iepi 



(B.3) 



PA 



q+p+k 



-k' 



Figure B.l: Example of a two loop parquet diagram. 
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The parquet diagram shown in Fig. p3.1| then equals 



Ja t. 



2C' + e 



(B.4) 



where is the density of states of the momenta p, such that p + k is a nesting vector. 

To the leading logarithmic order and in the absence of van Hove singularities, one can replace 
the density of states t'(^) and z/(±f ') by their value at the Fermi level, so the diagram of Fig. ^ 
is proportional to 



P(A) 



^de 

A ^ 
W 



w 



A 2C' + e 



w 



A 



d£ 1 

^ 2 °^?(2A + a 



A^O 1 



(B.5) 
(B.6) 



The diagram is indeed of the leading order log^ A. 

To get the contribution to the RG equation, we take the derivative with respect to A 



1 AW^ - 



A2 



w 



C(2A + e) 



(B.7) 



The second 
' 1/A and is 



The first term is obtained by putting the large loop variable f to the energy scale A. 
term comes from the derivative of the internal p-h loop. It is equal to j-^ log ^ 
thus of subleading order, as A ^ 0. 

This result can be generalized to the following statement. If a diagram is reducible in one channel, 
the dominant contributions to its A-derivative will be obtained by deriving only the propagators 
connecting the irreducible blocks and not the irreducible blocks themselves. This is what is used in 
the derivation of the one-loop flow equations in section 

The same is true in the case of a van Hove singularity. To consider the case, wh ere b oth, v and u 
are logarithmically diverging at the Fermi level, we multiply the integrand of Eq. B.5 by log^ log^' 
and repeat the same arguments. In this case, the leading term is given by log^ A. 



Appendix C 

Proof of the Polchinski equation 



Starting from the definition of the effective interaction, one obtains 



17=0 



For more details about the calculus with Grassmann variables, see ^al99| . The derivative of the 
last equation gives 



-WaM 



This is the Polchinski equation 3.98 
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Appendix D 



Estimate of the self-energy in 
second-order perturbation theory 



In order to decide whether self-energy effects are relevant within the leading-order approximation, 
I estimate the self-energy corrections to the quasi-particle weight z — {1 — dik^'^Aik))"^ and to the 
Fermi velocity vp = Vk(^k + S(0, k)), within perturbation theory. For simplicity, I assume a local 
interaction, i.e. a momentum-independent coupling g, but the result remains true for a general 
short ranged interaction. 



The first fc-dependent self-energy contribution comes from the second-order diagram Fig. D.l 
The interaction lines have been reduced to points, since I assume a local interaction. The analytic 
expression i^ 

= -.9'^E'CA('z)CA(g')CA(g + g'-fc) (D.l) 



2 1 , nqnq>(l - Tlq+q'-k) + (1 - "q)(l - ?^q')nq+q/_k 



where Up — {1 + e^^^) ^ is the Fermi occupation distribution. The sum in Eq. D.2 is restricted to 
the momenta allowed by the infra-red cutoff |^q|, |^q'|, |^q+q'_k| > A. 



D.l Quasi-particle weight 

Differentiating Eq. |D.2| with respect to ko leads to 



^.fcoSf (fc) 



ik(s—0 



nqriq/(l - riq+q'-k) (l-nq)(l 



(Cq+q'-k - ■Cq - ?q')^ 



(Cq-(-q'-k - Cq - Cq')^ 



(D.3) 



We are going to make a rather crude estimate of this expression. The Fermi function rtp reduces 
to a step function 6(— ^p) at zero temperature. It is then easy to check, that the expression 
^q+q'_k — ^q — ^q' in thc denominator is actually a sum of three terms of the same sign. Hence 



^In fact the diagram shown in Fig. |D.1| stands for two distinct diagrams with wavy interaction lines, one with a 
Fermion loop and the other without a Fermion loop. A factor of —2 is associated to the Fermion loop, hence the 
minus sign in Eq. 



D.l 
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q+q'-k 



Figure D.l: Second-order self-energy diagram. 



we can use (Cq+q'-k ^ Cq ^ Cq')^ — (Cq + ?q')^ to estimate Eq. [D.3| . Furthermore nq_|_q'_k and 
1 — ?^q+q'-k are both bounded by 1, thus 



^.fcoSf (^) 



ikn—O 



< 



(Cq + ^q' 



_ , (l-nq)(l~nq/) 
)^ (eq+Cq')^ 



(D.4) 



The estimates used so far are independent of the dispersion ^p. In the case of a van Hove 
singularity, where the density of states is logarithmically diverging, one has 



i/co— 



< const • g 

< const • g'' log 



' dd d^ 

J A J A 



log^log|a 



.2,„ J 



A 



(D.5) 
(D.6) 



D.2 Fermi velocity 

at zero 



For an estimate of the self-energy corrections to the Fermi velocity, we evaluate Eq. D.2 
frequency and zero Temperature, where it can be written as 

V^(2)^n U^ _ „2 1 Y^, "q"q'(l- "q+q'-k) -(l-»q)(l-«q')nq+q'-k „x 

Kq+q'-k| + |4q| + IWl 

We consider the nearest-neighbor tight-binding band and assume that k is close to one of the 
saddle points (Pi — (tt, 0) and P2 = (0,7r)), say k = Pi -I- k, where k is small. Furthermore from 
the sum Sq,q', we take into account only the contribution of momenta which are close to one of 
the two saddle points, because this is where the leading terms are expected to come from. I.e. 

q = Pi -I- q or q = P2 -I- q, where |q| is smaller than some patch size p and similarly for q'. Since 

(2) 

each of the two momenta are summed over two patches, E)^ (0, k) is a sum of four terms. 



A 

'^^(Ok)wg^^— '*q'^q'(l ^ "q+q'-k) - (1 - '^q)!! - '^qO'^q-t-q'-k 



ICq+q'-kl + ICql + ICq'l 



q=P.+q,q'=Pj+q', 

(D.8) 

Because of the exact nesting ^p^+p = — ^P2+p, the denominator is the same for all four terms, but 
the numerator changes according to npj^+p = 1 — np^+p. Summing the four contributions, one gets 

Ef(0,Pl-Hk)«g2ij y «ign(lq+q'-k)Bign(eqjqO (j3g) 
^^q,t^<p Wkl + leql + leq'l 



D.2. FERMI VELOCITY 



87 



where I have omitted the tildes on the momenta and introduced the abbreviation — ^p^+p 
instead. 

When taking the gradient of Eq. D.9, we have to remember the imphcit infra-red cutoff re- 
striction |Cq+q'-k| > A. The sign function sign(^q_|_q/_k) in Eq. D.9 has to be replaced by 
sign(^q_i_q/_k)0(|^q_i_q/_i(;| ^ ^) : i-S- a fuuctiou with two steps at '^q+q'-k = iA. Finally one 
finds 



VkSf (Pi +k) « .g2 Y.' sign(CqCq' ) 



q+q'-k| 



A) 



edlq+q'-kl-A) 

A + ICql + llq' I (llq+q'-kl + ICql + l^q' 



IqMq'Kp 

Because of the hyperbolic form of the dispersion near the saddle points, we have 

VkCk-q-q' = Vki^k ~ Vq^q — Vq'i^q'. 



VkCk-q- 

(D.IO) 
(D.ll) 



It is clear that we can neglect Vq^q and Vq'^q' for the leading logarithmic contribution to Eq. 
D.IO, since Vq^q vanishes at the saddle point. 



The first term in the parenthesis of Eq. D.10| is bounded as follows. 



r4 5Z 



IqMq'Kp 



, '^(llq+q'-kl - A) 

A+||q| + ||q'| 



- y 



^(llp|-A) 



|qMp|<p 



A+leql + ICp+k-ql 



< 



- y 



, -^dlpl-A) 



^lqi1;r<p^+i^<^i 



< const. • log — 
A 



A /-^f) los: — 
Ao t ^^ iog ^ 



A + e 



< const. • log' 



Ao 
A ' 



(D.12) 
(D.13) 



where I have introduced the summation variabl e p = q + q' — k (the approximation consists in 
replacing the constraint |q'| < p by |p| < p). In D.12 , I use the fact, that the density of states at 
the saddle point is proportional to log ^ . 



The second term in the parenthesis of Eq. |D.10| can be estimated in a similar way as we did in 

- y 



Eq. D.3 



|q|,|q'l<p 



e(|Cq+q'-k| - A) , 1 3A0 

^ ^ < const. • log — 



(|eq+q-k| + |^q| + |eq'|)^ 



(D.14) 



The estimate 4.2 follows immediately from Eqs. D.IO, D.ll, D.13 and D.14 
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Appendix E 

Analytic expression of the bubbles 
for momenta parallel to the Fermi 
surface 



Here I give an approximate analytic expression for the bubbles Bpp{A, k) and BP^{A, k) (defined in 
Eq. |3.102 ) at a momentum k — {k, k) (where k g [0, tt]) parallel to the square Fermi surface of the 
nearest-neighbor tight-binding band. Since BP'^'iA, k) = -BPP{A, Q - k), we only need to consider 



the p-p bubble, the quantity which is plotted in Fig. 4.5. 

The calculations are rather tedious, using the variables x± — tan for the integration over 

the Brillouin zone and distinguishing many cases. For the convenience of notation we introduce the 
variables 

u = tan k/2 

and 

V — tan k/ 4 

and we multiply the value of the cutoff by four [i.e. the formulas are given for Bpp{AA, k, k) instead 
of BPP{A, K, k)]. Finally for an approximate solution of the integrals, a parameter e was introduced, 
which is supposed to be 

A < e < 1. 



Case 1: < u < A, 
i. e. K very small. 

BPP(4A,k,k) = 



(27r)24A 



log 



A2 



A 

?log 



x/A2 - 
~i? + 1 - e 



arccot 



ue + u 



^(A2 -it2)(e2 _ A2) 



(VT+z? - 1 + e + u)(Vl + u2 + v-u) 



Case 2: A < u and v < 1 — e, 

i. e. K not too close to and tt. 



BPP{AA,k,k) = 



1 



(27r)24A 



log 



Ve^ - A2 + . 



A 



A2 



log 



Ve^ - A2Vu2 - A2 - ' 



{e + u)A 
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ANALYTIC EXPRESSION.. 



+ 



:l0g 



(Vl + m2 + 1 - e - + u2 - w + u) 



(Vl + - 1 + e + + + - w) 



Case 3: l-e<i;<l-A 
i. e. K close to tt. 

Sff(4A,K,«;) = 



1 + 



(27r)24A 



+ 



log 



y/{l-Vy - + 1-V 



A 



log 



V'(l - vy - AVu^ - A^ - «(1 - - A2 



(1 -t; + u)A 



Case 4: 1 - A < < 1. 

BPP(4A,k,k) = 0. 

(The approximate expression for Bpp{4A, k, k) given here vanishes at = 1 — A. However it is 
relatively easy to establish that the exact quantity B'pp^AA, k, k) vanishes at Vc = ^j^- The 
difference is due to the approximation used to solve the integrals.) 



The asymptotic behavior for k ^ 1 can be extracted from the expressions given in the cases 1 and 



Appendix F 

Thermodynamic relations for the 
charge compressibility 



Consider a thermodynamic system of N identical particles confined in a volume V in thermal 
equilibrium with a heat-bath of temperature T. The free energy is some function F{N,V,T). 
Because N, V, and F are extensive quantities, the free energy can be written as 



where fJ,{N, V, T) ~ dF/dN is the chemical potential and p{N, V, T) — —dF/dV is pressure. 
The isothermal compressibility is defined by 



(F.l) 



1 



(F.2) 



To relate this quantity to the uniform charge susceptibility xn = y (^) ' take the partial 
derivative of Eq. F.l with respect to N and obtain 



= 



V 



dp 
dN 



(F.3) 



Similarly we take the partial derivative of Eq. p.l with respect to V to obtain 

= N 



dV 



- V 



N.T 



dp 
dV 



(F.4) 



N,T 



The second term of Eq. p.3| and the first term of Eq. F.4 are related by the Maxwell relation 



dp \ 
dN J 



V,T 



av 



The relation 

V (dN 

is easily obtained from Eqs. F.3| , F.4 and F.l. 



N.T 



V 



(F.5) 



(F.6) 
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THERMODYNAMIC RELATIONS. 



The inverse compressibility is given, up to a multiplying factor by the second derivative of the 
free energy with respect to the particle number 

A diverging compressibility is thus a signature for phase separation (the system separates into two 
phases, if the free energy seizes to be a convex function of N) , whereas a vanishing compressibility 
signals a cusp of F as a function of the particle number. 

As an example, I calculate the charge compressibility of a system of non-interacting electrons. 
There, 

iV(M,y,T)=^2n(a), (F.8) 

k 

where n{C) is the Fermi distribution function. One immediately finds 

XN^-2j de HO^^^^ 2K0). (F.9) 



In the case, where the density of states diverges logarithmically at the Fermi level, Eq. p.9| is 
replaced by 

XN ' 2i.(r). (F.IO) 



Appendix G 

Calculation of the "uniform 
susceptibilities" 



The term "uniform susceptibilities" is precisely defined in Section 4.5 



itative behavior of each of them is predicted as a superconducting, 
instability is approached. 



In this appendix, the qual- 
density-wave, or flux-phase 



G.l Relation between uniform susceptibilities and the ver- 
tex function 

The uniform spin and charge susceptibilities at small momentum depend on the vertex function 
for forward- and exchange scattering. It is necessary to distinguish the so-called q- and w- limits. 
They are defined as 

r{{p,p') = hm r{p,p' + 1,0=0 (G-i) 

and 

ri(p,p') = lim r(p,p' + 9,p')lq=o ■ (G.2) 

The exchange vertices are defined as = -^r^ and FJ^ — Xr£. Similarly, we define two limits for 
the vertex function F**, which describes scattering of particle pairs with total momentum close to 
Q: 

r^(p,p')= liniF(p,g-p,g-p') 1,0=0. (G.3) 



p,q 



Q 



(G.4) 



The RG procedure with a cutoff A in the band energy yields F^^ rather than Fq. The reason is 
that at zero temperature, p-h pairs can only be created with a momentum q big enough, such that 
|v^- -ql > A, where Vi;- is the Fermi velocity (i.e. the gradient of ^ somewhere on the Fermi surface)0. 
Therefore, if not specified differently, I always refer to the co- limit in this thesis. 



iSee also CS95 or [Dup98 
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The two limits are related by the following Bethe-Salpeter equation^ 

KiP^P) = KiP,p') - ^J2^Zip,k)Dl''{knk,p'), (G.5) 

k 

where = -T", 2rf - or P', and DPl'ik) = C{k)C{k + g)|^„^o- The limit q ^ is to be taken 
after the summation over /cq- 

From Eqs. (3.61) - (3.63), we obtain the following identities for the abovementioned spin, charge 
and pairing susceptibilities 

^ = -Jfl E /'(P)^q' W + Tgjlyl E f{p)Dl\p)VX{py)Df{p')f{p'), (G.6) 

where /(p) is the form factor, i.e. / = 1 (for xs, Xn or Xns) ^ d-wave function (for xSdJ X-ZVd or 
XUd)- The vertex function is = — F^ (for xs or xSd)i 2F'^ — F"^ (for xn or xat^) or F'' (for x^s or 
XnJ. 

There is no infrared cutoff in Eqs. ( |G.E1 ) and (|G.6| ), but a finite temperature T. Nevertheless, 
we use the results of the zero temperature RG analysis and replace the cutoff by the temperature, 
since both are equivalent in the leading logarithmic approximation. 

We will neglect the dependence of F* (A;, k') on the frequencies fco and /cq from now on. A simple 
expression for x is readily obtained, if we manage to write the vertex functions F^ in the form 

F^(k,k') = ^7:e:;(k)e:(k'), (G.7) 

n 

where the functions (k) satisfy the orthogonality relation 

- ^ E <(P)em(P) = ^nn.. (G.8) 

p 

The e* (k) and 7^ are eigenfunctions and eigenvalues respectively of the operator defined by 

h{].)^-^Y.^Z(k,p)Dl'^{p)h{p), (G.9) 

p 

where h is any function of k. Summing the Matsubara frequencies po and taking the limit q — > 0, 
Eq. ( p.9| ) becomes 

Mk) ^--l^F^(k,p)n'(Cp)/i(p), (G.IO) 
p 

where n'(^) = dn/d^ is the derivative of the Fermi distribution function. 



Using the relations ( G.7 ) and ( G.8 ), the Bethe-Salpeter equation ( G.5 ) yields 

F^(p, p') = T^e„(p)e„(p'). (G.ll) 

We expand the form factor / in terms of the eigenfunctions as 

/(k) = ^aX(k), (G.12) 

n 

^Eq. (G.5) is nothing else than Eq. ( |3.19| ) (in the p-h channels) or (3.5) (in the p-p channel). They have an 
identical form because of particle-hole symmetry. plays the role of the irreducible vertex, since the lu- limit of the 
p-h bubble vanishes. 
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where are coefficients, and obtain the foUowing result for the susceptibihties 



If the form factor /(k) is chosen proportional to one of the eigenfunctions, the susceptibility will 
be given by 




where is the value of the susceptibility if the four point vertex is put to zero|^ 

The eigenvalue 7 depends on all the parameters of the system including temperature. At high 



temperature and weak coupling, the eigenvalues 7 are small numbers. In that case, Eq. (G.14) 
gives a small correction to the bare susceptibility. A positive value of 7 increases the susceptibility 
and a negative value suppresses it. The susceptibility diverges if 7 approaches unity. 

In the opposite case, if 7 ^ —00, the susceptibility is completely suppressed. One should be 
careful with this conclusion. It can not be justified within perturbation theory, since the infinite 
series (1 — 7)"^ = 1 + 7 + 7^^ + • • ■ does not converge for 7 < —1. The question is whether the 



Bethe-Salpeter equation (G.5) is non-perturbatively true (i.e. if is not small) or not. Although 
I am not aware of a non perturbative derivation of the Bethe-Salpeter equation, I will assume that 
X is (at least qualitatively) suppressed if 7 — )■ — cx). 

In the case /(k) = 1 and approximating the vertex by a constant r^(k, k') —U, one obtains 



the random phase approximation (RPA), where 7 = — x'-^-'J/. However the relation (G.14) goes 
beyond the RPA. It is exact apart from neglected self-energy effects and the neglected dependence 
of the vertex T^{k, k') on the two frequencies and k'^. The difficult part is of course to calculate 
the vertex functions and to write them in the form of Eq. ( p.7| ). In the following Sections, we show 
how our previous result on the behavior of the vertex near an instability can be used to predict the 
behavior of the susceptibilities. 



G.2 Diverging peaks in the vertex function 

In the spirit of Section f4.3| , I propose the following approximations. (The uj- index of the vertices 
are omitted from now on.) 

1. Due to the diverging density of states near the saddle points, only the momenta within the 
saddle point patches Pi and P2 are considered. 

2. The vertex functions r^(k, k') vary slowly with k and k' except close to the special configu- 
rations k -I- k' = and k — k' = Q, where they may develop strong peaks. 

For the special configurations k + k' = and k — k' — Q, each of the vertex functions , 
and F'' coincides with either F'', F"^ or F-^*-^"^. In factQ, 

F''(p, -p) = F-(p, -p), F''(p, p + Q) = F'^(p, -p), 
F-(p,-p)=F^^^^(p,-p), F-(p,p + Q)=F''(p,p), (G.15) 
F/(p,-p) = F^g^(p ,p), F/(p,p + Q) = F-(p,p). 

3-)^{0) jjiay in general contain self-energy corrections to the single-electron Green's function. In our approach, 
where self-energy corrections are neglected, x'*^' is the non-interacting susceptibility. 
■^See Appendix M for the notation. 
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For example the function r-''(k, k') shows the following behavior 

(G.16) 



\ r^ + r^,k-k' for kePi, k' 



eP2 



where T^^^, ^ develops a diverging peak at k + k' = in the case of a superconducting instability 
(s- or d-wave). The behavior of the peak close to the instability is given by 

ri^k+k' ~ 7777772' (G-17) 

where a and b are positive constantsn. This function has a strongly divergent maximum but is 
completely negligible for |k + k'| > ^/T. Similarly, the function Tq-q 2 diverges in the case of a 
spin instability (SDW or SF) as 

r2,q« 72- (G.18) 



T-Tr. 



q Q 

T 



We have thus separated two parts of the forward-scattering vertex function, the regular part (ap- 
proximated by the constants Fj and f{) and the possible peak close to the instability. 
Correspondingly, the vertex functions F"^ and F'' are approximated by 



and 



T^(l,^')^J fork,k'ePi 

i ik,k j ~ <; j^v^Yf k_k' for k e Pi, k' e P2 ' ^ ^ ^ 



where F5^_^, ^ behaves, close to a spin or charge instability, as 



j.d ^ a r ;SDWorSF 

^•'^-'^'^ , ./k-k'-QA^ ■ I 1 ;CDWorCF ■ ^^-^^^ 



On the basis of these approximations, two approximate eigenfunctions of the operator (G.IO) 
are given by 

s-wave: /(k) = 1 and d-wave: /(k) = ( _} foJ k ^ J • (G-22) 



2 



The approximate eigenvalues are given by 

_ - Ell EpgP. r^(k, P)n'(ep) /(P) 
^ /(k) 



(G.23) 



These eigenvalues are only approximate because the right-hand side of Eq. ( |G.23 ) is only approxi- 
mately independent of k. 

For example, to calculate the spin susceptibilities xs and xSdi '^6 consider the eigenvalue equa- 
tion for the operator ( p.lO| ), where F* = -F^ Using Eqs. ( |G.23| ) and ( p.l9| ), we obtain 

^This is the result of the ladder approximation at finite temperature (see Section 4.3.5). 



G.3. RENORMALIZATION OF THE FORWARD-, EXCHANGE- AND r]-VERTICES 
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pePi 


"'(Cp) r^^k+p ^ Yp. 


E "'(Cp)rlk-p- 

P6-P2 


(G.24) 


Analogous expressions are obtained for = 


2Vf - r^, namely 







2Ti^Tl±{2Ti-T'{) 



-y 

pePi 



-'(ep)ri:??pT^ E"'(^p) (2r; 



2, k-p 



r?,k-p) 



(G.25) 



peP2 



and for = T'', 




" 12 E "'(ep)r2,k+p+Q T E "'(^p)rlk-p- (G.26) 

pePi peP2 


The second and third terms 
ities very directly. Using Eqs. ( 


on t^ 
G.17 


le ] 

), ( 


"ight-hand s 
G.18) and ( 


ide 
G.21 


f these equations feel the different instabil 
) one obtains 



sSC or dSC instability: 


J2 E ^'i^p)^i,'i.+P 
pePi 


^ —00, 


(G.27) 


CDW or CF instability: 


^E-'(^p)rlk-p - 
PGP2 


-!■ —00, 


(G.28) 


SDW or SF instability: 


peP2 ^ i:k p 


— > +00. 


(G.29) 



The divergence is always logarithmic ^ log(r — Tc). In the case of a SDW or SF instability, F^ 

diverges twice as fa st as F* ^, suc h that F^ — 2T'^ — F^ remains finite. 

The equations ( G.24 )-( G.29| ) explain already most of the results mentioned in Section 4.5.2 . 
To ob ta in the full result, we still have to investigate the first terms of the right-hand side of Eqs. 
( |G.24 )-( G.26 ). The question is thus how the vertic es F f , F^, f{, F!j, f{ and F| behave as one of 
the instabilities is approached. We argued in Section 4.3 that general vertices are not dramatically 
influenced by the instability, i.e. they do not diverge. It turns out that the forward-, exchange- and 
77- vertices are not general in that sense but that they respond to F'', F^ and F^'-^'^ via a relatively 
subtle effect. 



G.3 Renormalization of the forward-, exchange- and 77- vertices 

The values of the vertices F'f, F^, Fj , . . . depend on the energy scale and are determined by the 



leading-order RG equations. We will use the same formalism as in Section 4.3, where the energy 
scale is represented by a sharp infra-red cutoff A. In the end of the calculation we identify, within 
logarithmic accuracy, the vertex at zero temperature in the presence of a cutoff A with the vertex 
at temperature T — A (without cutoff). 

As an example we choose two momenta k G Pi and k' e P2 with |k — k' — Q| > -s/A and 
consider the RG equation of the vertex F(k, k', k') = F2 . It is renormalized mainly from the p-h 1 
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p^^ vp + k'-k 



Figure G.l: The p-h 1 contribution to dAT-^O^,^')- 



diagram, since the exchanged momentum transfer k' — k is in the neighborhood of Q. The other 
contributions p-p and p-h 2 are neghgible 
The p-h 1 diagram is shown in Fig 



G.l 



If the internal momentum p approaches k 4- Q, both 
vertices of the diagram approach simultaneously the special vertex F^. This is in contrast to the 
case of a vertex with general external momenta (or the case of Fig. 49). Similarly, if p approaches 
— k', both vertices of the diagram approach simultaneously Tf-^^ . 



More quantitatively, the RG equation for Fj reads 



^Ar2^ = i^E ^k'-k(p) (r(k',p,k)f 



(G.30) 



where 



^f(p) 



dpo 

ei-CpCp+q) 



0(ICpl-A)e(|ep+q|-A) 



ICpI + ICp+ql 



(G.31) 
(G.32) 

d\D is zero except 



is the p-h pair propagator after the frequency-integration. Its derivative D 
on the one-dimensional curves |^p|, |Cp+q| — A, where it is positive. 

In the case of a charge- or spin- instability, the vertex F(k', p, k) develops a peak at p = k + Q, 
where it is equal to Fg (see Section 4.3.2). In order to estimate the contribution of this peak to 
the integral in ([G.30 ), we proceed as in Section 4.3.5 

of the order d/t (F(A, k, k)) , where F(A, q) is given by Eq, 
integral diverges like (A — Ac)^^. Similarly, in the case of superconductivity, the vertex F(k',p,k) 
develops a peak close to p = — k', where it approaches Tf-^^ . In both cases, we get 



The singular contribution to the diagram is 
It can be checked, that this 



4.36 



9aF{ 



1 



A- A, 
log(A - 



A.) 



(G.33) 
(G.34) 



We now use this kind of analysis to estimate the behavior of F| ± F'j , the quantity which enters 
Eq. (G.24). First, we choose k, k' e Pi with |k + k'| > VA and write the RG equations for 
F4 — F(k, k',k) and F^ = F(k, Q — k',k). F| is renormalized by the p-p diagram and Ff by the 
p-h 2 diagrams. The RG equation is 



^E^£%'(p)r(k',k,p)F(k,k',p) 
p 
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^J2 E ^Q-k-k'(p) r(p, Q - k', k)(2 - 2X)r{p, Q - k', k) (G.35) 
p 

^ E ^r+k' (p) [r(k', k, p)r(k, k', p) ± r(p, q - k', k)(2 - 2X)r{p, q - k', k)] , 



where DPP{p) is the frequency-integrated p-p pair propagator, Xr(pi,p2,P3) — r(p2,pi,P3) and 



we have used the nesting property DPP{p) — ~'D^_ (p). The vertex terms in the parenthesis have 



Q-qV 



two possible peaks. One at p = k+Q, where the term in the parenthesis equals T^T^±2T^(r^ — Tf) 
and one at p = k' + Q, where it equals T'lr^±2rf^^{Tf^^-Tf^^). As one can check, the function 
^k+k'(P) same value at both peaks. We can thus simply add the two contributions and 

write formally 



dA{n + ri) ^ -J 2(r^f 



ri) J 2r^ (2rf - r^) , 



(G.36) 
(G.37) 



where / . . . indicates integration over a peak. If the peak diverges as (A— Ac)~^, the integral diverges 
as (A — Ac)"^ (and if the peak diverges as (A — Ac)~^, the integral diverges only logarithmically, 
as pointed out in Section 4.3.5). We therefore obtain 



SDW or SF instability: 



log(A-A,) 



(G.38) 



whereas r| ± Tf remains finite in every other case^. 

There is a competition between the first and the last term of Eq. ( |G.24 ). Both terms diverge 
with opposite signs and our calculation is not precise enough to decide even qualitatively on the 
behavior of 75, as the SDW or SF is approached. Physically, one suspects that the uniform spin 
susceptibility is neither suppressed nor diverging in the SDW or SF phase, so that 75 should 
be a small number. I therefore assume that the two competing terms of Eq. (G.24) essentially 
compensate each other. But this remains of course to be shown. 

The analogous investigation of the quantities ± Fj, which enter Eq. (G.25), leads to 



9a(f{±f{)^-15]I)£%,(p) [F(k',k,pf Tr(Q-k',p,kf] 



(G.39) 



where k, k' G Pi with k + k' > \/A- The vertex f{ is renormalized by the p-p diagram, and Fj by 
p-h 1. We have used f{ = F(k,k',k'), F^ = F(Q - k',k, k) and I)^''_^„^,(p) = -i>£%,(p). The 
expression in the parenthesis on the right-hand side of Eq. ( G.39 ) has two possible peaks. One at 
p = k + Q, where the term in the parenthesis equals (Fg)^ =F (^3)^ and p = k' -f Q where it equals 
(rf)^ T (rf*"'^)^- As before, the prop agator I)^^, ^, (p) has the same value at both peaks. Adding 



the two contributions and using Eqs. ( |G.36D and ( |G.37D , one obtains 



5a 



9a 



2ri 



2ri - Ti 



2ri - Fl - 2F^ 



{2ri- 



3(F^)V2(Ff^^)^ 



We conclude that 



^Note that in the case of the spin instabilities SDW or SF, diverges twice as fast as F'*, such that F'^ : 
is finite. 



(G.40) 
(G.41) 
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dSC or sSC instability: 


2r{ - 




— > ~oo, 


(G.42) 


for any instability: 


2r{ - n - 


2r^ + r^ — . 


+00. 


(G.43) 



Because of Eqs. (G.42), (G.27) and (G.25), there is a competition between two terms contribut- 
ing to 7Ar in the case of superconductivity. The situation is similar to the spin instability discussed 
above. On physical grounds, the charge compressibility is neither suppressed nor diverging in the 
superconducting state, so there is a chance that the competing terms cancel each other and 7jv 
remains small. In contrast, the terms contributing to (Eq. ( G.25|) ) are never in competition 
but they reinforce each other and lead to a divergence of xn^ ^-ny instability is appro ached. 

An analogous analysis shows that the vertices and which enter Eq. ( G.26 ) are never 
divergent. 



Appendix H 

Overview of notations for the 
vertices 



This appendix gives an overview of the different notations for various vertices, which have been 
used in chapter ^. 

The main object in this thesis is the IPI vertex function r(A:i, fc2, fcs), which was introduced in 
Section 2.1.2| . Starting from this object, six different vertex functions of two momenta are defined 



as foUows 



r^^^(fc,/fc') 




(H.l) 




= T{k,k' + Q,k'), 


(H.2) 


r^{k,k') 


= T{k,k' + Q,k + Q), 


(H.3) 


r''{k,k') 


= lim T{k,q-k,q~k')\ a ^ 


(H.4) 


Tf{k,k') 


= hm r(/c,fc' + g,fc')lq=o, 


(H.5) 


r^ik.k') 


= hm T{k,k' + q,k + q)\ 

9o^U 


(H.6) 



where Q = (0, Q) = (0, tt, tt). The last three vertex functions sometimes carry the index uj to precise 
that the frequency is sent to zero after the momentum q has been set. The opposite succession of 
the hmits is denoted by an index q. If no index is given, I always refer to the w-limit. 

Very often, the two vertex functions and F** are transformed into the charge and spin vertices 

F^(fc,fc') = 2T\k,k')~T''{k,k'), (H.7) 
F"(fc,fc') = -F^(fc,fc')- (H.8) 

The similar combinations 2F^ — F'^ and — F*^ are also important, but no symbol has been introduced 
for them. 

For every vertex function of two momenta F^(fc, k'), I have introduced the following notation 

F^(/c, fc') = ^ if k e P,, k' e Pj, (H.9) 

where Pi and P2 are the two saddle point patches. Every time this notation is used, it is understood 
that k and k' are (^enera^ momenta in the saddle point patches (i.e. for example 1^ ^121 although 
F^(A:, fc + g) = V^ik.k + Q)). Obviously, F^j = T^^ and F^2 = ^21- The odd and even combinations 
are denoted 

F^=Fti±Ft2. (H.IO) 
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An alternative classification of the same special vertices T^^ 



is obtained as follows. If the mo- 
menta are restricted to the saddle point patches, the vertex function r(A:i, ^2,^3) "decays" into four 
disconnected functions (see Fig. 4.7) 



r(fci,fc2,A;3) = 



Ti{ki,k2,k3) ;ki,k3 e Pi and k2,k4 £ P2 
r2(fci, /c2, fcs) ;ki,k4 e Pi and k2,k3 e P2 
T3{ki,k2,k3) ;ki,k2 € Pi and k3,k4 G P2 
T4{ki,k2,k3) ;ki, . . . ,k4 e Pi 



(H.ll) 



In each of these functions I specify one general vertex and three special vertices which correspond 
to specific configurations of the momenta 



ri(fci,fc2,A:3) 



r2(fci, fc2, fcg) 



^3{ki,k2,h 



T4{ki,k2,k3 



Ti 

r'l' 



r2 
pi? 

1 9 



r'' 
^ 12 

^ 12 
^ 11 



^ 11 

^ 11 

^ 12 



if |k3-k2-Q|, |k3-ki|, |ki+k2-Q| >0(\/A) 

if ki + k2 = Q 



if k3 
if k3 



ki 
k2 





Q 



if |k3-k2|, |k3 

if ki + k2 = Q 
if ks ki = Q 
if ks - k2 = 

if |k3-k2- 



ki -Ql, |ki +k2 -Ql > 0(^/A) 



(H.12) 



(H.13) 



rf c's ^ rf2'^^ if ki + k2 

if ks - ki 
if ks — k2 



^3 — ^ 12 



^ 3 ^ ^ 12 



Ql, |k3-ki-Q|, |ki+k2| >0(VA) 


Q 
Q 



r4 if |k3 

T^CS ^ pBCS if ki - 

if k3- 
if k3- 



-k2| 

k2 = 



^ 4 — ^ 11 
^4 — 11 



|k3 


ki = 
k2 = 



kii, |ki + k2| >o(Va) 



(H.14) 
(H.15) 



The combinations, which are important for the different susceptibilities are the following 



SDW/SF 




'nx -r- T^x 
— 2 + ^ 3i 




CDW/CF 




= 2rf - rf ± (2r^ 


-rf), 


sSC/dSC 


YBCS 


_ pBCS _|_ pSCS 




S/Sd 


-r| 






N/Nd 


2r± - r| 


= 2r{ - r| ± {2tI 




7?s/IId 









{H.16) 
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